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Introduction 



When the electron was discovered in 1897, the first elementary particle had been found. This 
can be considered the birth of particle physics, the science of elementary particles and their 
fundamental interactions. The Standard Model of particle physics (SM) is in agreement with 
nearly all phenomena observed over the last decades, which have partly been measured with 
excellent precision. However, there are also strong arguments that it might need to be embedded 
in a more general theory. 

The Large Hadron Collider (LHC) at CERN, Geneva, is a proton-proton collider, which started 
operation in September 2008. It was built to search for the still undiscovered Higgs boson and 
for physics beyond the Standard Model, but also to perform precision tests of SM processes. In 
2011, protons were brought to collision at a centre-of-mass energy of 7 TeV with instantaneous 
luminosities of up to several 10^^ cm~^s~^. The data used for the analysis presented in this thesis 
were taken with the ATLAS detector, a general purpose detector, which has been designed to 
cope with the uniquely high collision rates, energies and instantaneous luminosities provided by 
the LHC. 

The top quark is the heaviest of the known elementary particles. It has been discovered in 
1995 [l|[2] at the Tevatron, Batavia (Illinois), and since then many of its properties have been 
measured at the Tevatron and at the LHC. Amongst others, the production, the mass, the decay 
and the spin properties of the top quark have been studied. Another fundamental quantity is 
the electromagnetic charge of the top quark, which in the Standard Model is predicted to be +| 
of the proton charge. Recently, an alternative charge hypothesis of — | the proton charge was 
ruled out (3j|6]. 

In the SM, the electromagnetic charge of a particle determines its electromagnetic coupling - 
the interaction with other electrically charged particles by the exchange of photons. Hence, the 
amount of photon radiation from a particle is directly sensitive to its electromagnetic coupling. 

Accordingly, top quark pair (tt) events with additional photons in the final state, denoted 
tt'j events in the following, are sensitive to the electromagnetic coupling of the top quark. The 
investigation of such events is an important test of the SM, because it provides a direct measure 
of the charge of the top quark and of a possible anomalous structure of its electromagnetic 
interaction. 

An important step towards precision tests of the coupling itself is the measurement of the 
production cross section of ttj events, 0"^^^, at the LHC. This thesis presents the first measure- 
ment of (Tff^ • BR at a centre-of-mass energy of 7 TeV, where BR is the branching ratio into tt'j 
decays with one or two leptons in the final state (single lepton and dilepton channel, respectively). 

The measurement was performed in the single lepton channel, which features a large variety of 
particles in the final state and, hence, also of experimental signatures, which need to be identified 
and distinguished: photons, electrons, muons, jets from partons, and missing transverse energy 
from the neutrino, which escapes the detector without interacting. In particular, the presence of 
the photon added complexity to the analysis: an analysis strategy which reduced the dependence 
of the measurement on simulations was set up and the contributions from the most important 
background processes were estimated from data. 
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Dominant background contributions are processes in which a hadron or an electron is misiden- 
tified as a photon. In addition, also processes with real photons in the final state were considered, 
such as the production of W bosons with additional jets and a prompt photon. 

A template fit to the sum of the transverse momenta in a cone around the photon candidate 
was performed, from which the number of tt'j events in the selected data sample and its sta- 
tistical uncertainty were estimated. This was translated into a measurement of a^ty ■ BR, and 
the statistical and systematic uncertainties of the measurement were evaluated. Finally, the 
statistical significance of the result was estimated. 

This thesis is divided into 14 chapters: the theoretical background and the experimental setup 
are introduced in chapters 1 and 2. In chapters 3-6, the modelling of the relevant processes 
and the selection of a data sample enhanced in ttj production are described. In chapter 7, the 
analysis strategy is outlined, which is then detailed in chapters 8-11. The derivation of the 
systematic uncertainties is described in chapter 12 and the final result is presented in chapter 13. 
In chapter 14, the results are summarised and an outlook for further studies is given. 

This analysis was first presented at the HCP conference 2012 in Paris |7|. The results pre- 
sented in this thesis include improved re-evaluations of several background contributions and of 
the systematic uncertainties, which resulted in a lower measured tt^ cross section and a lower 
statistical significance of the ttj signal. 

Natural units are used throughout this thesis {h = c = 1), and hence masses, energies and 
momenta are expressed in units of [GeV]. To avoid confusion, length and time are given in 
usual SI units ( [m] and [s] ) , because they do not refer to subnuclear but to macroscopic detector 
quantities. 

In many histograms, the last (first) bin contains also the sum of all entries above (below) 
the range of the histogram (overflow or underflow bin, respectively). This is noted for each 
histogram in the description. 
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1 tt7 production in the Standard Model 
and beyond 

In the Standard Model of particle physics (SM), the top quark is the weak isospin partner of 
the bottom quark. It is the most massive elementary particle known today. The top quark is 
special in the SM, because its Yukawa coupling of the order of 1 results in a large mass, which 
is much larger than the masses of all other particles. 

The production of top quark pair (tt) events with an additional photon in the final state (^7 
events) is predicted by the SM, and the amount of photons radiated from top quarks is given 
by the electromagnetic coupling of the top quark. Discrepancies in tt'j production with respect 
to the SM prediction would indicate an anomalous structure of the t7-vertex. 

In Sec. a brief summary of the SM is given. More detailed descriptions of the SM can 
be found in various text books, for example in Ref. (8|[9]. The properties of the top quark 
are discussed in more detail in Sec. |1.2[ together with an overview of experimental results. In 



Sec. 1.3, the production of ttj events in the SM and beyond is summarised. 



1.1 Brief summary of the Standard Model 

The SM is a theory of the interactions between elementary particles. It is based on a quantum 
field theory in which interactions are introduced by local gauge symmetries. The SM is very 
successful in describing a large variety of phenomena in particle physics. 



Fig. 1.1 shows a representation of the known elementary particles in the SM. Several properties 
of the particles are shown in an overview in Tab. |l.l| and |1.2[ The fermions - leptons and quarks - 
are spin-^ particles. The interactions between fermions are described by the exchange of spin-1 
gauge bosons: the electromagnetic force is carried by photons (7), the weak force by and Z 
bosons and the strong force by gluons (g). While photons and gluons are thought to be massless, 
and Z bosons are massive. Only particles which are electrically charged interact via the 
electromagnetic force. Weak charge is carried by all fermions, which therefore interact via the 
weak force. Colour charge, however, is only carried by quarks and gluons, which hence interact 
via the strong force. The gravitational force is not described by the SM. 

Leptons and quarks exist in three generations, where the masses of the particles increase from 
generation to generation. Additionally, for every fermion, there is an antiparticle with the same 
properties like the particle, but with opposite values of the additive quantum numbers, such as 
electric charge and the third component of the weak isospin, I3. 

Each quark generation consists of an up- type quark with I3 = +1 and its down- type quark 
partner with I3 = — ^. Due to the parity violating nature of the weak force, only left-handed 
quarks form doublets of up- and down-type quarks, while right-handed quarks form singlets. 
The quark doublets read: up (u) and down (d) quark, charm (c) and strange (s) quark, and 
top (t) and bottom (6) quark. Up-type quarks have an electric charg^of +|e, the charge of 
the down- type quarks is — ^e. 



is the absolute value of the charge of the electron. 
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I II III 

Generation 

Figure 1.1: Known elementary particles in the Standard Model: there are three generations of 
quarks and leptons. The interactions are mediated by four types of gauge bosons. 
The yet undiscovered Higgs boson is not depicted. 



Fermion 


Mass 


Electric 
charge 


h 


Colour 
charge 


up quark 


0.0015 




0.0033 


GeV 


+h 


+ 1 


yes 


down quark 


0.0035 




0.0060 


GeV 


-¥ 


1 

2 


yes 


charm quark 


1.27 


+ 


0.07 
0.11 


GeV 






yes 


strange quark 


0.104 


+ 


0.026 
0.034 


GeV 


-¥ 


2 


yes 


top quark 


173.2 


± 


0.9 


GeV 


+h 




yes 


bottom quark 


4.20 


+ 


0.17 
0.07 


GeV 


-h 


1 

2 


yes 


electron neutrino 




< 


0.000002 


MeV 







no 


electron 


0.510998910 




0.000000013 


MeV 


— e 


2 


no 


muon neutrino 




< 


0.19 


MeV 







no 


muon 


105.658367 




0.000004 


MeV 


— e 


2 


no 


tau neutrino 




< 


18.2 


MeV 







no 


tau 


1776.84 




0.17 


MeV 


— e 


2 


no 



Table 1.1: Overview of the masses 10 



11 



, electric charges, third components of the weak 
isospin, /3, and colour charges of quarks (upper part) and leptons (lower part). 



The lepton doublets consist of a lepton with electric charge — e (electron e, muon fi, tau r) 
and the corresponding neutrino (z^ej '^fj,, i^t), which is electrically neutral. The charged leptons 
form right-handed singlets. Since neutrinos are assumed to be massless in the SM, no right- 
handed neutrino singlets are foreseen. In neutrino oscillation experiments, however, it was 
shown that neutrinos have non- vanishing masses |10|. Although the neutrino masses have not 
yet been measured, the differences in the squares of their masses were measured in the oscillation 
experiments. However, the neutrino masses must be very small compared to the scales present 
in high energy physics experiments and can hence be ignored in this context. 
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1.2 The top quark in the Standard Model 



Gauge boson 


Mass 


photon 
gluon 

boson 
Z boson 


< I-IQ-^^ eV 


80.385 ± 0.015 GeV 
91.1876 ± 0.0021 GeV 



Table 1.2: Overview of the masses 10 ,12 of the gauge bosons 



Mathematically, the SM is formulated as a renormalisable, Lorentz invariant perturbative 
quantum field theory. Interactions are introduced by local gauge symmetries. The structure of 
the gauge groups is SUc(3) x SUl(2) x Uy(1), where SUc(3) is the gauge group for Quantum 
Chromodynamics (QCD) [l3[|l5] , which describes the strong interaction, and SUl(2) x Uy(1) is 
the gauge group for the unified electromagnetic and weak interactions [l6| - |18| . At the electroweak 
scale of 246 GeV, the symmetry between electromagnetic and weak interactions is spontaneously 

This mechanism describes the generation of particle 

to exist. 



masses in the SM, but also requires the yet undiscovered Higgs boson [22 23 

After breaking of the electroweak symmetry, the Lagrangian of the electromagnetic interaction 
only is described by Quantum Electrodynamics (QED): 



broken via the Higgs mechanism 19 -21 



1 



QED 



(1.1) 



The first two terms in Eq. ( |1.1[ ) are the kinetic terms for the fermion field ip and the photon 
field Afj^, respectively. The term (^—q'lp^'^Afj^il:^ describes the interaction between fermions and 
photons, which is determined by the electric charge q of the fermion. 

With ti^ events, the interaction between top quarks and photons can be studied. Hence, a 
measurement of ti^ events can be interpreted as a measurement of the charge of the top quark 
in the SM, or as a measurement of the structure of the interaction term. This is discussed in 
more detail in Sec. 11.3.21 

Although current data are in agreement with the SM, it is widely believed to be just an effective 
theory at low energies, because several questions remain unanswered in the SM. Examples are the 
hierarchy problem involving the fine-tuning of large corrections to the mass of the Higgs boson, 
missing explanations for dark matter and dark energy, the large number of free parameters in the 
SM, the unification of the electroweak and the strong interactions, and the inclusion of gravity 
in the model. 

In case the SM is just an effective theory, it would need to be embedded in a more general 
theory and measurements would differ from SM predictions starting at a certain energy scale. 
Several models, such as supersymmetry |24|, models with extra dimensions 25 ,26 or techni- 



color 127 -30 provide solutions to some of the issues of the SM and predict new phenomena to 



appear at a scale of the order of 1 TeV. 



1.2 The top quark in the Standard Model 

In this section, only some aspects of top quark physics can be highlighted. Detailed information 
on top quark physics at hadron colliders can be found in Ref. (sTl. The current status of top 



quark physics is summarised in Rcf. [10 
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1.2.1 Top quark production in pp collisions 

In proton-proton (pp) collisions, top quarks are dominantly produced in pairs via the strong 
interaction, but also the production of single top quarks via the weak interaction is possible. 



The contributing diagrams for production in leading order (LO) are shown in Fig. 1.2 They 
can be categorised into s-channel quark-antiquark annihilation (upper diagram) and gluon-gluon 
fusion in the t-, u- and s-channels (lower diagrams). 




Figure 1.2: Diagrams for ti production via the strong interaction in leading order: the upper 
plot shows the quark-antiquark annihilation, the lower plots show the different 
contributing diagrams from gluon-gluon fusion. 



In contrast to e^e~ colliders, at high-energy hadron colliders the interacting particles are par- 
tons, which are confined in the colliding hadrons. Hence, the parton-parton cross section a^^^^^ 
for partons i and j needs to be convolved with the parton distribution functions (PDFs) [32-34 



of the hadron, fi{xi, fi'jp), which describe the probability to find a parton i within the hadron with 
a fraction Xi of its momentum. The PDFs are evaluated at a factorisation scale fip, which sepa- 
rates perturbative QCD from non-perturbative effects. Hence, in pp collisions, the tt production 



cross section reads 31,35 



crPP^"(\/s,mt) = ^ / dxidxj/i(xi,//|)/j(xj,//|) • <T*^"'**(mt,\/5,a;i,Xj-,as(/u|),//|) , 
i,j=g,q,9 

- (1-2) 

where rrit is the top quark mass, ^/s is the centre-of-mass energy of the pp collision, is the 
centre-of-mass energy of the parton-parton system and /i/j is the renormalisation scale, which is 
introduced to allow for finite-order calculations in QCD perturbation theory. A typical choice 
for fj,R as well as for fip is the energy scale for the process under study. In processes with top 
quarks, the top quark mass is often chosen as this energy scale. While ij*-?^** can be calculated 
in perturbative QCD, the PDFs can not and need to be extracted from data measured in ep, pp 



and pp collisions and from other sources 36 - 391 . 



An approximate next-to-next-to-leading order (NNLO) calculation for (j^p^** at y/s = 7 TeV 



40 



41 



_ig pb using HATHOR 142] and CTEQ PDFs 



36 



at a top 



yielded a prediction of 165 
mass of 172.5 GeV, where the uncertainty is due to the uncertainties on the PDFs, and on the 
factorisation and renormalisation scales. 
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1.2 The top quark in the Standard Model 



1.2.2 Top quark decay 

In the SM, about 99.8% of the top quarks decay into a boson and a 6-quark with a hfetime 
of about 0.5 • 10-24 s llOj. Decays to a W+ boson and a d- or s-quark are strongly suppressed 



by the CKM matrix elements Vtd = 0.00359 ± 0.00016 and Vts = OMWto fon \}^- Since 
the lifetime of top quarks is shorter than the typical time for hadronisation, they do not form 
bound states before they decay. Hence, the top quark is the only known quark for which the 



quark properties are accessible, such as for example its electromagnetic coupling (cf. Sec. 1.2.5 



and 1.3). 

The boson from the decay of the top quark can either decay into a pair of up-type quark 
{u or c) and down-type antiquark (d, s or b), or into a charged lepton and the corresponding 
neutrino. This is illustrated in Fig. |1.3[ The dominant decay of the antitop quark analogously 
reads t — )■ W~b. 




Figure 1.3: Top quark decay: almost all top quarks decay into a bottom quark and a W'^ 
boson. The boson can then either decay into a charged lepton and the cor- 
responding neutrino (left), or into an up- type quark and a down-type antiquark 
(right). 



1.2.3 Experimental signatures of tt production 

Since W bosons can decay into quarks (hadronic decay) or leptons (leptonic decay), the decay 
of top quark pairs can be categorised into three channels, depending on the number of leptons 
and quarks in the final state. Since quarks hadronise and form jets, bundles of particles, the 
categories are labelled accordingly: the all-hadronic (or alljets) channel, where both W bosons 
decay hadronically, the dilepton channel, where both W bosons decay leptonically, and the 
single lepton (or lepton+jets) channel with one W boson decaying leptonically and the other 
one decaying hadronically. The branching ratios of the respective channels are illustrated in 
Fig.fLil 

While the all-hadronic channel has the largest branching ratio, at hadron colliders it suffers 
from large background contributions from multijet production. From an experimental point of 
view, electrons and muons provide a clear signature, while the identification of r-leptons is more 
involved, because of the different leptonic and hadronic decay modes of the r-lepton. 

The single lepton channel is sometimes called the golden channel, because of its clear signature 
and relatively high branching ratio (about 30%). It is characterised by a high-energetic electron 
or muon, four high-energetic jets, out of which two are originating from 6-quarks, and a large 
imbalance of the momentum in the transverse plane due to the neutrino, which escapes detection. 
Depending on the lepton type, the final state is called single electron or single muon channel. 
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Since electrons and muons from r-decays are experimentally indistinguishable from electrons and 
muons from the direct decay of the W boson, in the analysis presented in this thesis, r — )• eVeVT 
and T — )• [iv^iVt decays were included in the electron and the muon channel, respectively. 

Top Pair Branching Fractions 




"dileptons" 



e+jets 15% 



alljets" 46% 



T+jets 15% 



jets 1 5% 



'iepton+jets' 



Figure I.4: Experimental signatures of top quark pair decays and their branching ratios 43 
in about 30% of the cases, top quark pairs decay into either the single electron or 
the single muon channel. 



1.2.4 Background processes for single lepton decays of tt events 

The background processes for the single electron and single muon channels can be divided into 
reducible and irreducible background processes. Irreducible background processes feature the 
same final state as the signal. This is the case for the production of leptonically decaying W 
bosons in association with jets (VF+jets). These jets arise from QCD corrections to the W boson 
production and can originate from gluons and light quarks, but also from 6-quarks. Moreover, 
events with electroweakly produced single top quarks may feature additional jets from QCD 
corrections and may be misreconstructed as tt events. 

Reducible background processes have a different experimental final state, but one or more par- 
ticles are not correctly reconstructed or are just outside of the detector acceptance. Momentum 
imbalance may always originate from miscalibrations of jet energies. Hence, Z+jets events with 
Z — 7- ee or Z — 7- fifi, where one lepton is not reconstructed correctly, may be misreconstructed as 
tt events. The same holds true for WW, WZ and ZZ production with additional jets from ini- 
tial or final state radiation (ISR/FSR). However, WW — t- li'qq' and WZ — )• li/qq decays are part 
of the irreducible background, strictly speaking. Finally, jets in multijet events may give rise to 
signatures similar to electrons or muons from W boson decays: hadrons from jet fragmentation 
may be misidentified as electrons, and electrons and muons within jets coming from S-meson 
decays may be misidentified as isolated electrons or muons from a W decay, respectively. Due to 
the large cross section for multijet events at hadron colliders, such events are expected to have 
a non-negligible background contribution, even if the misidentification rates are small. 
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1.2 The top quark in the Standard Model 



1.2.5 Top quark properties and experimental results 



In addition to the production cross section for tt and single top events in pp 144 - 50 and pp 51 - 54 



collisions, various properties of the top quark have been studied. In particular, the mass of the 
top quark has been measured with high precision to 173.2 it 0.9 GeV at the Tevatron [II] . 
Measurements at the LHC are in agreement with this result |55[|56| . 

The knowledge about the couplings of the top quark is still limited, although the SM provides 
clear predictions for them. The structure of the Wtb-vertex has been studied in the parity 

and the strength of the coupling has been tested 
To date, results are consistent with SM 



59 



53,54 



violating weak decay of the top quark 57 
in the production of single top quarksj49 
expectations. 

Since the top quark decays before it can form bound states, the spin of the two top quarks from 
ti production is transferred to the decay products. A certain correlation between the top quark 
spins is expected in the SM, which may be altered by additional production processes. To date, 
measurements of the spin correlation at y/s = 7 TeV are consistent with SM expectations 60 - 62 



Also the measurement of the charge asymmetry in tt production is a test of the production 
process. The asymmetry originates from interference effects in NLO between initial and final 
state gluon radiation as well as from interferences between Born and box diagrams. It is small in 
the SM, but significantly enhanced in several alternative models 63 . Measurements performed 



by the CDF and D0 collaborations show discrepancies of up to more than 3a from SM expec- 
tations (64-66 . Since in pp collisions the asymmetry is less pronounced than in pp collisions. 



ATLAS and CMS have not yet gained the sensitivity of the Tevatron experiments. To date, 
measurements of the asymmetry in pp collisions are consistent with SM expectations (67||68] . 

Alternative top quark decays have been searched for and exclusion limits have been set for 
example on flavour changing neutral current decays, such as t — >• qZ 69-72 and t ^ [73] . 
Also, the branching ratio t — )• Wb 74 -76 has been measured. To date, all results are consistent 
with the expectations from the SM. 



The measurement of the number of jets produced in association with tt events 77 is a first 
measurement towards tests of the strong coupling of the top quark. 

In the SM, the electromagnetic coupling is given by the electric charge of the top quark, 
which is predicted to be +|e. The pair production of exotic top quarks with an electric charge 
of - ■ 



78 



79 



results in the same final state as SM tt production. Analyses in which the top 
quark charge is measured from the charge of the lepton and the associated 5-quark could rule 
out the exotic top quark scenario with up to 99% confidence level (3||6j. 

ttj events directly probe the electromagnetic coupling of the top quark without necessarily 
assuming the structure of the t7-vertex predicted in the SM. Hence, they provide a test of the 
SM which is complementary to analyses considering only the exotic top quark scenario. In pp 
collisions, a first measurement of the ttj cross section and the ratio of the cross sections for 
tt7 and tt production at ^/s = 1.96 TeV has been performed by the CDF collaboration 80 . 
The ttj cross section was measured to 0.18 it 0.08 pb, and the signal significance was estimated 
to 3.0 standard deviations. The ratio of the tt^ and tt cross sections has been measured to 
0.024 lb 0.009, consistent with SM expectations, ti'j production is discussed in more detail in 
Sec.Ol 

A variety of additional measurements of the properties of the top quark itself, of its decay 
and of processes with top quarks in the final state has been performed. These measurements 
are not listed in this thesis, but an overview can be found in Ref. [lOj. 

Due to its large mass and its large Yukawa coupling, the top quark may play a special role 
in electroweak symmetry breaking and for the discovery of phenomena beyond the SM. For 
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instance, new heavy particles could manifest themselves in modifications of the SM predictions 
for top quark properties, for example in the production process, which would then be observed 
in measurements of the charge asymmetry. Moreover, since top quarks decay before they form 
bound hadronic states, the top quark is a unique opportunity to study the properties of a bare 
quark, such as the couplings or the spin properties. In particular, top quarks are the only quarks 
for which the electromagnetic coupling is directly accessible. 



1.3 tt'y production in pp collisions 



1.3.1 Production of ttj events in the Standard Model 



As outlined in Sec. 1.2.5 the production of ttj events is sensitive to the electromagnetic coupling 
of the top quark. In the SM, it is solely given by the electric charge of the top quark - cf. 



Eq. (1.1). Hence, within the SM, a measurement of the tt^ cross section can be interpreted as 



a measurement of the top quark charge. 

In ti events, photons can be radiated from all charged particles, including the top quark, 
but also from incoming quarks and the charged decay products of the top quarks. Artificially, 
the contributing diagrams can be divided into diagrams from radiative production and radiative 
decay, but the tt'y final state is only well-defined when all interference terms between the different 
processes are taken into account. 



In LO, radiative tt'y production can occur in quark- antiquark annihilation (Fig. 1.5) or gluon 



gluon fusion (Fig. 1.6) |81| - similarly to tt production (Sec. 1.2.1). Photons can be radiated 
from the incoming quarks or from the top quarks. The LO processes for the radiative decay are 



sketched in Fig. 1.7 After a top quark pair is produced, a photon can be radiated either from 
the decaying top quark, the W boson or the 6-quark. Photons can also be radiated from the 
charged decay products of the W boson, the charged lepton or the quarks, respectively, which 



is not illustrated in Fig. 1.7 



It has been shown that interferences between the different diagrams are not negligible 82 and 



hence, a separation between processes with photons radiated from top quarks and processes with 



photons radiated from other particles is not well defined. From Eq. (1.1) it could be concluded 



that the cross section for tt^ production is proportional to the square of the top quark charge. 
Since also diagrams with photons radiated from other particles need to be added to the definition 
of the inclusive ttj cross section, the dependence on the top quark charge is not just simply 
quadratic. 

A calculation of the inclusive cross section in pp collisions at ^/s = 14 TeV in next-to-leading 
order (NLO) QCD is available 83 . In this calculation, infrared divergencies were avoided by a 



minimum photon transverse momentum of 20 GeV. In order to get rid of collinear divergencies. 



photons were required to fulfil a Frixione-type parton isolation 84 with a width parameter 



of 0.4. The resulting cross section in NLO was found to be larger than the calculation in LO 
and, hence, the ratio of the two calculations, the so-called NLO k-factor, differs from unity. 

For ^/s = 7 TeV, this fc-factor was estimated to 2.6 it 0.5 |j85|. For this calculation, the cut 
on the photon transverse momentum was adjusted to 8 GeV in order to be applicable to the 



simulations used in the measurement presented in this thesis (Sec. 4.1). The effect of different 
isolation criteria used in the calculation and the simulation was expected to be well covered 
by the uncertainty on the /c-factor, which was estimated by changing the renormalisation and 
factorisation scales from the nominal value of 2mt to rut- 

The /c-factor was applied to the ti'j LO calculation obtained with the WHIZARD Monte Carlo 
generator |86l|87|, which is described in more detail in Sec. 4.1 For ^/s = 7 TeV an estimate 
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Figure 1.6: Leading order diagrams for tt'^ production in gluon-gluon fusion. 
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7 




Figure 1.7: Top quark decay with additional photon radiation: photons can be radiated from 
the top quark, the bottom quark or the W boson. 



of cr^f^ • BR = 2.1 lb 0.4 pb was obtained, where BR is the branching ratio into the single lepton 



and dilepton decay modes (see Sec. 4.1). 



1.3.2 Production of ttj events beyond the Standard Model 

In Sec. |1.2.5[ it has been mentioned that the top quark may be a window to physics beyond the 
SM. In particular, it may be of interest for an understanding of electroweak symmetry breaking, 
because its mass is of the same order as the electroweak scale. In order to observe deviations 
from SM predictions, precision tests of the properties of the top quark are necessary, where 
tests in the electroweak sector are especially interesting. In ti'j events, deviations from the SM 
electromagnetic t7-vertex = —iQiCy^ can be searched for. 

The following generalised form of the t7- vertex is a natural extension of the SM |88l|89|: 



7- il') + Ff [q^] 75) + {Fr id') + Fi {q^) 75) 



pi/ 



where q is the photon four-momentum, Qt is the value of the top quark charge in units of e, and 
mt is the mass of the top quark. and F^ are the form factors for vector and axial-vector 
couplings, respectively. The form factors F^ and F^ represent magnetic and electric dipole 
moments of the top quark. In the SM, all form factors vanish at tree level except for F^^ , which 
is equal to unity. F^ and receive non-zero contributions only when higher loop corrections 
are considered. Deviations from these predictions, for example enhanced dipole moments, would 
indicate the presence of phenomena beyond the SM. 

At hadron colliders, a precise determination of the form factors will be challenging due to 
the limited precisions in the measurement of the four-momenta and in the identification of 
the particles involved. However, with growing statistics at the LHC, the sensitivity to the 
electromagnetic form factors will increase. Previous studies in the ti^ topology [81,82 have 



focused on the discrimination of two different scenarios for the electromagnetic charge of the top 
quark {+ie and — le, cf. Sec. 1.2.5). However, in Ref. \82i it has been estimated that the charge 



of the top quark can be measured with a precision of about 10% with an integrated luminosity 
of 10 fb"^ at = 14 TeV. 

Electron-positron colliders would provide a cleaner environment for studying the electromag- 
netic coupling of the top quark. For an integrated luminosity of 10 fb~^ at y/s = 500 GeV, a 

Even better limits would 



precision of 5 - 10% on the axial form factors has been predicted |90| 
be achieved at future photon-photon colliders p3l|. 
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The ATLAS detector is one of the main experiments located at the Large Hadron CoUider (LHC) 
at CERN. It has been built to cover a large physics program with collisions of protons and heavy 
ions at unprecedented energies and extremely high rates and instantaneous luminosities. 

In Sec. 2.1, the accelerator is introduced. In Sec. 2.2 the detector and its different subcom- 
ponents as well as the detector readout are described. 



2.1 The Large Hadron Collider and its experiments 



The LHC [92] is a circular pp collider at CERN near Geneva, Switzerland. In addition to 
protons, also heavy ions can be brought to collision. The accelerator is located in the tunnel of 
the dismounted Large Electron-Positron Collider, which has a circumference of 27 km. Fig. |2.1| 
shows a sketch of the LHC and its preaccelerators. Protons from the ionisation of hydrogen 
atoms are accelerated to 50 MeV in a linear collider (LINAC2) before entering the Proton 
Synchrotron Booster (BOOSTER). In the Proton Synchrotron (PS), the protons gain a total 
energy of 25 GeV. Finally, they enter the Super Proton Synchrotron (SPS), in which they gain 
an energy of 450 GeV, which is the nominal injection energy for the LHC. 




BDDSTER 



LINAC2 ». 




Figure 2.1: Sketch of the Large Hadron Collider, its preaccelerators and the four main exper- 



iments at the interaction points (not to scale) 94 
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Bunches of protons are injected into separate beam pipes in both directions around the LHC 
ring and are accelerated with a system of radio cavities to an energy of 3.5 TeV, corresponding 
to ^/s = 7 TeV (2010 - 2011). This is roughly 3.5 times more than the highest energies achieved 
at any other particle collider before. Superconducting NbTi magnets with a field of up to 8.33 T 
keep the protons on their circular track. In 2012, the proton energy has been increased to 
4 TeV, corresponding to ^/s = 8 TeV. In the following years, the design centre-of-mass energy 
of ^/s = 14 TeV is aimed for. 

Each proton bunch consists of about 10^^ particles. With a minimal bunch spacing of 25 ns 
and a maximum of 2808 bunches, a pp design luminosity of 10^^ cm~^s~^ can be achieved, which 
results in an average number of more than 20 inelastic collisions per bunch crossing depending 
on the beam focusing. 

In 2011, the largest part of the data was taken with a bunch spacing of 50 ns |95|. With 



this setting, up to 1380 bunches were circulated. Luminosities of up to 3.65 
achieved with up to 1.5 • 10^^ protons per bunch. 



10^3 



cm 



were 



There are four main experiments at the LHC, which are also depicted in Fig. |2.1 

CMS 



which is discussed in more detail in Sec. 2.2 



98 



ATLAS, 



are located 



1 96 1, ALICE [97] and LHCb 
at the four interaction points of the two beams. ATLAS and CMS have been designed to cover 
a wide variety of physics in pp collisions, ranging from precision tests of the SM to searches for 
the Higgs boson and searches for new phenomena. The design choices for the ALICE detector 
were guided by the needs of a detector for ion-ion collisions with strong emphasis on the tracking 
detectors rather than on calorimetry. LHCb is a detector optimised for the study of 6-physics, 
in particular measurements of CP violation in the 6-sector. Since bb systems are highly boosted 
at the LHC, LHCb has an asymmetric design to focus on the forward region on one side of the 
pp interaction point. 



44m 



25m 




Tile calorimeters 



LAr iiadronic end-cap and 
forward calorimeters 



Toroid magnets 
Muon ctiambers 



Pixel detector 
LAr electromagnetic calorimeters 
Solenoid magnet | Transition radiation tractser 
Semiconductor tracl<er 



Figure 2.2: Drawing of the ATLAS detector showing the different detectors and magnet sys- 



tems 99 
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2.2 The ATLAS experiment 



Fig. 2.2 shows a drawing of the ATLAS detector [99, 100 with its solenoid and toroid magnets 
and the different subdetector systems. The subdetectors - inner tracking detector, calorimeters 
and muon spectrometer - are briefly described in Sec. 2.2.1 - 2.2.3 The magnet, trigger and 



data acquisition systems are introduced in Sec. |2.2.4| and |2.2.5[ respectively. 

All detector systems have been designed to cope with two main challenges set by the LHC: 
on the one hand, the high event rate puts special requirements to the detector. Fast and 
radiation-hard electronics and sensor elements are necessary as well as an efficient trigger and 
data acquisition system. On the other hand, very good particle identification is the key to an 
efficient suppression of the large background contribution from multijet production at a hadron 
collider, as well as from additional inelastic interactions from the same bunch crossing (pile- up). 



The design of the ATLAS detector follows the structure illustrated in Fig. 2.3, which shows 



schematically the interactions of different types of particles with the detector material: the 
innermost detector layer is a tracking detector, in which electrically charged particles are traced 
for example by creating electron-hole pairs in semiconductors or by ionising gas. Typically, the 
tracking detector is embedded in a magnetic field so that the momenta of the particles can be 
measured from the curvature of the tracks. 



tracking electromagnetic hadronic muon 
chamber calorimeter calorimeter chamber 



photons 
electrons 
muons 
pions 




innermost layer 



outermost layer 



Figure 2. 3: Schematic view of the interactions of different types of particles in a general pur- 
pose detector at a modern high energy particle collider. 



The next detector layers consist of electromagnetic (EM) and hadronic calorimeters. The 
calorimeters are massive and therefore induce electromagnetic showers from electrons and pho- 
tons by pair production and bremsstrahlung, and hadronic showers by various processes between 
hadrons and matter. Electron and photon showers are typically contained in the electromagnetic 
calorimeter while hadronic showers range into the hadronic calorimeter. 

Since almost all muons produced at the LHC are minimum ionising particles, they are the only 
electrically charged particles which may pass the calorimeters and reach the outermost layer, 
which is made of tracking chambers. Hence, these detectors are called muon chambers. 

In principle, it is desirable to cover the full solid angle with sensitive detector material. How- 
ever, in the very forward region some space must be kept open for the beam pipes. Moreover, 
support structures, cables, cooling systems etc. need to be included in the detector design and 
reduce the sensitive volume. 

2.2.1 Inner detectors 

The ATLAS Inner Detector (ID) consists of three subdetector systems: the Pixel detector 
and the Semiconductor Tracker (SCT), which use silicon semiconductor technology, and the 
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2.1m 




Figure 2-4: Overview of the Inner Detector 99 



the left figure shows a longitudinal section 
of the Inner Detector with the different subdetectors. The right figure shows a 
transverse section and illustrates the distances of the different detector layers from 
the beam line. 



Transition Radiation Tracker (TRT), which exploits the transition radiation produced in a gas 
mixture of Xe, CO2 and 02- Fig. |2.4| shows a longitudinal and a transverse section of the ID. In 
particular, the distances of the different subdetector layers from the beam line are illustrated. 
The whole ID is embedded in a 2 T solenoidal field (Sec. 2.2.4). 

With three concentric cylinders (barrel part) and three endcap disks, perpendicular to the 
beam axis, the Pixel detector covers a rang^of \r]\ < 2.5. Each of the 1744 sensors consists of a 
segmented silicon wafer with pixels of minimum area 50 x 400 /im^ and 46080 readout channels. 
The innermost pixel layer, the so-called 6-layer, is as close to the beam line as 50.5 mm and 
allows for a precise extrapolation of tracks to the vertices. This is crucial for any 6-tagging 
strategy based on impact parameters and the identification of secondary vertices. 

The SCT consists of four layers in the barrel and nine endcap disks. It covers the range 
\7]\ < 2.5. The sensors use silicon microstrip technology with a strip pitch of 80 fim. In the 
barrel, the strips are arranged parallel to the beam line, while in the disks, the strips are oriented 
radially. Modules are arranged back-to-back with a small stereo angle of 40 mrad to allow for 
a measurement of the azimuth angle in each layer. A typical track yields three space-points in 
the Pixel detector and eight in the SCT. Together, the silicon trackers ensure the measurement 
of the track momenta and the identification of primary and secondary vertices. 

In the barrel part of the TRT, there are 73 planes of straw tubes filled with gas, which 
are arranged parallel to the beam axis. In the endcap, there are 160 straw planes, oriented 
radially. The TRT covers a range of \rj\ < 2.0, in which the separation of electrons from charged 
pions is improved by exploiting transition radiation. Although the TRT does not provide track 
information in the direction along the beam line, pattern recognition and the measurement of 
the track momenta become more robust by using the signals from the TRT. 



ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre 
of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC 
ring, and the y-axis points upwards. Cylindrical coordinates (r, (j>) are used in the transverse plane, where (j) 
is the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the polar angle 9 as 
r? = -ln [tan(f)]. 
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The total amount of material of the ID is as large as roughly 0.5 electromagnetic radiation 
length^ Xq for \rj\ < 0.6. For 0.6 < |r/| < 1.37 and 1.52 < \rj\ < 2.5, the amount of material 
reaches up to 1.5Xq. In the barrel-to-endcap transition region at 1.37 < \ri\ < 1.52, the amount 
of material is even larger. Electrons and photons in this region were not taken into account 
in this analysis. A particular consequence of the sizable amount of material in front of the 
calorimeters is that a large fraction of photons convert into e~^e~ pairs in the ID volume. 



2.2.2 Calorimeters 



Fig. 2.5 shows an overview of the different electromagnetic and hadronic calorimeters of the 
ATLAS detector. All calorimeters are sampling calorimeters consisting of alternating layers 
of dense absorber material and active material, where only the active material is used for the 
energy measurement. This design allows for a compact size of the calorimeter system. 

The hadronic calorimeter in the barrel (Tile) uses steel as absorber and scintillators as ac- 
tive material. All other calorimeters use liquid argon (LAr) technology with different types of 
absorbers: lead in the electromagnetic barrel (EMB) and the electromagnetic endcap calorime- 
ter (EMEC), copper in the hadronic endcap calorimeter (HEC) and the electromagnetic part 
of the Forward Calorimeter (FCal), and tungsten in the hadronic part of the FCal. The LAr 
calorimeters are placed in three cryostats: one for the barrel and one for each endcap. 




The technologies have been chosen to provide fast readout, radiation hardness and high con- 
tainment of electromagnetic and hadronic showers to ensure a precise measurement of their 
energies. The energy flux is varying in the different detector regions. Especially in the very for- 
ward region, which is covered by the FCal (3.1 < |r/| < 4.9), extremely high fluxes from minimum 
bias events drove the design towards dense absorber material and small LAr gaps. 

All calorimeters are finely granulated and also segmented longitudinally to allow for a precise 
determination of the position of the showers and to distinguish different shower types by the 
use of shower shapes. This is particularly important for the central region, which is devoted to 

^The radiation length is defined as the typical amount of material traversed by an electron after which it has 
lost i of its original energy by bremsstrahlung. 
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precision measurements of electrons and photons: the EMB {\rj\ < 1.475) is segmented into three 
longitudinal layers, where the first layer, the so-called LAr strips, provide a very fine granularity 
in 7] of 0.0031. To ensure continuous coverage in azimuth and to enable fast readout, the lead 
absorbers are folded into an accordion shaped structure. 

A similar design as for the EMB has been used for the EMEC, which is divided into two 
wheels covering the ranges 1.375 < \r]\ < 2.5 and 2.5 < \ri\ < 3.2. The inner wheel has a coarser 
granularity in ij and (p, limiting the region devoted to precision physics to |r/| < 2.5. A thin 
LAr layer (presampler) is placed in front of the EMB and the EMEC for It^I < 1.8 to correct for 
energy lost in front of the calorimeter. 

The Tile calorimeter is located behind the EMB and the EMEC and is divided into three 
longitudinal layers. It consists of a central barrel (|r/| < 1.0) and an extended-barrel part 
(0.8 < \r]\ < 1.7). The radial depth is about 7.4 nuclear interaction length^ i^^i)- 

The HEC is a traditional LAr sampling calorimeter covering the region (1.5 < |?7| < 3.2), which 
is placed behind the EMEC in the same cryostat. It consists of two independent wheels, each 
of which is divided longitudinally into two parts. 

Altogether, the calorimeters cover the range \r]\ < 4.9 and, thus, provide good hermiticity to 
ensure also a precise measurement of the imbalance of the transverse momentum. Over the 
whole range in r], the total thickness of the calorimeter system is approximately lOA/, ensuring 
good containment of electromagnetic and hadronic showers and limiting punch-through effects 
to the muon spectrometer. 



2.2.3 Muon detectors 

The ATLAS muon system covers the range \r]\ < 2.7 and is designed to measure the momenta 
of muons exiting the calorimeter system starting at energies above ~ 3 GeV. The tracks of the 
muons are bent by the magnetic field of the air-core toroid system in the barrel and in the 



endcaps (Sec. 2.2.4). The fields in the barrel and in the endcaps are oriented such that muon 



tracks in both regions are mostly orthogonal to the field lines. 



An overview of the different subsystems is shown in Fig. 2.6: the muon system consists of 
high- precision tracking chambers as well as trigger systems. In the barrel part. Monitored 
Drift Tubes (MDTs) are used for tracking and Resistive Plate Chambers (RPCs) for triggering 



(Sec. 2.2.5 ). In the endcaps, tracking information is provided by Cathode Strip Chambers (CSCs) 
and Thin Gap Chambers (TGCs) are used for triggering. In the barrel as well as in the endcaps, 
muons typically cross three longitudinal layers of the muon spectrometer. The muon system is 
divided into eight octants with overlaps in (p to avoid gaps in the detector coverage. 

The technologies for the tracking systems have been chosen such that high precision can be 
achieved given the level of the particle fiux. The MDTs in the barrel part follow a robust and 
reliable detector design. Since each tube contains only one sense wire, the simple geometry 
allows for the prediction of deformations as well as for a precise reconstruction. As the particle 
flux is increasing with \rj\, the CSCs are more suited for the endcap region: the higher granularity 
of the multiwire proportional chambers facilitates to cope with the increasing rates. 

The choice of the technologies for the trigger chambers was driven by the requirement for 
fast and highly efficient trigger capabilities given the different conditions present in the barrel 
and endcap regions during data taking. Additionally, an adequate resolution of the transverse 
momentum of the tracks was required. In the barrel, RPCs provide good spatial and time 
resolution. In the endcap regions, however, higher particle fiuxes as well as the need for a higher 



^The nuclear interaction length A/ for hadrons is defined in analogy to the electromagnetic radiation length Xo 
for electrons and photons. 
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granularity required a different technology: TGCs are used for the region 1.05 < |r/| < 2.4. They 
are based on the same principle as multiwire proportional chambers and fulfil the needs in terms 
of rate capability and granularity. With RPCs and TGCs, a time resolution of 15 — 25 ns can be 
achieved, which is sufficient for fast trigger decisions and a good association of tracks to bunch 
crossings. 

The benchmark for the tracking performance of the muon spectrometer is set by a 10% 



resolution on the transverse momentum of 1 TeV muons 99 . To achieve this goal, the position 
of the MDT wires and the CSC strips must be known with a precision better than 30 ^m. 
Therefore, a high- precision optical alignment system was set up to monitor the relative position 
of the MDT chambers and their internal deformations. 



2.2.4 Magnet system 

The ATLAS magnetic system consists of four large superconducting magnets: a central solenoid 
and three toroid magnets in the barrel and the two endcaps. A sketch of the magnet system is 



shown in Fig. 2.7 The solenoid and the toroids are shown, as well as the tile calorimeter. 



The central solenoid provides an axial field with a strength of 2 T. The magnetic flux is 
returned by the tile calorimeter and its girder structure. The solenoid was designed to be 
particularly lightweight and to minimise the amount of material in front of the calorimeter 
system to which it contributes only a total of 0.66 electromagnetic radiation lengths. 

The toroid systems provide magnetic fields with a bending power of 1.5 — 5.5 Tm in the 
barrel and 1 — 7.5 Tm in the endcap regions. Each system consists of eight air-core coils placed 
in aluminium housings. The toroidal fields contain non-uniformities which need to be known to 
high precision to allow for an accurate measurement of muon momenta. Hence, 1800 Hall sensors 
were installed in the muon spectrometer volume to enable the monitoring of the magnetic field. 
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Figure 2.7: Sketch of the magnet sys- 

the solenoid and the 



tern 99 



toroids are shown, as well as 
the tile calorimeter. 
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Figure 2.8: Sketch of the trigger chain 
and the different trigger lev- 
els [Too] 



The indicated rates 



are orders of magnitude only. 



2.2.5 Trigger and data acquisition 

Assuming a bunch spacing of 25 ns and approximately 20 inelastic interactions per bunch cross- 
ing the event rate at the ATLAS detector is of the order of 1 GHz. A three-level trigger 
system was set up to reduce this rate to about 200 Hz. In 2011 data taking, the real trigger 



rate was indeed of the order of 300 Hz 95 . The triggers need to suppress minimum bias events 
very strongly while efficiently selecting rare physics events. The data acquisition system (DAQ) 
gathers the data from the different detector subsystems and buffers them until a trigger decision 
is received. When the event is not rejected by one of the trigger levels, the data are recorded per- 



manently. Fig. 2.8 shows a sketch of the ATLAS trigger chain indicating the order of magnitude 
of the trigger rates at the different trigger levels. 

The first trigger level (LI) is a hardware-based trigger, which reduces the event rate to approx- 
imately 75 kHz. Muons, electrons, photons, jets and hadronically decaying r-leptons with high 
transverse momenta are searched for as well as a large momentum imbalance in the transverse 
plane and a large total transverse energy. The muon trigger chambers are used as well as the 
calorimeter system with reduced granularity. Within less than 2.5 //s, Regions-of-Interest (Rol) 
are identified in ?7-(/)-space, which serve as seeds for the decision at the second trigger level (L2). 

The high level trigger is composed of the L2 and the Event Filter (EF), both of which are 
software-based trigger systems. At L2, the energy and direction of the Rols are further investi- 
gated and also the types of the trigger objects are analysed. Within 40 ms a decision is made, 
and the event rate is reduced to below 3.5 kHz. The EF further decreases the rate down to 
roughly 200 Hz. Events passing the EF are stored permanently. The full event information is 
available at the EF level and, hence, energies and directions of the trigger objects are estimated 
with higher precision than at LI and L2. In particular, the discrimination between the different 
particle types is enhanced by the use of the ID tracking system and calorimeter shower shapes. 

Selections of different trigger signatures are collected in so-called trigger menus. For triggers 
with very high rates, only a fraction of the triggered events can be selected on a random basis in 



As mentioned in Sec. 2.1 in 2011 most of the data was taken with a bunch-spacing of 50 ns. 
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order to perform cross-checks and studies of less rare physics processes. The trigger menus are 
adjusted to the data taking conditions, in particular to the instantaneous luminosity, in order 
to make optimal use of the band width available for storage. 

2.2.6 Performance in pp collisions 

The performance of the different subsystems of the ATLAS detector was studied with the data 
taken in 2010 and 2011. Of particular importance is the agreement with Monte Carlo (MC) 
simulations (Ch. |4]). 

The left plot of Fig. 2.9 shows the invariant mass of pairs around the mass of the Z 

boson in 0.70 fb^^ of data, where the muon momentum was measured using ID tracks only. An 
early and an improved version of the alignment of the ID subdetectors were used to measure 
the distribution (full and open circles, respectively). In shaded grey, the expectation from MC 
simulations is shown. The width of the reconstructed Z mass distribution is a measure for the 
ID track momentum resolution. The resolution in data with improved alignment is generally 
well reproduced by the MC simulations. However, the resolution is slightly worse in data, so that 
the distribution in the MC simulations needs to be broadened to match the data distribution. 

The right plot of Fig. |2.9| shows the estimated resolution of the vertex position in x-direction 
as a function of the squared sum of the transverse momentuirj^ of the tracks associated to 

the vertex y^'S^T ™ ^-^ million minimum bias events 101 . As expected, the vertex resolution 
improves with increasing \r^p\- The general trend of the data is well described by the MC 



simulations although there are discrepancies in particular for low ^J^p\- The distributions for 
the resolution in y- and z-direction show a similar behaviour 102] . 
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Figure 2.9: Performance of the Inner Detector in pp collisions |102|: the left plot shows the 
invariant mass of fJ-^fJ-^ pairs around the mass of the Z boson from ID tracks 
only in 0.70 fb~^ of data. In shaded grey, the expectation from MC simulations is 
shown. The full and open circles show different versions of the relative alignment 
of the ID subsystems. The right plot shows the estimated resolution of the vertex 
position in x-direction as a function of the squared sum of the of the tracks 
associated to the vertex in 1.5 million minimum bias events. Below the plot, the 
ratio of data and MC simulations is shown. 



The transverse momentum is defined as pT ~ psinO = p/ cosh 77. 
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Figure 2.10: Performance of the electron energy measurement in pp collisions ||103]: the left 
plot shows the correction to the electron energy scale derived from Z — )• ee , 
J/^ — 7- ee and W ^ ev events with the associated uncertainties for < |?7| < 0.6. 
The right plot shows the di-electron invariant mass distribution around the mass 
of the Z boson for data and MC simulations. For both plots, 40 pb~^ of data 
were used. 



The left plot of Fig. 2.10 shows the correction to the electron energy scale derived from Z — )• ee, 
J/'ip — )• ee and W ^ eu events with the associated uncertainties for < |7?| < 0.6 in 40 pb^^ of 
data. The measurements for the different processes are consistent within uncertainties and the 
electron energy scale is known to a level better than 1%. For higher the uncertainties are 
slightly larger |103j. 

The right plot of Fig. 2.10 shows the di-electron invariant mass distribution around the mass 
of the Z boson in 40 pb^ of data and in MC simulations. The width of the distribution is a 
measure of the electron resolution. The MC simulations describe the data well, although the Z 
mass resolution is slightly narrower in data, which can be corrected for by enlarging the electron 
resolution in the simulations. 



Fig. 2.11 shows the total uncertainty on the jet energy scale (JES) as well as different con- 
tributions to the uncertainty as a function of the pi of the jet. The left plot shows jets in the 
region 0.3 < < 0.8, the right plot shows jets in the region 2.1 < \r]\ < 2.8. For the derivation 
of the uncertainties the data taken in 2010 as well as MC simulations were used. 

Contributions to the total JES uncertainty are shown for the underlying MC model (up- 
pointing triangles and solid circles), the noise threshold of the clustering algorithm (down- 
pointing triangles), the in-situ calibration using the balance in j?x of dijet events, called rj- 
intercalibration (crosses and open circles), the response to single hadrons (open squares), and 
the knowledge of the detector material (full squares). 

The total JES uncertainty is of the order of 4% (6%) for low-px jets in the range 0.3 < |r/| < 0.8 
(2.1 < |ry| < 2.8) and decreases to roughly 2% (2.5%) for jets with 60 GeV < pt < 800 GeV. The 
uncertainty in the low-px-region is dominated by MC modelling uncertainties and uncertainties 
due to the ?7-intercalibration. For higher the JES uncertainty increases again to 4% (3%) 
due to the limited knowledge of the response to high-energetic single hadrons. 

The left plot in Fig. 2.12 shows the invariant di-muon mass around the mass of the Z boson 
in 205 pb~^ of data. The expectation from MC simulations is also depicted and shows a slightly 
narrower distribution than observed in data. 
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2.2 The ATLAS experiment 



The right plot in Fig. 2.12 shows the width of the distribution around an invariant mass of 
90 GeV for different regions of the muon 77. The resolution is similar for all ry-regions. Only 
in the barrel-to-endcap transition regions at 1.05 < \rj\ < 1.7, the resolution is slightly worse. 
The resolution is broader in data than in MC simulations consistently for all r/-regions, which is 
taken into account by correcting the muon momentum resolution in the simulations. 
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Figure 2.11: Performance of the jet energy measurement in pp collisions 104 : the total uncer- 
tainty on the jet energy scale as well as different contributions to the uncertainty 
are shown as a function of the pT of the jet. The left plot shows jets in the region 
0.3 < \r}\ < 0.8, the right plot shows jets in the region 2.1 < \r]\ < 2.8. For the 
derivation of the uncertainties the data taken in 2010 as well as MC simulations 
were used. 
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Figure 2.12: Performance of the muon momentum measurement in pp collisions 

left plot shows the invariant di-muon mass around the mass of the Z boson in 
205 pb~^ of data together with the expectation from MC simulations. The right 
plot shows the width of the distribution around an invariant mass of 90 GeV for 
different regions of the muon pseudorapidity. The resolutions measured in data 
are compared to the resolutions from MC simulations. 



23 



3 Description of the analysed data set 



The ATLAS experiment started data taking at = 7 TeV in 2010 and recorded roughly 
45 pb~^ of data in the first year 



106 



In 2011, the instantaneous luminosity of the LHC 
was very much increased and a total of 5.61 fb~^ of data were provided by the LHC throughout 
the year. The data taking efficiency of the ATLAS experiment was very good: 5.25 zt 0.19 fb^^ 



of data were recorded (left plot in Fig. 3.1) 
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Figure 3.1: Left: integrated luminosity in 2011 106 



the luminosity delivered by the LHC is 
shown as well as the luminosity recorded by the ATLAS experiment. Right: peak 
luminosity in 2011 delivered by the LHC ll06j. 



For the analysis presented in this thesis, 1.04 fb~^ of data were analysed. The data were 
taken between March 2011 (data taking period B) and June 2011 (period H). The instantaneous 
luminosity was increased in several steps during this time (right plot in Fig. 3.1 ): the number of 
protons per bunch and the number of bunches present in the LHC were increased, the focusing 
of the proton beams at the interaction point was improved and the bunch spacing time was 
reduced. 

The increasing instantaneous luminosity lead to a larger average number of interactions per 
bunch crossing (pile-up). The increase in the number of protons and the improvements in fo- 
cusing enhanced the contribution from in-time pile-up, i.e. additional interactions in the same 
bunch crossing. Shorter time differences between the bunches and longer bunch trains cre- 
ated additional out-of-time pile-up, which consists of overlaid interactions from different bunch 
crossings. Out-of-time pile- up effects originate from the finite readout time of the subdetector 
systems, which may lead to a wrong assignment of detector signals to bunch crossings and to 
an intrinsic integration of signals from different bunch crossings. 

The average number of interactions per bunch crossing for the data analysed in this thesis is 
shown in Fig. 3.2 It varies between three and eight interactions depending on the data taking 
period. 
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3 Description of the analysed data set 




Average number of interactions per bunch crossing 



Figure 3.2: Average number of interactions per bunch crossing for the 1.04 fb of data which 
were analysed in this thesis. 



Not only the configuration of the LHC proton beams was changing during data taking, but 
also the conditions of the ATLAS detector and readout itself. Although the ATLAS detector was 
operated in a stable mode in general, it was unavoidable that certain failures in the subdetector 
systems temporarily deteriorated the quality of the data taken. Examples are short periods of 
increased noise in the LAr calorimeters and modules in the pixel detector which were temporarily 
unreachable for readout. The effect on the quality of the data caused by each of these defects 
was studied. Based on the defects present in each short time period of data taking (luminosity 
blocks of the order of 1 min), a list of good luminosity blocks in each run of data taking was 
compiled (good runs list). Only data fulfilling the good runs list criteria were taken into account 
for the analysis as well as for the calculation of the integrated luminosity. 

In about 84% of the data considered in this analysis, six front-end boards in the barrel of 
the EMB were not operational due to broken optical links. This issue was permanently present 



until the end of period H and required a special treatment, which is described in Sec. 5.4 It 



was also corrected for in the MC simulations and was considered as a source of systematic 



uncertainty (Sec. 12.3) 



The data were reconstructed with release AtlasPhysics-16 . 6 . 5 . 5 . 1 of the ATLAS recon- 



struction framework ATHENA 107 . The analysis of the data was then performed using C++ 
and the ROOT framework 108]. 

Data were bundled in three different data streams according to the triggers that fired: Egamma, 
Muons and JetTauEtmiss. In order to reduce the amount of data to be processed in each analysis, 
certain preselections (skims) of the different data streams were provided. The following skims 
were used: in the Egamma stream at least one electron candidate with loose shower shape criteria 
and a transverse energ}|^ of -Et > 20 GeV was required, and in the Muons stream at least one 
muon candidate with j?t > 18 GeV was required. In the JetTauEtmiss stream, at least four 
jets with > 20 GeV and two jets with > 40 GeV were required, or at least five jets with 
a PT > 20 GeV. 

The object definitions and requirements on the transverse momenta in the skims of the Egamma 
and Muons streams are looser than the requirements which were actually used for the analy- 
sis (Ch. [5] and [6]). The skim of the JetTauEtmiss stream was intended for tt analyses in the 
all-hadronic channel, but was also used as a control region in this analysis (Ch. [Tjandjo]). 



^The transverse energy is defined in analogy to the transverse momentum; Et — Es'md = iJ/ cosh 77. 
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4 Signal and background modelling 



Simulations of physics processes and detector responses are crucial for the modelling of signal and 
background events in modern high energy physics. As the processes involved can be described 
by probability density functions, Monte Carlo (MC) simulation techniques are a natural choice, 
but also methods to extract contributions from certain processes from data are established. 



According to the factorisation theorem, Eq. (1.2), the hard scattering process can be sepa- 



rated from non-perturbative QCD effects. There are specific programs for the generation of hard 
processes, which can be interfaced with programs which provide models for the non-perturbative 
evolution of the final state including parton showering and hadronisation. There are two main 
models for hadronisation: the Lund string model, which is implemented in the PYTHIA genera- 
tor [109|| , and the HERWIG cluster model, which is implemented in the HERWIG generator [110 
which is commonly used together with the JIMMY generator for multiple parton scattering [111 
These generators also provide models for the underlying event, which adds contributions from 
ISR and FSR, multiple parton interactions, beam remnants and pile-up contributions to the 
final state. 

All particles which are (meta-)stable on time scales of the order of At = Z/c, with / the distance 
of the first detector layer from the interaction point and c the speed of light, were passed to a 



detector simulation using Geant4 [112 , 113 . Geant4 is a general framework for the simulation 
of the interaction of particles with matter based on MC techniques. A detailed description of 
the ATLAS detector was used as an input for the Geant4 simulation [114 . The output of the 



detector simulation was then passed through the same reconstruction software which was used 
for the reconstruction of data. 



The signal process (Sec. 4.1) and a part of the background processes (Sec. 4.2) were modelled 
using MC simulations. However, the MC modelling of certain background processes is not 
sufficiently reliable and methods for the estimation of these contributions from data needed to 



be applied, as discussed in Sec. 11 



4.1 Simulation of tt^ events 

The simulation of the ttj signal process was performed with the WHIZARD MC genera- 



tor [86 , 87 . The full seven-particle final state was calculated in the single lepton and dilepton 
ttj decay channels: Iviqq/bbj and luilujbb-), with l/l = e^, /i^, and I'l/i'i the corresponding 
antineutrino. 

The WHIZARD MC generator was developed especially for automated calculations of matrix 
elements (MEs) in LO. For a given initial and final state it calculates the full ME, taking into 
account all possible contributing diagrams. Hence, all diagrams shown in Fig. 1.5 and 1.6 were 
included together with the different decay modes of the W boson, and interference effects were 
properly taken into account. 

The masses of the u-, d-, s- and c-quarks, and the masses of electrons are much smaller 
than the typical scale of the hard ti'^ process, which is of the order of the top quark mass. In 
order to simplify the calculation of the MEs, the masses of these particles were approximated 
by zero. The masses of muons and r-leptons were set to m„ = 105 MeV and = 1776 MeV, 
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4 Signal and background modelling 



respectively. The masses of b- and top quarks were set to mi, = 4.2 GeV and mj = 172.5 GeV, 
respectively. 

Since the photon is massless, the radiation of photons from charged particles is collinear and 
infrared divergent. Hence, the phase space for photon radiation needs to be reduced in finite 
order calculations in order to obtain a finite cross section for the tt'j process: infrared divergencies 
were avoided by requiring a minimal pT of 8 GeV for photons from the hard process. To avoid 
collinear divergencies, the invariant masses of pairs of massless particles were required to be 
larger than 5 GeV. For consistency, the invariant mass cuts were not only applied for electrons, 
but also for muons and r-leptons. 

The following invariant masses were considered in the single lepton decay mode: m{qi,q2), 
m(gi,7), m(g2,7), "i(^7)> rniQi^j), m{Q2,-f), m{gi,qi), 771,(51,^2), m{g2,qi), and m{g2,q2), 
where qi and q2 are the quarks from the decay of the hadronic W boson, / is the charged lepton 
from the decay of the leptonic W boson, Qi and Q2 are the incoming quarks in the case of 
quark-antiquark annihilation, and gi and 52 are the incoming gluons in the case of gluon-gluon 
fusion. In addition, for each incoming quark i, the invariant mass m{Qi, qj) was considered if qj 
is the antiparticle of Qi. The requirements on Qi and Q2 were dropped for 6-quarks, because 
the latter were not assumed to be massless in the event generation. In the dilepton decay mode 
the invariant mass criterion on m{l,^) was required to hold for both leptons. 

The renormalisation and factorisation scales for the tt'j process were set to 2mt and a LO cross 
section of 0.84 pb was calculated with WHIZARD. As discussed in Sec. 1.3.1 , the fc-factor for the 
tt'j process was estimated to 2.6 it 0.5, which leads to a prediction for the tt'j cross section times 
branching ratio (BR) into the single lepton and dilepton channels of att^ ■ BR = (2.1 it 0.4) pb. 

For the ME calculation, CTEQ6L1 PDFs were used. Modified LO PDFs based on the 



MRST 115 set, called MRST20071omod, available from the LHAPDF package 116 , were used 



for the parton shower generation. Parton shower and underlying event were added to the tt^ 
events using HERWIG (version 6.510) and JIMMY, respectively. 



4.2 Background modelling 

Traditionally, many of the background processes to ti events are estimated from MC simulations, 
because they provide a good description of the data. Backgrounds from multijet events, however, 
are known to be poorly modelled by MC generators and need to be estimated from data. A 
region which is strongly enhanced in multijet production is defined, and the contribution in the 
signal region is obtained with extrapolation methods. Details on the treatment of the multijet 



background in this analysis are given in Sec. 11.2 



At hadron colliders, additional jets are frequently produced by ISR and FSR of quarks and 
gluons. Hadrons in these jets can be misidentified as photons (hadron fakes) and therefore 
analyses with photons in the final state typically feature sizable background contributions from 
this source. 

For two reasons, the use of MC simulations for the hadron fake contribution is disfavoured 
with respect to estimates from data: firstly, the exact simulation of geometric shapes of the 
electromagnetic clusters in the calorimeter (shower shapes) requires very detailed detector un- 
derstanding and is hence challenging jl03| . Secondly, the description of jet fragmentations with 
a leading neutral hadron (vr'', tj, . . .) by MC simulations is known to be difficult. Such neutral 
hadrons are likely to give rise to photon-like signals via the decay into two photons. 
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4.2 Background modelling 



Hence, the strategy of this analysis was set up in a way to minimise the dependence on MC 
simulations (Ch. [7]). In particular, the background from hadrons misidentified as photons was 
estimated from data. 

The estimates from MC simulations only make up a small part of the final background es- 
timate, but MC simulations were also used for cross-check studies. The MC modelling of the 
different background contributions is described in the following. 



Top quark pair production 



tt events were produced with the MC@NLO |117| generator (version 3.41) using CTEQ6.6 [118 



PDFs. For the simulation of the parton shower and the underlying event, MC@NLO was inter- 
faced to the HERWIG (version 6.510) and JIMMY generators, for which the AUETl tune to 
ATLAS data |119| was used. The ti cross section was calculated to 165]tig approximate 
NNLO with HATHOR [42 1. 

Photons are also produced in the ti simulation, which may lead to a ti'y signature: two kinds 
of processes can occur: HERWIG produces real photons in the fragmentation processes and 
also allows for QED corrections in the production and decay of the top quark pairs using the 



PHOTOS package 120 . In order to avoid double-counting of tt'^ events in the samples generated 



with WHIZARD and MC@NLO, contributions which fulfil the requirements for the signal phase 



space, as defined in Sec. 4.1, were removed from the MC@NLO sample. 



The definition of the signal phase space involves the invariant mass cuts described in Sec. 4.1 



Due to the different handhng of photon radiation in WHIZARD and HERWIG+PHOTOS, the 



application of these cuts is not trivial as illustrated in Fig. 4.1 in WHIZARD, photon radiation 
is handled as part of the ME calculation and the invariant mass cuts are applied to the seven- 
particle final state (left plot). In the MC@NLO tt sample, photon radiation is added a posteriori. 
Photons can be radiated from the incoming quarks, the top quarks, the W boson or the decay 
products of the W boson (not illustrated in the figure). 

Particularly different with respect to WHIZARD is the treatment of the radiation from quarks, 
which is part of the parton shower process in HERWIG (right plot): a quark from the W decay 



(HERWIG status code 123 or 124) is translated to a jet four- vector (status code 143 or 144) [110 



The constituents of the jet, which may comprise photons, are listed as the decay products of the 
jet four-vector. Thus, photon radiation is not treated in a single step, but as part of the parton 
shower process, and the definition of the quark to be considered for m{q, 7) is ambiguous. 




q (123) 



W 




Figure Illustration of the different treatment of photon radiation from quarks in tt events 
in the WHIZARD (left) and MC@NLO (right) MC generators. 
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4 Signal and background modelling 



The following approach is considered closest to the treatment in WHIZARD: the four- momenta 
of all particles in the parton shower are added up except for the photon. The combined four- 
momentum is then used as an estimate for the four-momentum of the quark after photon 
radiation in the evaluation of m{q,^). However, a systematic uncertainty is assigned to this 
ambiguity as described in Sec. |11.2[ 

For the study of systematic uncertainties related to MC generators, additional samples for 
the ti process were available. POWHEG [121 was used as an alternative NLO generator. It 
was interfaced with HERWIG and JIMMY as well as with PYTHIA (version 6). For studies 
of the amount of ISR and FSR, different samples with modified values for the ISR and FSR 
parameters were generated with AcerMC 122] interfaced to PYTHIA. The parameters were 
varied in ranges currently not excluded by experimental data 123] , comparable to those used in 



the Perugia Soft/Hard tune variations 124 



Single top production 

Single top processes were simulated with the MC@NLO generator interfaced to HERWIG and 
JIMMY. In order to avoid an overlap of the M^t-channel with production, the diagram-removal 
scheme described in Ref. [125 was applied. 

The single top cross sections were taken from the approximate NNLO calculations in t he t- , 
s- and PFt-production channels, which yielded 64.6120 P*-"' ^-^ ^ ^-^ P'-' 15.7 ± 1.1 pb 
1281. 
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W boson production in association with jets (VV^-|-jets) 



VF+jets events were simulated with the ALPGEN [129] generator using the CTEQ6L1 36 PDFs 



For the simulation of the parton shower and the underlying event, ALPGEN was interfaced to the 



HERWIG and JIMMY generators, for which the AUETl tune to ATLAS data for CTEQ6.1 [119 
was used. 

The process was simulated for different parton multiplicities in the ME, and the matching 
to the parton shower was applied exclusively for parton multiplicities smaller than five and 
inclusively for the 2 — )• 7 process with five partons in the final state. In order to enhance 
statistics for W boson production processes with additional heavy flavour jets (c- or 6-jets), 
additional samples for the processes W^+c+jets, W+cc+jets, and W+bb+iets with up to five 
partons in the final state were produced. The overlap with the fiavour inclusive samples was 



removed in order to avoid double-counting [130 



The cross sections of the different parton multiplicity samples were calculated in NNLO using 
the FEWZ |131j and ZWPROD |132j codes. One inclusive M^+jets sample was constructed by 
weighting the individual samples corresponding to their integrated luminosities. 



The uncertainty on the cross section was evaluated using Berends-Giele scaling 133 , 134 : the 
uncertainty increases with the number of jets and was estimated to 24% for each jet in the final 
state, to be treated uncorrelated. Hence, for events with at least four jets, an uncertainty of 
48% was obtained. An additional uncertainty on the amount of the fraction of b- and c-quarks 
produced in association with the W boson was added in quadrature: the fraction of events 
with a bb or cc pair was scaled by a factor of 1.63 it 0.76 with respect to the generation with 
ALPGEN. The fraction of events with only one c-quark was scaled by a factor of 1.11 it 0.35. 
The fraction of events without c- or 6-quarks was corrected correspondingly in order to preserve 
the total predicted number of VF+jets events. 
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4.2 Background modelling 



Z boson production in association with jets (Z+jets) 

Z+jets production was treated similarly to W^+jets production: ALPGEN, HERWIG and 
JIMMY were used with the same PDFs and generator tunes. Also, the same scheme for the 
generation of samples with the different parton multiplicities was applied. The statistics for 
Z+jets production with heavy flavour jets was enhanced by the use of Z+66+jets samples with 
up to five partons in the final state. The overlap with the fiavour inclusive samples was removed. 

As for 1^+jets production, the cross sections were calculated in NNLO with the FEWZ and 
ZWPROD codes, and one inclusive Z+jets sample was constructed by a reweighting of the 
individual samples. 

The uncertainty on the cross section was estimated to be the same as for T^+jets production, 
that is 48% for events with at least four jets in the final state. 



W boson production in association with jets and photons (Vl^+jets+7) 

events with additional jets in the final state were generated with ALPGEN. Parton shower 
and underlying event were added using HERWIG and JIMMY, and the CTEQ6L1 PDFs were 
used - as for the VF+jets sample. As for the VF+jets generation, the matching of the parton 
shower to the ME was performed exclusively for the zero to four parton samples, and inclusively 
for the five parton sample (2—7-8 process). To avoid infrared divergencies, a minimum px of 
10 GeV was required for the photon. 

The relative weights of the different parton multiplicities were taken from the ALPGEN event 
generation. The absolute cross section for the Ty+jets+7 sample was used for the measurement 
presented in this thesis. 



Diboson production {WW, WZ, ZZ) 

WW, WZ and ZZ events were generated with HERWIG. The cross sections calculated with 
HERWIG were corrected by fc-factors obtained with the MCFM code (135|, which read 1.48 for 
WW, 1.30 for WZ, and 1.60 for ZZ production. The uncertainties on the diboson cross sections 



were estimated to 5% following the approach in Ref. 136 



For the simulation of pile- up effects, minimum bias events generated with PYTHIA were 
overlayed to the hard processes. The pile- up configuration corresponded to a bunch spacing of 



50 ns, which represents most of the data taken in 2011 95 . The pile-up rate was kept variable 
in the simulation of the minimum bias events, and the MC samples were reweighted so that the 
distribution of the number of interactions per bunch crossing was the same in MC simulations 



and in the data analysed (Fig. 3.2) 
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5 Object definitions 



Particles traversing the detector produce characteristic signatures in the different subdetectors 
(Sec. 2.2), which are used to identify the particle type. Combined objects were constructed 



using information from the Inner Detector (ID), the calorimeter system and the muon spec- 
trometer (MS), which were then compared to the object definitions for electrons (Sec. 5.1), 
muons (Sec. |5.2[ ), jets (Sec. 5.3) and photons (Sec. 5.6). These definitions combine a high prob- 
ability for the identification of real electrons, muons, jets and photons with a low probability of 
misidentifying another object as one of these particles. Moreover, 6-tagging was used to identify 
jets from 6-quarks (Sec. 5.5). 

Since the initial momentum of the colliding partons in the transverse plane is small, an imbal- 
ance of the transverse momentum of the whole event (missing transverse energy ^t) indicates the 
presence of high-energetic, undetected particles, such as neutrinos. The energy and momentum 
measurements in the whole detector were used to measure (Sec. 5.4). 



5.1 Electron definition 



Electron candidate objects 103 were selected by searching for large energy deposits in the elec- 



tromagnetic (EM) calorimeter using a fixed-size window in ry-c/i-space with an ID track pointing 
in its direction. The energy deposited in the EM cluster was corrected for energy losses in front 
of the calorimeter, lateral leakage outside of the cluster window and energy deposited behind 
the calorimeter. The electron four-vector was built from the cluster energy and the direction 
measurements from the ID track. 

Only the central part of the calorimeter was used and | ?7ciuster I was required to be smaller than 
2.47. The transition region from the barrel to the endcap calorimeter, 1.37 < |77ciuster| < 1-52, 
was not considered. The i?T of the electron was required to be larger than 25 GeV in the 
analysis, but also electrons with an Eic down to 15 GeV were used for the estimation of the 
photon isolation properties (Ch. [s]). 

In order to suppress backgrounds from other particles misidentified as electrons, cut-based 
sets of quality criteria (menus) provided increasing background rejection: the so-called loose, 
medium and tight menus. In the following, the tight menu is described, which yields an overall 
efficiency of roughly 75%. 

Electron clusters tend to be smaller in size than clusters from hadrons within jets. Hence, 
several observables constructed from the geometrical shape of EM clusters (shower shapes), such 
as their lateral width or the energy in the highest-energetic calorimeter cells, were used in the 
tight menu. Also, the fraction of the energy deposited in the hadronic calorimeter, which 
is typically very small for electrons, was exploited to suppress backgrounds from jets. The 
transition radiation in the TRT was used in addition to discriminate against charged hadrons. 
In order to assure that tracks are not accidentally associated to clusters, quality criteria on the 
number of hits in the silicon trackers, a good geometrical matching of the track direction and 
the cluster position as well as of the track momentum and the cluster energy were required. The 
track also had to point back closely to the primary vertex. Backgrounds from converted photons 
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5 Object definitions 



(7 —7- e^e~) were suppressed by requiring a hit in the Pixel 6- layer, because most photons do 
not convert before they reach this first detector layer. 

Electrons which were also reconstructed as photon candidates were not considered, but treated 
as photons (cf. photon recovery procedure in Sec. 5.6). To further suppress backgrounds from 
jets, the energy of the electron candidate in a cone of AR = 0.2 around the electron directioi|^ 
(isolation energy) was required to be less than 3.5 GeV. The isolation energy was corrected for 
average energy deposits from pile-up events. Since the electron identification and jet finding 
algorithms (Sec. 5.3) are independent of each other, most electrons were also reconstructed as 
jets. In order to avoid double-counting, jets which were closer than 0.2 in r]-(j)-spsLce to an electron 
were disregarded. Electrons close to a region in the EM calorimeter which was known to feature 
a broken optical link (Ch. [s]) or a dead high- voltage channel were ignored. MC simulations were 
corrected for this effect. 

Electron signatures were used in the trigger menus at LI, L2 and EE (Sec. 2.2.5). At LI, 
trigger towers above certain £'t thresholds were searched for. At L2, a simplified version of the 
offline reconstruction algorithm was used. The final trigger decision at EE level used the full 
offline algorithm with slightly looser requirements. Eor the data analysed, the EF_e20_mediuin 
trigger was used for the selection of candidate events in the electron channel. This trigger 
required a minimum Eic of 20 GeV at EE level and the medium trigger shower shape menu. The 
thresholds at LI and L2 had to be adjusted during data taking to maintain manageable event 
rates given the increasing instantaneous luminosity. 

Trigger, reconstruction and identification efficiencies were measured in data using tag-and- 
probe methods in Z — )• ee and W ^ ev events jl03j : while one electron was used to define the 
sample (tag), the efficiencies of the other electron (probe) were measured. In case of W events, 
was used as a tag. Discrepancies between data and MC were accounted for by scale factors 
(SEs): 'SF[Ei,ri) = eAataiETirj) / emc^Ei^ti), where ffdata and emc are the efficiencies in data 
and MC, respectively. MC events were weighted with the SEs, so that the efficiencies in the MC 
simulations yielded those in data. 

The SEs for the EF_e20_medium trigger were derived in 18 bins in rj and vary between 0.97 
and 1.00 with uncertainties smaller than 0.01. The SEs for the electron reconstruction were 
derived in three bins in \r]\ and vary between 0.98 and 1.01 with uncertainties smaller than 0.02. 
Eor the electron identification, the SEs were divided into 18 bins in r] and five bins in E^. The 
range for these SEs is 0.95 - 1.12 with uncertainties smaller than 0.04. 



The energy scale and resolution of electrons was studied in data using Z — )• ee events 103 
In order to match the distribution of the di-electron invariant mass in data, the electron energy 
resolution was corrected in MC simulations. Additionally, small corrections to the energy scale 
in data were applied. 



5.2 Muon definition 

Muons were reconstructed using tracks in the MS and in the ID using the Muld algorithnj^ 137 



138 . The pattern recognition started from track segments in the MS, which were extrapolated 
to the ID. If a matching ID track for a MS track segment was found, a combined fit to the ID and 
MS muon hits was performed and a combined muon object was built. The muon momentum was 
measured from the combined fit. Only muons with a px of at least 20 GeV fulfilling |r/| < 2.5 
were considered in this analysis. 



^A7? = -J IS.rj^ + A(/!>2, with Ary = jyceii — rj and Ai^ = 0coii — 4>- 

■^The Muld algorithm is also called "chain 2" reconstruction in Ref. 137 and 



138 
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Several isolation criteria were applied in order to suppress backgrounds from muons produced 
in jets, in particular muons from leptonically decaying 5-mesons in 6-jets. The energy deposited 
in the calorimeter in a cone of AR = 0.3 around the muon candidate was required to be smaller 
than 4 GeV. Also the sum of the transverse momenta of the tracks in the ID in a cone of 
AR = 0.3 around the muon track had to be smaller than 4 GeV. Additionally, if a jet with min- 
imum of 20 GeV was closer than 0.4 in r/-(/>-space to the muon, the muon was not considered 
isolated and was disregarded. 

Muon signatures were used in the trigger menu at LI, L2 and EF (Sec. 2.2.5). The triggers 
at LI were based on coincidences in rj and (p of hits in the different layers of the RFC and TGC 
systems requiring a minimum p^^ of the trigger object candidate. At L2 and EF level, track 
finding algorithms using information from the MS and the ID were used. In this analysis, the 
EF jttulS trigger chain was used which required a combined muon with > 18 GeV at EF level. 

Trigger, reconstruction and identification efficiencies were measured in data using tag-and- 
probe methods in Z — )• jU^//~ events. Discrepancies between data and MC were accounted for by 



SFs 138 . MC events were weighted with the SFs, so that the efficiencies in the MC simulations 
yielded those in data. Events in which the muon object had a pT larger than 150 GeV were 
disregarded, because muon trigger SFs were not available for larger due to an issue with the 
trigger modelling in MC simulations in this region. 

The SFs for the EFjiuIS trigger were derived in three bins in p^, seven bins in rj and three 
bins in (p. For px < 120 GeV, they vary between 0.88 and 1.21 with uncertainties up to 0.09 in 
the barrel, and between 0.98 and 1.04 with uncertainties smaller than 0.02 in the endcap region. 
For larger transverse momenta, the measurement from Z boson decays was limited by statistics 
and SFs down to 0.66 and up to 1.48 were derived in certain regions of the phase space with 
correspondingly large uncertainties. The SFs for the muon reconstruction were derived in 20 bins 
in rj and vary between 0.98 and 1.01 with uncertainties which do not exceed 0.005. For the muon 
identification, one inclusive SF was derived, which reads: 1.0008 ± 0.0003 (stat.) ± 0.0003 (syst.). 

The momentum scale and resolution of muons were measured in data using Z — t- fJ,^fJ,~ events 



139 . In order to match the distribution of the di-muon invariant mass in data, the muon 



momentum scale and resolution were corrected in MC simulations. 



5.3 Jet definition 



Jets were reconstructed with the anti-fc^ algorithm 140, 141 with a distance parameter of 



R = 0.4. Jets were built from topological calorimeter clusters based on the significance of the 



energy deposit in the calorimeter cells with respect to their noise level 1142 . The clusters, and 



hence the jets, were calibrated to the EM scale, which is smaller than the hadronic energy scale 
due to the non-compensating nature of the ATLAS calorimeters. The hadronic energy scale was 
restored using MC based correction factors depending on the jet px and rj fl04]. In this analysis, 
jets were required to have a minimum px of 25 GeV. This requirement was lowered to 20 GeV 
for the muon isolation criterion. 

In order to ensure that the selected jets originated from the hard scattering process, quality 
criteria 1041 were applied to all jets with a minimum px of 20 GeV. These criteria strongly 



reduce beam-gas and beam-halo events arising from the interactions of the proton beam with 
gas in the beam pipe and with the collimators upstream the detector, respectively. They also 
suppress backgrounds from cosmic ray muons and from calorimeter signals which were due to 
large electronics noise. Since all of these background contributions tend to create artificial high- 
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Pt signals in the detector, the rejection of events with jets faihng the criteria improves the 
quahty of the measurement. 

The uncertainty on the jet energy scale (JES) was derived from the response of single hadrons 
from test beam measurements, in-situ techniques applied in collision data, and input from MC 
simulations [104 . Dijet events were used to test the relative energy scale of the different regions of 
the detector given the varying design of the calorimeter and material in front of the calorimeter. 
The JES was validated in-situ comparing the jet energy to the momentum carried by the tracks 
associated to the jet, and by exploring the balance in of 7+jet events, photon events with 
hadronic recoil, and multijet events with one high-px jet 104 . The total JES uncertainties in 
the ranges 0.3 < \r]\ < 0.8 and 2.1 < |?7| < 2.8 were presented in Fig. 2.11 



The jet reconstruction efficiency was measured comparing jets built from calorimeter clusters 

The jet energy resolution was estimated using the px-balance in 



to jets built from tracks 104 
dijet events jl43j . 



5.4 Definition of the missing transverse energy 



The measurement of 144 was based on the energy measurement in the calorimeters and the 
muon momentum measurements in the MS. In order to calibrate the different energy depositions 
in the calorimeter correctly, different objects were identified and their energy was taken into ac- 
count at their correct energy scale. Firstly, electrons and photons with a minimum i^x of 10 GeV 
were identified requiring the respective tight shower shape menus (^E^^^'^"^^ and Ex]^y^°^^^ . 
Secondly, anti-kt jets with a pT of at least 20 GeV were taken into account at the hadronic 
energy scale (^E^x^y^ ■ Jets with transverse momenta between 7 and 20 GeV (soft jets) were 
not calibrated to the hadronic scale, but were included in the definition at the EM scale 
[Elf^'^'Y The energy in topological clusters not used in any of the reconstructed objects was 

also included at the EM scale (£;^^^'°"*). 

Muons reconstructed with the Muld algorithm were taken into account up to \rj\ < 2.7. They 
were required to feature a combined ID and MS track for < 2.5. The energy deposition in 
the calorimeter around the muon track was included in case the muon was separated from any 
anti-fcf jet by a minimum distance of 0.3 in 7/-(/)-space (-E™^""^). Hence, the full definition of the 
read^ 

771 / TTielectrons i i^photons i ipiets i iTTSoft jets i ^muons i iricell out i 

= v^7+^7- (5-1) 



As mentioned in Sec. 5^, the quality of the measurement was assured by rejecting events 
where a jet with p^ > 20 GeV did not fulfil the jet quality criteria. Also events with jets with 
Pt > 20 GeV closer than 0.1 in r/-(/>-space to a region which was known to feature a broken 
optical link (Ch. [s]) or a dead high-voltage channel in the LAr calorimeter were disregarded. 
MC simulations were corrected for this effect. Events with large noise in the LAr calorimeter 
were also rejected. 

The object definitions used for the calculation were chosen consistent with those described 



in Sec. |5.1[ |5.2[ 5.3 and 5.6, and the uncertainties of the energy and momentum scales of the 



different contributions were propagated consistently to the ^t. 



Ex and Ey are defined in analogy to the transverse energy Et- 
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5.5 6-tagging 

6-tagging was used to identify jets originating from the 6-quarks of the top quark decays. The 
relatively long lifetime of i?-mesons leads to displaced secondary vertices and impact parameters 
of their decay products. Algorithms were constructed to discriminate b- from light jets exploiting 



these properties. Two 6-tagging algorithms were combined in this analysis 145 ,146]: I3PD and 
JetFitter. 

JetFitter makes use of the topology of weak b- and c-hadron decays inside jets to discrim- 
inate between b-, c- and light jets. Different properties of the reconstructed decay vertices, 
such as their masses, momenta, flight-length significances and track multiplicities, were used 
to build a likelihood discriminant. The second algorithm, I3PD, uses the transverse and lon- 
gitudinal impact parameter significances do/crfo ^o/f^zo of the tracks associated to the jet. 
Since both algorithms use different observables sensitive to the jet fiavour, their combination is 
advantageous. A single discriminating observable for 6-tagging was constructed using a neural 
network. 

For this analysis, a working point for the neural network output was chosen such that ap- 
proximately 70% of the 6-jets in simulated ti events were correctly 6-tagged. This corresponds 
to misidentification rates of roughly 20% for c- and 1% for light flavour jets. 



The 6-tagging efficiency was measured in data using different techniques [147 , 148 . These ex- 



ploit in particular the relative of the muon from semileptonic 6-decays and the reconstruction 
of D*~^ mesons from b — )• Xfj,D*^ decays. Mistagging rates were measured using the invariant 
mass of the tracks associated to a secondary vertex and events with negative impact parame- 
ters. From these measurements, SFs were derived to reproduce the data tagging and mistagging 
efficiencies in MC simulations. 

The SFs for the 6-tagging efficiency were derived in five bins in p^, and vary between 0.94 
and 0.99 with uncertainties between 0.06 and 0.15. The mistagging SFs were estimated in eight 
bins in for jets with |?7| < 1.2 and 1.2 < \r]\ < 2.5, and they vary between 0.98 and 1.27 with 
uncertainties between 0.11 and 0.22. 



5.6 Photon definition 

The reconstruction of photon and electron candidate objects is done by one single egamma 



algorithm, so that objects are unambiguously identified as either photons or electrons 149 
The egamma algorithm takes into account that photons may convert into e^e~ pairs by the 
interaction with the material in front of the calorimeter, and different reconstruction paths are 
foreseen for unconverted and converted photons. While unconverted photons do not feature a 
track pointing to their EM clusters, converted photons typically have two tracks pointing to the 
cluster (two-track conversions). There are also asymmetric conversions into an e~^e~ pair with 
two very different transverse momenta, so that the low-pT track is likely not to be reconstructed, 
which leads to one-track conversions. The tracks of one-track conversions often miss hits in the 
first ID layers, depending on where the conversion into the e^e~ pair took place, and also energy 
and momentum measurements tend to be not consistent with an electron object. 

Unconverted photons were reconstructed by finding fixed-size tower clusters with significant 
energy in the EM calorimeter and no ID track matched to it. Clusters with tracks associated to 
them were treated as electrons/conversion candidates. The next step was the photon recovery 



procedure 149 , which resolves the ambiguity between electrons and photons. This is particularly 
important for converted photons, but also unconverted photons with erroneously associated 
tracks need to be recovered. 
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Conversion vertices were reconstructed by either fitting two tracks under the assumption that 
they originated from a massless particle or by identifying single tracks that did not feature hits 
in the inner layers of the ID. EM clusters were then checked for conversion vertices matching 
the cluster centre in r]-(p-space when extrapolated to the calorimeter surface. 

Converted photons were identified by comparing the tracks from the associated conversion ver- 
tices to the track which matches the cluster best. Additionally, clusters with the best matching 
high-pT track without hits in the silicon trackers (TRT-only track) were considered as conver- 
sions candidates. Refined requirements were applied using the ratio of calorimeter energy and 
track momentum (E/p), and the presence of hits in the Pixel 6-layer. In addition, unconverted 
photons with low-px TRT-only tracks and low-px tracks with a large E/p ratio were identified. 

A special treatment was applied to converted photon candidates with a track missing a hit in 
the Pixel 5- layer, which would have been identified as an electron candidate if the hit had been 
present: if the track passed through a module in the 6-layer which was known to feature readout 
problems during data taking, the candidate was not considered a photon. 

The photon energy was calibrated using the energy deposited in the EM calorimeter including 
the presampler. Cluster sizes of 3 x 5 and 3x7 cells were used in the second layer of the barrel 



calorimeter for unconverted and converted photon candidates, respectively 123 . A larger cluster 



size was used for converted photons to account for broader showers in cp due to bremsstrahlung 
from the e~^e~ pair. The cluster energy was corrected for energy losses in front of the calorimeter, 
lateral leakage outside of the cluster and energy deposited behind the EM calorimeter. These 
corrections were parametrised as a function of the energy depositions in the presampler, the 
three calorimeter layers, and the rj of the photon object. 

Photon four- vectors were constructed from calorimeter information only. Assuming a photon 
mass of zero, the cluster energy was used together with the rj- and (/>-position of the cluster in 
the second calorimeter layer, in which the bulk of the photon energy is typically deposited. 

Only photons with \r]\ < 2.37 were considered, thus limiting the acceptance to the pseudora- 
pidity range with particularly fine granularity in the first calorimeter layer (LAr strips). The 
transition region from the barrel to the endcap calorimeter, 1.37 < \r]\ < 1.52, was not consid- 
ered. Photons were required to have a minimum E't of 15 GeV. 

A loose and a tight cut-based menu using photon shower shapes and hadronic leakage were 
used to purify the sample of selected photon candidates. In the following, only the tight menu 
is described. Since the clusters of unconverted and converted photons had slightly different 
properties, the cuts on the various observables used in the tight menu were tuned separately 
for these two kinds of photons. The cuts are binned in seven |r/|-regions accounting for the 
varying amount of material in front of the calorimeter and the different granularities: [0.0, 0.6), 
[0.6,0.8), [0.8,1.15), [1.15,1.37), [1.52,1.81), [1.81,2.01), and [2.01,2.37). An overview of the 
different observables used in the tight menu including the definitions and symbols used further 



on is given in Tab. 5.1 



As for electrons, the energy in the hadronic calorimeter behind the EM cluster (hadronic 
leakage) and the cluster widths in the second layer of the EM calorimeter were used to discrim- 
inate real photons from hadrons produced in jet fragmentation. The latter tend to deposit a 
significant amount of energy in the hadronic calorimeter and produce broader clusters. 

In order to suppress backgrounds from jets which fragment with a high-energetic hadron 
decaying to two photons (vr*^ mesons, t] and rj' mesons ...), the first calorimeter layer (LAr 



strips) was used. The LAr strips are finely granulated in rj, as illustrated in Fig. 5.1 The 



photons from the meson decays are typically very close to each other and therefore give rise to 
one single cluster in the EM calorimeter. However, this cluster originating from two photons 
tends to be slightly broader than single photon clusters, which is exploited by the observables 
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Category 


Symbol 


Description 


Hadronic leakage 


•'iiaai 


Ratio of tlie Ej' in the first layer of the hadronic 
calorimeter to the E^ of the EM cluster (used for 
\r]\ < 0.8 and |r?| > 1.37) 




-Rhad 


Ratio of the E^ in the whole hadronic calorimeter to 
the E-r of the EM cluster (used for 0.8 < |r/| < 1.37) 


Second calorimeter layer 


Rri 


Ratio of the cell energies in 3 x 7 and 7x7 cells in ri-(j) 






Lateral width of the shower 




R4> 


Ratio of the cell energies in 3 x 3 and 3x7 cells in rj-cp 


First calorimeter layer 


Ws3 


Shower width for three cells around the maximum cell 


(LAr strips) 




in the first layer 




Ws tot 


Total shower width in rj in the first layer 




-fgide 


Fraction of the energy outside of the core of the three 
central cells but within seven cells in the first layer 
(in rj) 




AE 


Difference between the energy in the second maximum 
cell and the energy reconstructed in the cell with the 
minimal value found between the first and second max- 
imum cells in the first layer 




-E'ratio 


Ratio of the energy difference of the two largest energy 
deposits along r/ in the first layer over the sum of these 
energies 



Table 5.1: Overview of the shower shapes used in the tight menu 149 



WsS, yJs tot and .Fside, and also tends to have a second maximum within the cluster, to which AE 
and -Bratio are sensitive. 

No isolation criterion was included in the photon definition, although real photons are expected 
to feature significantly less activity in a small cone around the photon candidate than fake 
photons from hadrons inside jets. However, the track isolation p^^^"^^, as defined in Ch. [Tj was 
used in a template fit to estimate the amount of hadrons misidentified as photons directly from 
data. 

Photons close to a region in the EM calorimeter which was known to feature a broken optical 
link (Ch. [3]) or a dead high-voltage channel in the LAr calorimeter were ignored. MC simula- 
tions were corrected for this effect. Moreover, EM clusters with very narrow energy deposits 
and a large contribution from untypical electronics pulse shapes in the LAr calorimeter were 
disregarded, because they were most probable to originate from large electronics noise. 

As mentioned already in Sec. 5.1, the eganima algorithm and the jet finding algorithm are 
independent of each other. Hence, not only electrons are double-counted as jets, but also pho- 
tons. In order to avoid this effect, the corresponding jets needed to be removed, which was done 
based on a geometrical matching in 7/-(/)-space of jets and photons: jets which were closer than 
AR = 0.1 to photons were disregarded. 

The left plot in Fig. 5.2 shows the distance in r/-0-space between photons and the closest 
reconstructed jet in simulated tt'j events. The photon objects were required to be real photons, 
that means originating from true simulated photons. For roughly 85% of the photons, the 
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closest jet centre differs by less than 0.1 in ry-0-space, which indicates that the photon was also 
reconstructed as a jet and was hence double-counted. 



The right plot in Fig. 5.2 shows the average ratio of the transverse energy of real photons and 
the transverse momentum of the jet closest to the photon as a function of the distance between 
the photon and the jet. The jet was taken on the EM scale to be comparable to the photon 
i?T- In cases where the jet axis was very close to the photon direction, the jet pic was found 
to be similar to the photon E^, which means that the jet does not contain additional particles 
and mainly consists of the photon object. If the jet algorithm picks up additional particles close 
to the real photon, the pT of the jet exceeds the photon Eic and also the jet axis differs from 
the photon direction. Larger distances between photons and jets indicate that more additional 
energy was included in the jet on average. 

Hence, there are two effects which needed to be addressed: the double-counting of photons 
as jets was avoided by removing jets closer to good photons than AR = 0.1. In these cases, the 
photon was found to be still the main component of the jet comparing photon and jet px- 
The second effect was the presence of particles in the proximity of the photon, which happened 
when a real jet was so close to the photon that both particles were reconstructed as one single 
jet. This was indicated by larger AR values for the closest jet. Since such events featured two 
overlapping objects, just removing the jet would have biased the jet reconstruction. Hence, the 
whole event was removed if it featured a AR value between 0.1 and 0.5. 



The upper limit of AR = 0.5 was motivated by the left plot in Fig. 5.3 , which shows the photon 
identification efficiency for real photons in tt'j simulation as a function of the AR between the 
photon and the closest jet after jets with AR < 0.1 were already removed. The identification 
efficiency is defined with respect to reconstructed photons. It drops significantly for small 



It is emphasised tiiat both, jet pr and photon Et, are based on energy measurements in the calorimeters. 
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Figure 5.2: Left: the distance in 77-(/)-space between real photons and the closest reconstructed 
jet in tt7 simulation. The effect of double-counting real photons by the jet re- 
construction is clearly visible. The last bin includes the overflow bin. Right: the 
average ratio of real photon E't and closest jet on electromagnetic scale as a 
function of the distance between photon and jet. Jets with centres that differ from 
the photon direction include energy deposits from additional particles. The cut 
applied at Ai? = 0.1 is indicated by a vertical line. 



distances in ?7-(/)-space, which indicates a distortion of the photon shower shape variables used 
in the tight identification menu by the energy deposits of surrounding particles from a close- by 
jet. The efficiency was found to be close to constant for Ai? > 0.5, which was hence chosen as 
the upper limit for events to be disregarded due to the proximity of jets to photons. 

Photons and electrons as well as photons and muons were implicitly separated in ?7-(/)-space by 
their respective object definitions: the double-counting of electromagnetic clusters as electrons 
and photons was avoided by the photon recovery procedure described above. Muons were not 



considered if they were closer than 0.4 in 7/-(/)-space to a jet (Sec. 5.2). Since photons were also 
reconstructed as jets, this translated to a minimal Ai? requirement between muons and photons. 
Fig. |5.4| shows the normalised distributions of the distance in r/-(^-space between photons and 
electrons (left), and between photons and muons (right), respectively, in simulated ti^ events. A 
full event selection as described in Ch. [6]was applied. It can be seen that photons and electrons 
are separated by a Ai? of at least 0.2. The minimal distance between photons and muons is of 
the order of 0.4. 



The right plot in Fig. 5.3 shows the identification efficiency as a function of the photon 
i?T after the requirements on the photon-jet distance were applied. The efficiency rises with 
increasing photon i?T and approaches a constant value of roughly 80%. For Iow-E't photons, the 
efficiency was found to be significantly lower and only as large as roughly 55%. The reason is the 
decreasing discriminating power between showers from real photons and from fake photons from 
vr" — )• 77 decays, because the shapes of low-energetic showers are subject to stronger statistical 
fluctuations than high-energetic showers. 

Since there is no easily accessible reference process for photons in the energy range used in this 
analysis, the measurement of the photon identification efficiencies as well as the energy scale and 
resolution is not trivial. For the correct simulation of the identification efficiencies, the shower 
shapes used in the tight menu were shifted in the MC simulations so that their efficiencies 
matched the ones in data. This is a procedure which was adopted before in the isolated prompt 
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Figure 5.3: Left: photon identification (tight menu) efficiency for real photons in tt^ sim- 
ulation as a function of the AR between photon and closest jet with a minimal 
distance of AR = 0.1. The efficiency drops significantly for small distances in 
77-(/>-space. The cut at AR = 0.5 is indicated by a vertical line. Right: photon 
identification efficiency for good real photons in ^7 simulation as a function of the 
photon E'Y- Only photons with a minimal distance of AR = 0.5 to the closest jet 
were taken into account. 
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Figure 5.4: The left plots shows the normalised distribution of the distance in r/-0-space be- 
tween photons and electrons after a full event selection (Ch. [g]) in simulated ti^ 
events. The right plots shows the distance between photons and muons. Photon 
and electron objects are separated by a AR of at least 0.2. The minimal distance 
between photons and muons is of the order of 0.4. In both plots, the last bin 
includes the overfiow bin. 



photon [153 and the isolated diphoton cross section measurements It accounts for the 

average mismodelling of the shower shapes in the simulation. 

The shifts were obtained by comparing the shower shape distributions for photons in data to 
the distributions obtained from true photons in MC simulations, where the photons originated 
from jet fragmentation or 7+jet processes. For this comparison, the energy in a cone of 0.4 in rj- 
0-space around the photon candidate was required to be less than 3 GeV, so that a purer sample 
of photon candidates was selected. Fig. 5.5 shows the identification efficiencies as used in the 
diphoton measurement 151 for the tight menu in the isolated photon sample for unconverted 
and converted photons in bins of E't and ry. The values are larger than those given in the right 
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Figure 5.5: Identification efficiency for unconverted and converted photons in bins of i?T and 
r] derived for the isolated diphoton cross section measurement 151 . The isolation 
energy in the calorimeter in a cone of size Ai? = 0.4 around the photon was required 
to be less than 3 GeV. The efficiencies themselves were not used in this analysis 
but the relative systematic uncertainty was used. 
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plot in Fig. 5.3, because of the correlation between calorimeter isolation and the shape of the 
cluster in the calorimeter. 

In this analysis, the same shower shape shifts as in Ref. [m] were used to correct the efficien- 
cies in MC simulations, but the actual efficiencies shown in Fig. 5.5 were not used. Fig. 5.5 also 
shows the systematic uncertainty on the efficiency measurement. The systematic uncertainties 
used in this analysis were taken as the relative uncertainties shown in the plots. 

The photon energy scale and energy resolution were taken from the measurements of Z — )• ee 
events [103 mentioned in Sec. 5.1 Slight modifications were applied to the energy scale correc- 
tion to account for the differences in the interaction of electrons and photons with the material 
in front of the calorimeter and with the presampler. 
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The signature of tt-f events in the single lepton decay channel is defined by a high-i?T electron 
or a high-pT muon, large ^t, four high-px jets, and a high-Sx photon. Two of these jets are 
6-jets originating from 6-quarks. 

The background processes for tij production are the same as for tt production, but with one 
or more additional photons in the final state, which can be radiated from all charged particles. 
Moreover, events may feature hadrons and electrons misidentified as photons. 

Some of the background processes to tt production (Sec. 1.2.4) have cross sections which 
largely exceed the tt production cross section at hadron colliders. This holds in particular for 
multijet production and for W and Z production in association with jets. 

In a first step, the preselection, a pure sample of events containing top quark pairs was obtained 
by applying selection criteria which suppress different background components to tt production 
while maintaining a reasonable efficiency for tt events (Sec. 6.1). These requirements were 
motivated by previous analyses in the single lepton tt decay channel at the ATLAS experiment 
(cf. for example Ref. fl36]). In particular, the use of 6-tagging information enhanced the signal- 
over-background ratio. 



In a second step, the final event selection, photon-specific criteria were applied (Sec. 6.2) 



Sec. 6.3 presents the event yields for data, which are compared to the expectations for the signal 
and the different background contributions. 



6.1 Preselection 

Events in the electron and muon channels were required to be triggered by the EF_e20jiedium 



and EFjnulS trigger chain, respectively (Sec. 5.1 and Sec. 5.2). Non-collision background was 



rejected by selecting only events with a good primary vertex with at least five tracks associated 
to it. 

Exactly one good lepton was required in each event: in the electron channel, an electron 
with an of at least 25 GeV was required. In the muon channel, exactly one muon with a 
minimum px of 20 GeV had to be present. Events in which the muon px exceeded 150 GeV were 



disregarded (Sec. 5.2). In both channels, no lepton of the other type beyond the respective 
or px threshold was allowed to be present. Electrons were required to match the trigger object 
that fired the EF_e20_mediuni trigger: the distance in r]-(j)-spa,ce between these objects had to 
be smaller than 0.15. Muons were required to match the trigger object that fired the EPjiulS 
trigger at L2 as well as at EE level in MC simulations: a requirement of 0.15 on the distance 
between the two objects was imposed. In data, trigger matching was not applied in the muon 
channel due to a software issue at the L2 level. A systematic uncertainty was assigned to this 



feature (Sec. 12.3) 



The offline requirements on the electron £^x and the muon px of 25 GeV and 20 GeV, respec- 
tively, were sufficiently high to ensure that effects due to the trigger threshold were avoided: in 
this region, the trigger efficiencies were constant as a function of the lepton E^ or pi, respec- 
tively. Backgrounds from multijet production with jets misidentified as electrons or with muons 
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from the decay of heavy flavour hadrons were strongly suppressed by the electron and muon 
definitions presented in Ch. [5} 

In order to avoid the rare cases where a muon was also reconstructed as an electron, events in 
which a good electron and a muon shared a track were discarded. For this check, muons were 
required to fulfil all identification criteria apart from the distance requirement between muons 
and the closest jet. 



Events with jets that failed certain quality criteria (Sec. 5.3) were disregarded, as well as 
events with potentially incorrect data in the LAr calorimeters (Sec. 5.4). By these means, biases 
of the measurement and the jet reconstruction efficiency were avoided. 

Since multijet events are typically balanced in the transverse plane, the presence of large 
provides strong discrimination between multijet events and processes with from neutri- 
nos in the final state, as in M^+jets and ti production. However, may originate from the 
miscalibrations of jet energies. Moreover, the transverse mass of the W boson, as defined by 



Tot = y'2 • (pt • - • -^x - • -^j/) , l = e,^i, 

discriminates between events with true W bosons and multijet or Z+jets production. 

A of at least 35 GeV and 20 GeV was required in the electron and muon channel, respec- 
tively. The requirement was stronger in the electron channel, because multijet events feature 
more electron than muon candidate objects. Additionally, in the electron channel a transverse 
W mass of 25 GeV was required. In the muon channel the sum of and was required to 
be larger than 60 GeV, hence exploiting the correlation between and . 

All background processes feature strongly decreasing cross sections with an increasing number 
of jets in the final state. Since ti events typically possess many jets with high p^, a minimum 
of four jets with px larger than 25 GeV was required, which enhanced the total signal-over- 
background ratio. 



Finally, at least one of the jets had to be 6-tagged using the algorithm described in Sec. 5.5 
This requirement provided additional discrimination power between ti events and the back- 
ground processes, which feature only a small amount of b- and c-jets. 



Fig. 6.1 and 6.2 show a selection of kinematic distributions after the preselection in the electron 



and muon channels, respectively. The data distributions are compared to the expectations from 



MC simulations and a data-driven estimate of the multijet background (see Sec. 11.2). The 
yields and their uncertainties for the different contributing processes are described in Sec. 6.3 

In the electron (muon) channel, the electron (muon) Et (px) and r] (0) as well as the and 
r) {(/)) of all selected jets are shown. Moreover, the total number of jets, the total number of 
6-tagged jets, the ^x and m}^ are depicted. Despite a slightly overall lower expectation in the 
number of events than observed in data, well covered by the uncertainty on the expectation, 
the agreement with the observed data distributions is good in both lepton channels. This holds 
also for the slightly asymmetric distribution of the muon and the significantly lower number 
of jets at the value of (j) which corresponds to the position of the broken optical link in the LAr 
calorimeter around (j) = ~1 (Ch. [3]). 



6.2 Final event selection 

The preselection, as detailed in the previous section, improves the tt-over-background ratio and 
defined a sample dominated by ti production, ti^ candidate events were identified in this sample 
by requiring the presence of a good photon with i?x > 15 GeV. The photon i?x requirement was 
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Figure 6.1: Kinematic distributions after the preselection in the electron channel. Shown are 
the distributions from data, and the expectations from MC simulations and a data- 
driven estimate of the multijet background: electron E't and r/, jet pi and r/, the 
number of jets, the number of 6-tagged jets, and ni^ . The hashed uncertainty 
band represents the uncertainty on the expectation. In all plots, except for the 77 
distributions, the last bin includes the overflow bin. 
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Figure 6.2: Kinematic distributions after the preselection in the muon channel. Shown are the 
distributions from data, and the expectations from MC simulations and a data- 
driven estimate of the multijet background: muon and (/>, jet pT and (p, the 
number of jets, the number of 6-tagged jets, and m}^ . The hashed uncertainty 
band represents the uncertainty on the expectation. In all plots, except for the (j) 
distributions, the last bin includes the overflow bin. 
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Figure 6.3: Distribution of the true photon E'^ in tt^ simulations. The last bin includes the 
overflow bin. 



chosen as low as possible, because the distribution of the photon £'t in ify events is dropping 
strongly towards high values of -Eti as shown in Fig. |6.3[ 

In order to suppress Z+jets events with one of the electrons misidentified as a photon, the 
invariant mass of the electron and the photon in the electron channel was required to be outside 
a window of 5 GeV around the Z boson mass: m{e^) ^ [86 GeV, 96 GeV]. As discussed in 
Sec. 5^, the distance in r]-(j)-space between good photons and the closest jet had to be larger 
than 0.5. 



Fig. 6.4 and 6.5 show the photon E^^, rj and the fraction of non-converted and converted 
photons of the ti'j event candidates in the electron and muon channel, respectively. The amount 
of real ttj events within these candidates was estimated from a template fit to the p^^^'^^ distri- 
bution as described in Sec. [7) This estimate is data-driven, because simulations are not trusted 
for the description of the background rates. Consequently, no expectations from simulations are 
shown in these plots. 

It is worth pointing out that the tt^ candidates contain a sizable fraction from hadrons 
misidentified as photons. However, the E^ spectrum of the photon candidates drops towards 
large transverse momenta, as expected for true photons as well as for fake photons. More photon 
candidates are found in the very central part of the detector than at larger values of |ry|, and 
photon candidates are non-converted and converted to approximately equal parts. 



6.3 Event yields 



Tab. 6.1 shows the event yields for the preselection for data and the expectations from MC 
simulations as well as the estimation of the multijet background from data in the electron and 
muon channels, respectively. The yields before and after requiring at least one 6-tagged jet in 
the event selection are shown. 

For the uncertainties on the expectation of the different processes, only the uncertainty on the 
cross section calculations were taken into account for ti'j, tt, single top and diboson production 
as mentioned in Ch. [4| The procedure for the estimation of the multijet background is described 
in Sec. 11.2, The uncertainty on the multijet background was estimated conservatively to 100% 



if at least one 6-tagged jet was required in the events, and to 50% otherwise. 

Additional systematic uncertainties on the expectations were not taken into account, because 
this table is only intended to illustrate that the yields observed in data are well under control 
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Figure 6-4: Photon distributions for the ttj event candidates in data in the single electron 
channel. Shown are the photon and rj, and the fraction of non-converted and 
converted photons. In the E't distribution, the last bin includes the overflow bin. 
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Figure 6.5: Photon distributions for the ti^ event candidates in data in the single muon chan- 
nel. Shown are the photon E'y and ry, and the fraction of non-converted and 
converted photons. In the E't distribution, the last bin includes the overflow bin. 



for the preselection. The yields observed in data are roughly 100 events larger than expected in 
both lepton channels after the full preselection, which is well covered by the uncertainty on the 
expectation. Overall, the yields observed in data are consistent with the sum of the expectations 
within the considered uncertainties before and after the 6-tagging requirement. 



Tab. 6.2 shows the event yields for data and the expectations for signal and the different 
background sources for the final event selection. The ti^ signal expectation from the WHIZARD 
MC simulation yields 22 it 4 and 28 it 6 events in the electron and muon channels, respectively. 
This corresponds to a combined efficiency and acceptance of 0.97% and 1.27% with respect to 
the total generated signal. The efficiency of the final event selection for ti"f events was found to 
be roughly 20% with respect to the preselection (cf. Tab. 6.1 and 6.2). 

The yields for ti'^ events outside of the signal phase space, ti^ background as described in 
Sec. |4.2[ are discussed in Sec. |11.1[ The contributions from VF+jets+7, multijet+7, Z+jets+7, 
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single top+7 and diboson+7 production are discussed in Sec. 11. 2| - 11.4 The contributions 



from electrons misidentified as photons in dileptonic tt decays, Z+jets, single top and diboson 
events was estimated by measuring the misidentification rate in data (Ch. [To| ). 

Hadrons misidentified as photons can occur in all of the background processes by the presence 
of additional jets. However, a prediction for the yield cannot be obtained from simulations as 
discussed in Sec. |4.2| In order to estimate the number of hadrons misidentified as photons in 



the selected data events, a template fit to the pjpii'^^o distribution was performed (Ch. [7|. The 
final result of the fit including systematic uncertainties is presented in Ch. 13 

As a cross-check for the final fit result, the amount of fake photons from hadrons was estimated 
from photon candidates with p^rie20 ^ 3 QeV. While only 2% of true photons have p^fi'^^o values 
larger than 3 GeV, this holds true for 42% of the hadron fakes. Hence, p^^^"^^ > 3 GeV defines 
a control region (CR) largely dominated by hadron fakes. 
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Figure 6. 6: Template fits to the p'^^^'^^ distribution for photon candidates with p'^^'^'^^ > 3 GeV 
in the electron (left) and the muon channel (right) with templates for true photons 
and for hadrons misidentified as photons. The data distribution is compatible with 
being only due to hadron fakes. In both plots, the last bin includes the overflow 
bin. 



Fig. 



6.6 



shows template fits, as described in Ch. [7| to the p™'^'^^^ distribution for photon 
candidates in the CR in the electron (left) and the muon channel (right). Templates for true 
photons and for hadrons misidentified as photons were used. The data distribution is compatible 
with being only due to hadron fakes. 

The fits yield 10 it 3 and 9 it 3 hadron fake events in the CR in the electron and muon channel, 
respectively, which translates into expectations of 24 ± 8 and 22 it 7 events for the whole p^p'^^'^o 
spectrum. The uncertainty quoted is the statistical uncertainty from the fit only. The numbers 
from the final template fit (cf. Ch. 13 ) read 20 (stat.) it 3 (syst.) and 27 (stat.) it 4 (syst.) in 
the electron and muon channel, respectively. The final estimates in both channels are consistent 
with the cross-check in the region p™"^^*^ > 3 GeV within the statistical uncertainties. However, 
it needs to be pointed out that the statistical uncertainties of the final and the cross-check 
estimates are partly correlated. The final estimate is believed to be closer to the real hadron fake 
contribution, because it takes into account more information from the whole p™"'=20 gpectrum. 
The uncertainties on the other background contributions quoted in Tab. 6.2 are the total 



uncertainties including statistical and systematic uncertainties as derived in Ch. 10 and Ch. [TT| 
In data, 52 and 70 events in the electron and muon channels were identified, respectively. The 
sum of the signal and background expectations yield 60 it 10 and 67 it 10. Hence, the number 
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of observed events in data is compatible with the expectation within the uncertainties in both 
lepton channels. 
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Table 6. 1 : Event yields for data and expectations for signal and the different background con- 
tributions for 1.04 fb~^ for the preselection in both lepton channels. 
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Table 6.2: Event yields for data and expectations for signal and the different background con- 
tributions for 1.04 fb~^ for the final event selection in both lepton channels. 
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Background processes to tt'j production feature either real photons, such as VF+jets+7 pro- 
duction, electrons misidentified as photons, or hadrons misidentified as photons (hadron fakes). 
The strategy for the analysis was set up to cope in particular with the sizable background from 
hadron fakes, r-leptons may be misidentified as photons if they decay into an electron or into 
hadrons. Both cases are covered by the treatments of electrons and hadrons misidentified as 
photons. 

The treatment of the background from events with hadrons misidentified as photons needs 



particular attention, because it is not well modelled by MC simulations as discussed in Sec. 4.2 
The analysis strategy takes this into account and the amount of signal tt'j events as well as the 
amount of hadron fakes were estimated from a template fit to the photon isolation distribution 
of the selected tt'j candidates. While prompt photons are generally isolated, hadron fakes 
are typically surrounded by other particles from the fragmentation process and hence a good 
discrimination between prompt photons and hadron fakes can be achieved by the use of isolation 
observables. 

Backgrounds with real photons or electrons misidentified as photons cannot be distinguished 
from tij events by considering the photon isolation, and they were estimated separately: the 
probability for electrons to be misidentified as photons was measured in data using a sample 



largely enriched in Z — )• ee events (Ch. 10). Background contributions with real photons were 



estimated partly from control regions in data and partly from MC simulations (Ch. 11). 



Photon isolation 

Isolation observables are constructed from the surrounding energy in the calorimeter or from 
the tracks close to the photon candidate. The transverse isolation energy in the calorimeter 
in a narrow cone around the photon depends on the photon rj, because of the varying amount 
of material in front of the calorimeter. Given the limited amount of tij candidate events, 
calorimeter isolation was hence disfavoured as a discriminating variable with respect to the 
pcone20 observable, which is independent of the photon t] to first order (Ch.|8|. p™°s20 jg (jgfined 
as the sum of the transverse momenta of the tracks in a cone of AR = 0.2 around the photon 
A small cone size of 0.2 was chosen in order to avoid a bias to p^j?"*^^^ due to 



149 



candidate 

the presence of nearby particles from the tt'j final state. 

Tracks were required to have a of at least 1 GeV, at least seven hits in the Pixel and SCT 
detectors, and a hit in the Pixel 6-layer. Tracks associated to conversion vertices closer than 
0.1 to the photon in ri-(p-space were not considered. The transverse and longitudinal impact 
parameters of the tracks with respect to the primary vertex had to be smaller than 1 mm to 
reduce biases from tracks that originated from pile- up interactions. The longitudinal impact 



parameter requirement was not included in the definition of p^^^^^o £qj, photons 149 , although 
it was used for electrons. This requirement was added for photons for a consistent treatment of 
electron and photon isolation, given that the photon distribution was estimated using electrons 
from Z — )• ee decays (Ch. In order to add the requirement on the longitudinal impact 
parameter, the photon observable was recalculated. 
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7 Analysis strategy 



Fig. |7.l| shows the p™°s20 (distributions normaUsed to unity (templates) for prompt photons 
and for hadrons misidentified as photons. The derivation of the templates is described in detail 
in Ch. |8] and |9j The bin sizes were chosen in studies using simulations such that the expected 
bin contents for the hadron fake distribution were similar. 
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Figure 7.1: p^^^'^^ templates for prompt photons and fake photons from hadrons, normalised 
to unity. The last bin includes the overflow bin. Details about the derivation of 
the templates can be found in Ch. [8] and [9j 



Description of the template fit 

The 

^cone20 (distribution of the photons from the selected tt'j candidate events was fitted using 
template distributions of real photons and hadron fakes in order to estimate the expected number 



of signal events (s). The template distributions are those presented in Fig. 7.1 A binned 
likelihood fit was performed in five bins: [0 GeV, 1 GeV), [1 GeV,3 GeV), [3 GeV, 5 GeV), 
[5 GeV, 10 GeV) and [10 GeV, oo). 
The number of events in each bin i of the signal template distribution relates to s by 

Si = ei-s, (7.1) 

where describes the acceptance and selection efficiency for ti^ events, and the probability 
to end up in bin i. For each background j, templates were built to describe their respective 
contribution 6^ in the i-th bin of the isolation distribution. Hence, the sum of all contributions 
in each bin reads: 

iVbkg 

= + ^ 6i . (7.2) 

i=i 

The following likelihood was then maximised in the fit using Markov Chain Monte Carlo 



implemented in the Bayesian Analysis Toolkit 154 



A^bins ^bkg 

L = n P(iV,|A,) • n ni^) ■ Pi^) > (7-3) 

i=l j=l 

where Ni is the number of observed events in bin i of the isolation distribution. P(Ni\Xi) is 

the Poissonian probability to observe Ni data events given an expectation of Aj. P{V) is the 
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probability for the j-th background contribution, and P{s) is the probability for the signal 
contribution. 

The background probabilities were either chosen to be constant in a range [b^^^^, 6max], if the 
background yield was treated free parameter 



0, 



else 



or fixed to a background estimate V: 

P{h>) = 5{V -V) , 

where 6{x) is the delta distribution. The uncertainty on the background estimate V was then 



treated as a source of systematic uncertainty (Ch. 12.2) 



Tab. 7.1 gives an overview of the different parameters of the template fit and their respective 
probabilities: as already mentioned, the hadron fake contribution was treated as a free parame- 
ter, and a constant background probability was assigned to it covering the whole range of hadron 
fake contributions between 0% and 100%. 

The tt'j signal contribution was also treated as a free parameter: 



Pis) 



0, 



else 



with s covering a range of signal fractions between 0% and 100%. It is worth stressing that 
the parameter describing the number of signal events accounts for the acceptance and selection 
efficiency and therefore represents the total number of ttj events extrapolated to the whole 
signal phase space. Acceptance and efficiency were estimated in MC simulations, and systematic 



uncertainties were evaluated accordingly (Sec. 12.1 and 12.3). 



The contributions from background processes with electrons misidentified as photons and 
from processes with true photons were fixed to the estimates derived in Sec. 10. 2| and Ch. [TT| 
respectively. Systematic uncertainties were evaluated by up- and down- variations of the different 



contributing processes (Sec. 12.2) 



Process 


Parameter 


Parameter 
probability 


Note 


ti^ signal 


free 


constant 


accounts for acceptance and selec- 
tion efficiency 


True photon bkg. 


fixed to estimate 


delta 
function 


systematic uncertainties by up- 
and down-variation 


Electron fakes 


fixed to estimate 


delta 
function 


systematic uncertainties by up- 
and down-variation 


Hadron fakes 


free 


constant 





Table 7. 1 : Parameters and parameter probabilities for the signal and background contributions 
used in the template fit. 
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7 Analysis strategy 



The template fit was performed in botli lepton ctiannels simultaneously. One combined like- 
lihood was constructed in order to estimate the expected number of signal events s, which, 
combining Eq. (7.1), (7.2) and (7.3), explicitly reads: 



L = 11 P\ iV*,c+jcts 
i=l 

1=1 



A^bkg 



Ai,e+jets — Ei.e+jets "Si ^ ^j^e+jets 

iVbkg 



9 ■ 

e+jets 



A^bkg 
A'bkg 



■Pis) 



As a cross-check, also separate fits in the electron and muon channels were performed. 



56 



8 Derivation of the prompt photon template 



In order to reduce the dependence on MC simulations, the template distributions for prompt 
photons were derived using a data-driven approach. Since the definition of a data sample which 
contains mainly prompt photons is not straightforward, the prompt photon template was derived 
from electrons from Z — t- e^e~ data, because the isolation properties of electrons and photons 
are very similar given the similar detector signature. Small differences in the pSj?"^^*^ distribution 
between photons and electrons were corrected for using simulated Z — )• e^e~ and ti7 events. 

A data sample dominated by Z — )• e^e~ events was selected by the following criteria: two 
oppositely charged electrons had to be present in the event, and their invariant mass had to be 
within a 50 GeV window around the mass of the Z boson. Since the background from multijet 
events is small close to the Z mass peak, the number of selected events was increased by loosening 



the requirements on the electron objects with respect to the definitions presented in Sec. 5.1 
only medium shower shape criteria were applied, which still lead to a high signal-over-background 
ratio. In order to avoid a bias on the ■p^'^^'^^ distribution from the calorimeter isolation criterion, 
the latter was disregarded. 

Events were triggered by the EF_e20_mediuin trigger, which is the same trigger as used for 
the selection of the it7 event candidates. Consequently, the first electron was required to have 
an E'j^ of at least 25 GeV. The E't threshold for the second electron, however, was lowered to 
15 GeV, which corresponds to the lower threshold for the photon E't in this analysis. 

Selected events were required to have a good primary vertex with at least five tracks associated 
to it. Events with large noise in the LAr calorimeter were rejected. 



|Ldt = 1.04fb' 



opposite sign 
ft □ same sign subtracted 



S 80000 ^ 

■g 70000 - ATLAS work in progress 

S 60000 - 
w - 

50000 ^ 

40000^ a a * ™S = 6.00 GeV J 

30000- ° RMS = 5.79 GeV ^ 

a ^ 

i ■ : 



20000 : 
10000: 



70 75 80 85 90 95 100 105 110 115 

rriee [GeV] 



Figure 8.1: Distribution of the di-electron invariant mass in data. Electrons fulfil the medium 
shower shape criteria. The distribution for opposite-sign pairs is shown with the 
solid circles. When the contribution from same-sign electron pairs was subtracted, 
the distribution with the open squares was obtained. The first bin includes the 
underflow bin, and the last bin includes the overflow bin. 



1.04 fb ^. The solid circles show the distribution for electron pairs with opposite charge. A 



Fig. 8.1 shows the distribution of the di-electron invariant mass in the whole data set of 
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8 Derivation of the prompt photon template 



small contribution from multijet events with two jets misidentified as electrons was accounted 
for by subtracting the contribution from electron pairs with the same electric charge, as shown 
by the distribution with the open squares. Since no real electron is present in multijet events, the 
measurement of the electric charge is arbitrary and pairs with same and opposite charge occur 
equally often. Charge misidentification in Z ^ e~^e~ events also contributes to the opposite-sign 
distribution, but does not introduce a bias to the pJj?'^^^^ distribution of the electrons. Hence, 
subtracting the opposite-sign contribution removes the background from multijet events from 
the same-sign p^^^"^^ distribution. 

When the contribution from same-sign electron pairs is subtracted, the width of the invari- 



ant mass distribution decreases (Fig. 8.1): the root mean square is reduced from 6.00 GeV to 
5.79 GeV. The fact that the resolution of the Z boson mass peak was found to be improved 
after the subtraction of the same-sign contribution indicates that the background from multijet 
events was indeed suppressed. 

The p^'^^'^^ distributions for the electrons from this data sample are shown in Fig. 
different values of (left plot) and different regions in |r/| (right plot): [15 GeV, 20 



8.2 



for 



GeV), 

[20 GeV, 30 GeV), [30 GeV, 50 GeV), [50 GeV, 100 GeV), and [0,0.60), [0.60,1.37), [1.52,1.81), 
[1.81, 2.37), respectively. The upper part of the plots show the regions 15 GeV < E't < 20 GeV 
and < rj < 0.6, respectively. The lower part of the plots show the difference with respect to 
this distribution for the other regions in E't and rj. The observable is stable with respect 

to the electron and rj. Discrepancies between different Et and |?7| bins are smaller than 2% 



and were ignored given the size of the systematic uncertainties discussed in Ch. 12 
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Figure 8.2: p'^'^^'^^ distributions for electrons from Z ^ e~^e~ decays for different regions 
of Et (left) and |??| (right). The upper part of the plots show the regions 
15 GeV < -Et < 20 GeV and < rj < 0.6, respectively. The lower part of the plots 
show the difference with respect to this distribution for the other regions in E'j' 
and rj. In both plots, the last bin includes the overflow bin. 



Differences between the p^^'^'^^ distributions of electrons and photons were studied in simulated 
events. The same event selection as for the selection of the Z — )• e~^e~ data events was applied to 
Z — )• e~^e~ simulations in order to obtain electron distributions. These distributions were then 
compared to those from real photons from simulated tt'j events. 

Fig. 8.3 and 8.4 show the electron distribution in the upper part of each plot. The lower 
part of each plot shows the difference of the distribution of photons from simulated tt'j events 
(solid line) with respect to the electron distribution. Additionally, distributions for unconverted 
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Figure 8.3: pJj?'^^^'^ distributions for electrons from simulated Z — )• e^e~ decays (upper part 
of each plot) in different bins of £Jt for < Iryl < 0.60 (four upper plots) and 
0.60 < |?7| < 1.37 (four lower plots) normalised to unity. The lower part of each 
plot shows the difference of the distribution of photons from simulated ti'^ events 
(solid line) with respect to the electron distribution. Additionally, the distributions 
for unconverted (dotted grey line) and converted photons (dashed black line) from 
tfy simulations are depicted. In all plots, the last bin includes the overflow bin. 
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8 Derivation of the prompt photon template 



;o.5 

3 
) 
j 



, 0.1 

) 

^-0.1 



ATLAS work in progress 
simulation — 



1 .52 < hi < 1 .81 
15GeV<Ej<20GeV - 

■ e(Z ^ ee) 

■ y (tty) 

unconverted y (tty) 
converted y (tti) 



2 4 6 



10 12 14 16 18 20 

pcone20 [QgV] 



:0.5 





, O.lh 



2 



'-0.1 





ATLAS work in progress 


1 .52 < |ti| < 1 .81 

30 GeV < < 50 GeV - 




simuiation 


e(Z ^ ee) 

y (fty) - 
unconverted y (tty) - 
converted y (tti) 




' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 


' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' 








J 


^ 











2 4 6 



10 12 14 16 18 20 

pcon820 [QgV] 



iO.5 



, 


'-0.1 



ATLAS work in progress 



simuiation 



1.81 < hi < 2.37 
15GeV<Ej<20 GeV 

■ e(Z ^ ee) 

■ y (tty) 

unconverted y (tty) 
converted y (tty) 



_ ^iffi ililifc 



2 4 6 

ATLAS work in progress 
simuiation — 



10 12 14 16 18 20 



;o.5- 



0- 

0, 0.1- 
o 

£ 0- 



1.81 < hi < 2.37 

30 GeV < < 50 GeV 

■ e(Z ^ ee) 

■ y (t'ty) 

unconverted y (tty) 
converted y (tty) 



M ill 



'-0.1- 



2 4 6 



10 12 14 16 18 20 

pcons20 [QgV] 



= 0.5 



ATLAS work in progress 



simuiation 



1.52 < hi < 1.81 
20 GeV < < 30 GeV 

e(Z^ ee) 

— y (tty) 
unconverted y (tty) 

— converted y (tty) 



10.5 


> 0.1 

\ ° 

^-0.1 



2 4 6 8 10 12 14 16 18 20 

pcone20 [QgV] 

._, , . 1.52 < hi < 1.81 

ATLAS work m progress gp QeV < < 1 00 GeV- 

e(Z^ ee) 

simuiation y (tty) _ T 

unconverted y (tty) - 
— converted y (tty) " 



10 12 14 16 18 20 

pcone20 [QgV] 



= 0.5 
) 



> 0.1 

\ ° 

^-0.1 



ATLAS work in progress 



simuiation 



1.81 < hi < 2.37 
20 GeV < E^ < 30 GeV 

e(Z^ ee) 

— y (tty) 
unconverted y (tty) 

— converted y (tty) 



1^ 



= 0.5 


> 0.1 

\ ° 

^-0.1 



ATLAS work in progress 



10 12 14 16 18 20 

pcone20 [QgV] 

1.81 < hi < 2.37 
50 GeV<E^<100 GeV- 

e(Z^ ee) 

simuiation y (tty) _ ~ 

unconverted y (tty) - 
— converted y (tty) i 



10 12 14 16 18 20 

pcone20 [QgV] 



Figure 8.4: p^j?'^^^^ distributions for electrons from simulated Z — )• e~^e~ decays (upper part 
of each plot) in different bins of for 1.52 < \ri\ < 1.81 (four upper plots) and 
1.81 < |?7| < 2.37 (four lower plots) normalised to unity. The lower part of each 
plot shows the difference of the distribution of photons from simulated tt'j events 
(solid line) with respect to the electron distribution. Additionally, the distributions 
for unconverted (dotted grey line) and converted photons (dashed black line) from 
tt'f simulations are depicted. In all plots, the last bin includes the overflow bin. 
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(dotted grey line) and converted photons (dashed black line) are shown in four different regions 
of and \r]\. 

The differences between the shapes of the electron and photon p™"'^^'^ distributions are small 
for low E't and increase up to 0.1 in the first two bins with increasing electron and photon E't 
in the different |r/|-regions. The |r/|-region [1.81,2.37) is special, because the acceptance of the 
TRT ends at \r]\ < 2.0, and the differences between electron and photon p^^<^'^^ are smaller than 
in the other |r/|-regions. No significant differences between unconverted and converted photons 
were observed, so that both photon types were treated together. 

A priori no differences between electrons and photons are expected for the pSj?"*^^*^ distributions. 
If there were inherent differences between electrons and photons, they would unlikely contribute 
to the differences observed in the comparison of Z ^ e~^e~ and ti^ simulations: one-track 
conversion photons are by construction very similar to electrons, and the differences between 
electrons and photons are significantly larger than those between the p^^'^'^^ distributions of 
unconverted and converted photons. Hence, the difference in the topologies of Z — )• e~^e~ and 
ti'y events is expected to cause the observed electron-photon differences. 

shows the p^'^^'^^ distributions for simulated electrons and photons, as already shown 
id 8.4 However, only photons that had a minimum distance of 0.2 in r]-(j)-space 
from a true electron or muon were considered. Such leptons only needed to be present among 
the generated particles before detector simulation, but did not need to fulfil all quality criteria 



Fig. 
in Fig 



8.5 



of the lepton definitions (Sec. 5.1 and 5.2) and did not even need to be reconstructed at all. The 



track of such leptons, however, is likely to have a large contribution to the p^'^^'^'^ value of the 
close-by photon. This effect is particularly relevant for dileptonic tt decays. 
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Figure 8.5: p^^^'^^ distributions for electrons from simulated Z — ?■ e'^e decays (upper part of 
the plots) normahsed to unity for 0.60 < \7]\ < 1.37 and 30 GeV < Et < 50 GeV 
(left plot) as weh as for 50 GeV < < 100 GeV (right plots). The lower part 
of each plot shows the difference of the distribution of photons from simulated 
tij events (solid line) with respect to the electron distribution. Additionally, the 
distributions for unconverted (dotted grey line) and converted photons (dashed 
black line) from ttj simulations are depicted. In both plots, the last bin includes 
the overflow bin. 

Photons closer than 0.2 in r/-(/)-space to a true lepton were not considered and the 
agreement with the electron distribution is improved with respect to the plots in 
Fig. [Ol and [8^ 
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8 Derivation of the prompt photon template 



Only two example plots for the region 0.60 < |?7| < 1.37 for transverse momenta between 30 
and 50 GeV (left plot) and 50 and 100 GeV (right plot) are shown. For reference, the remaining 
plots are presented in App. [Aj Compared to Fig. |8.3[ the agreement between the electron and 
photon distributions is significantly improved. In particular, the increased tail in the photon 



distributions, as seen in the higher bins in Fig. 8.3 and Fig. 8.4 is suppressed. The differences 
between the electron and photon distributions were hence concluded to be mainly due to the 
different topologies of Z ^ e^e~ and tt'j events. 

The differences between the electron and photon distributions in the MC simulations were 
used to derive a prompt photon template from the electron template Sg measured in data: 



s^ = Se + As^^{ET,\v\) with As^^ 



MC 



The derivation was done in 16 bins, four in E'y and four in \rj\, as detailed above. The rela- 
tive weights of the different extrapolation bins were taken from the i?T and \ri\ distributions of 
real photons in tt'j MC simulation, presented in Fig. 8.6 The resulting prompt photon tem- 



plate was presented in Fig. 7.1 The systematic uncertainty assigned to this electron-to-photon 



extrapolation is discussed in Sec. 12.2 
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Figure 8.6: Distribution of the i?T and |r/| of real photons in ^7 MC simulation in 16 bins, as 
used for the derivation of the photon p^'^^'^^ templates from the electron templates 
extracted from Z — )• e~^e~ data. 
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9 Background events with hadrons 
misidentified as photons 



As discussed in Sec. 4.2 and Ch. [7j simulations cannot be trusted for the modelling of the 
properties of hadrons misidentified as photons, and the contribution from hadron fakes was 
hence estimated using a template fit to the p^^°'^^ distribution of the photon candidates in data. 
The 

^cone20 template distribution for the hadron fakes was derived from a control region (CR) 
in data largely enriched in hadron fakes, as described in this chapter. 

The CR was defined by requiring that at least one out of the following four shower shape 
variables from the LAr strips failed the cut from the tight photon menu: AE, inside! and 



-Erratic (cf. Tab. 5.1). All other criteria of the good photon definition detailed in Sec. 5.6 needed 



to be satisfied. Such photon objects are called hadron fake candidates in the following. 

Variables constructed from the LAr strips were chosen, because they are sensitive to the core 
of the photon cluster and are hence expected not to be correlated to photon isolation variables, 
such as pSj?"^^*^. This procedure was previously used in the prompt photon 
cross section measurements 



153 



and diphoton 



151], where the CR was shown to model the properties of hadron 
fakes well, although these analyses used calorimeter isolation instead of p^^^'^^. 
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Figure 9. 1 : Normalised p^^^"^^ distributions for true hadron fakes (solid line) and hadron fake 
candidates (dashed line) from tt MC simulation. The last bin includes the overflow 
bin. 



It was therefore checked that also the p^^'^'^^ distribution in the tt environment is well modelled 
by the choice of this CR. Fig. 9.1 shows the p^^^'^^ distribution for true hadron fakes (solid line) 



and hadron fake candidates (dashed line) from tt MC simulation. The distributions for true 
hadron fakes and the background model agree within uncertainties. The four shower shapes 
used for the deflnition of the CR have hence been shown to be sufficiently uncorrelated with the 
^cone20 variable, and the hadron fake candidate selection therefore provides a good estimator for 
the p^'^^'^^ distribution of hadron fakes in the tt environment. 
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9 Background events with hadrons misidentiBed as photons 
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Figure 9. 2: Normalised p^'^^ distributions for 
region for different regions of -Et (1* 
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In order to derive the p^'^'^^o (Jigtributions from data, events from the skim of the JetTauEtmiss 
stream (Ch. [3| were used. Fig. 9.2 shows the p'^'^'^'^^ distributions for hadron fake candidates in 



different regions of and \ri\ as derived from data. In contrast to the electron distributions, 
the distributions for hadron fake candidates show a clear dependence on the i?T of the object. 
This behaviour is expected, because hadrons with large i?T are likely to originate from the jet 
fragmentation of high-pT partons, which features more track activity than the fragmentation of 
lower-pT partons. 

Two regions in can be distinguished: very similar p'^'^'^'^^ distributions are observed for 
photons in the region < |r/| < 1.81, while they differ significantly from the p5j?°s20 (Jigtribution 
in the region 1.81 < \r}\ < 2.37. This is explained by the acceptance limit of the TRT at |r/| = 2.0. 

In order to correctly construct the fake hadron template, the templates for the different E't 
and r/-regions needed to be reweighted according to the £^t and r] spectra of the hadron fakes 
within the selected ti^ candidates. In the following, the procedure for the estimation of the 
E'Y spectrum as well as of the fraction of photons in the high- |r/| -region 1.81 < \ri\ < 2.37 is 
described. 

First, it is shown that the E't spectra of fake photons are very similar in the low- and the high- 
|7y|-regions. Fig. 9.3 shows the photon E't distributions from a full event selection in tt simulations 
for true hadron fakes and hadron fake candidates. In both cases, the £Jt distributions are 
similar in the two |r/| -regions. Hence, the reweighting of the p'^'^^'^^ templates can be performed 
independently in and ij. 

The idea for the estimation of the fake photon i?T spectrum and the high-lryl fraction is the 
following: the good photon definition is replaced by the hadron fake candidate definition in the 
event selection. The resulting photon E't and rj spectra in data are then used to estimate the 
distributions for hadron fakes passing the actual good photon definition. 

This approach was validated using tt MC simulations: the left plot of Fig. 9.4 shows the Ej' 
spectra for hadron fakes and hadron fake candidates from a full event selection. The right plot 
shows the fraction of photons with < |r/| < 1.81 and 1.81 < \r]\ < 2.37, respectively. 

The Et spectrum of the hadron fakes is modelled by an exponential function and the left plot 
of Fig. 9.4 shows the exponential fits to the photon E^ spectra: the exponential model is fitting 
well the MC distributions. The mean lifetimes extracted from the fits read 27.1 it 0.9 GeV and 
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Figure 9.3: £'t spectra of true hadron fakes (left) and hadron fake candidates (right) in ti 
simulations for < jr^l < 1.81 (solid line) and 1.81 < |r/| < 2.37 (dashed line). In 
both plots, the last bin includes the overflow bin. 
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Figure 9.4: The left plot shows the E't spectra of true hadron fakes (solid line) and hadron 
fake candidates (dashed line) in tt simulations. An exponential fit to both spectra 
is also depicted. The right plot shows the fraction of photons with \ri\ smaller than 
1.81 and larger than 1.81. In the E't distribution, the last bin includes the overflow 
bin. 



21.7 lb 0.7 GeV for true hadron fakes and hadron fake candidates, respectively. The difference of 
25% is taken into account as a systematic uncertainty on the exponential estimated from data. 

The event selection with the hadron fake candidate definition yields 17 events in data in the 
single electron channel and 26 events in the single muon channel. Since no differences between 
the photon spectra in the two lepton channels are expected, the events in both channels were 
combined in order to increase the available statistics. The left plot of Fig. |9.5| shows the 
spectrum of the hadron fake candidates from this selection. The exponential fit is also shown. 
The entries in the last bin, which contains the overflow bin, are neglected. The fit yields a 
mean lifetime of 31 it 11 GeV, which is compatible within uncertainties with the mean lifetime 
observed in MC simulations. Considering an additional 25% uncertainty as discussed above, 
results in an estimate of 31 it 13 GeV. Due to the limited statistics, the uncertainty on the 
mean lifetime of the exponential is large. This uncertainty is considered in the evaluation of the 



systematic uncertainties in Sec. 12.2 
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9 Background events with hadrons misidentiBed as photons 
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Figure 9.5: The left plot shows the estimated spectrum of hadron fake candidates in data 
in both lepton channels together. An exponential fit to the spectrum is also shown. 
The right plot shows the fraction of photons with |r/| smaller than 1.81 and larger 
than 1.81. In the E't distribution, the last bin includes the overflow bin. 



The effect of the last bin, which contains the overflow bin, on the estimation of the -Bt spec- 
trum is negligible: taking the last bin into account did not change the result of the exponential 
fit significantly. Even when the weight for the high-£^T-region was enhanced by 0.07, which 
corresponds to the three events in the last bin, this was found to be covered by the systematic 
variations of the i?T spectrum, for which the weight of the high-i?T-region was varied by +0.12 
and —0.19, respectively. 



The right plot of Fig. 9.5 shows the fraction of hadron fake candidates with < |?7| < 1.81 and 



with 1.81 < \r]\ < 2.37. The estimated fractions are in agreement with the MC estimates shown 



in Fig. 9.4 within the statistical uncertainties. 

The model for the hadron fake spectrum and high-|r/| fraction was used to reweight the 
eight data templates b^^'^ (p^j?"^^^ | E^^, \r]\^ from the four different regions in E't and the two 
regions in |r/| to one single template for the background from hadrons b^^'^ (p™"'^^'^): 

^had (pconc20) = ^ Wexp (E^) ■ Wfrac {\v\) ' t"^"^ {pT"^' I ^T, H ) , 

Et bins I'*?! bins 

where Wgxp (-Et) and Wfj-ac (|??|) are the weights derived from the exponential fit to the -Et 
spectrum and the high-|r/| fraction for the given bin in E't and Ir^l, respectively. The resulting 



template b^'^'^ was presented in Fig. 7.1 
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10 Background events with electrons 
mis identified as photons 



Electron and photon objects are very similar: since their electromagnetic clusters have nearly 
identical properties, electrons and photons are only distinguished by the tracks associated to 
the cluster. Unconverted photons do not feature associated tracks, and converted photons are 
discriminated from electrons by requiring that the associated tracks belong to a conversion vertex 
(Sec. [5^. 

However, electrons may be misidentified as photon objects when either the electron track was 
only reconstructed with poor quality or not at all. Electromagnetic clusters from electrons may 
also feature additional tracks pointing to them, which may originate from close-by jet activity, 
from the underlying event, or from another collision in the same bunch crossing (pile- up). Such 
an electron may be misidentified as a converted photon. 

The electron-to-photon misidentification rate /e-s>7 is defined as the probability for a true 
electron to be misidentified as a photon object. /e-s>7 was measured in data using Z — )• ee events 
using a tag-and-probe method as described in Sec. 10.1 Scale factors (SFs) with respect to the 
rate observed in MC simulations were derived in 16 bins of the photon and rj. The main 
backgrounds with misidentified electrons are dileptonic tt decays and Z+jets production, but 
also single top and diboson production contribute. MC simulations of these backgrounds were 



corrected using the SFs as described in Sec. 10.2 



10.1 Estimate of the misidentification rate 

Z ^ ee decays provide a clean environment for the study of electron properties. The mea- 
surement of /e-s>7 was based on the ratio of Z — t- ee and Z — t- e7fake events in data, where in 
Z — )• e7fake one electron was misidentified as a photon. In Z — )• e7fake events, the electron was 
required to have a larger E't than the photon. Therefore, the electron with the larger in 
Z — 7- ee events and the electron in Z — >• e7fake) called tag electron in the following, originated in 
most cases from the higher- E't electron in Z — t- ee events. 

Let £i denote the combined trigger, reconstruction and identification efficiency for the tag 
electron. Let £2 denote the combined reconstruction and identification efficiency for the second 
electron in the Z — )■ ee selection. Then, the numbers of events in the Z ^ ee sample (A'^ee) and 
in the Z e7fake sample (N^'y) are given by 

Nee = El ■ 82 ■ N , and 

-^67 = £1 • /e-s>7 • N , 

where is the total number of true Z — t- ee events in the acceptance region. 
Hence, /e-s>7 is given by 

_Ne^ 
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10 Background events with electrons niisidentiBed as photons 
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Figure 10.1: Misidentification rate fe^'y from MC simulations as a function of the true electron 
rj (left) and i?T (right). 



Fig. 10.1 shows /e-s>7 as a function of the true electron 77 (left) and i?T (right) in simulations. 
The dependence on r] is explained by the varying granularity of the calorimeter and the changing 
amount of material in front of it. /e-!>7 may be affected by the altering shape of calorimeter 
showers and the varying performance of the track reconstruction. While for high electron E't 
the misidentification rate is constant at roughly 6%, for low E't, a turn-on is observed, because 
in this analysis photons were required to have a minimum E'y of 15 GeV. 

A scale factor was derived in order to correct simulations so that they describe the misiden- 
tification rate observed in data: 



SF 



Ne 



data _ ^QJ^ 



MC 

£2 cancels out, because the simulations were corrected for discrepancies with respect to data 
using SFs for electron reconstruction and identification (Sec. 5.1). 

For the measurement of the misidentification rate, events triggered by the EF_e20jiedium 
chain were selected which featured also a good vertex with at least five associated tracks. At 
least one good electron (Sec. |5.1[ ) with E't > 25 GeV was required, which was closer than 0.15 
in 77-i;^-space to a trigger object (tag electron). Only events with good data taking conditions of 
the LAr calorimeters were considered. 

For the selection of Z — )• ee events, a second good electron with E't > 15 GeV was required 
{probe electron). In order to avoid a bias from the trigger, the first electron had to fire the trigger, 
and the second electron was used for the measurement of the misidentification rate, where the 
second electron was chosen to be the electron with the lower Eic ■ The invariant mass of the two 
electrons was required to lie within a ±50 GeV window around the Z boson mass. This wide 
mass range was chosen in order to include side bands with large background contributions from 
multijet production. 

For the selection of Z ^ e7fake events, a photon with Ei^ > 15 GeV was required {probe pho- 
ton). The Et of the photon had to be lower than the -Et of the tag electron. The invariant 
mass of the electron and the photon also had to lie within a ±50 GeV window around the Z 
boson mass. 



The left plot in Fig. 10.2 shows the invariant mass distributions from 1.04 fb ^ of data for the 



two electrons from the Z — ?• ee selection. The right plot shows the invariant mass distribution 



68 



10.1 Estimate of the misidentification rate 



for the electron and the photon from the Z e^f^ke selection. Both distributions feature a 
distinct Z boson mass peak, which indicates that the photon candidates from the Z e7fake 
selection were indeed dominated by misidentified electrons from the Z boson decay. 
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Figure 10.2: Invariant mass distributions for the two electrons from the Z — )• ee selection (left), 
and the electron and the photon from the Z — ?• e7fake selection (right) in data. 



Estimation of the background contribution from multijet events 

Hadrons from jet fragmentation produced in multijet events may be misidentified as electrons 
or photons. Since the cross section for multijet production is much larger than the cross section 
for Z boson production, multijet events with two misidentified hadrons may contribute sizeably 
to the Z — )• ee and Z — )• e7fake samples. 

A broad invariant mass window around the Z boson mass was chosen in order to include 
regions of very high and low invariant masses, which feature a larger fraction of multijet events 
than the central peak. These side bands were used to estimate the contribution from multijet 
events in a narrow window of 10 GeV around the Z mass. The purer Z — )■ ee and Z — )• e7fake 
samples in the narrow mass window were then used for the measurement of the misidentification 
rate /e_^^. 

In order to estimate the multijet background, the invariant mass distributions presented in 



Fig. 10.2 were fitted with a signal and a background model: while the multijet background was 



modelled with an exponential, the signal Z peak was modelled with the convolution CB*BW 
of a Breit-Wigner distribution (BW) for the Z boson mass and a Crystal-Ball function CB in 
order to describe the mass resolution: 



CB (m; a, n, m, a) = norm. • < 



(m— m)^ 



< -a 



else 



Fig. 10.3 shows the fits to the invariant mass distributions for the two electrons from the 
Z — 7- ee selection (left) and for the electron and the photon from the Z — )• e7fake selection (right) 
in data. The combination of the background fit (dashed line) and the signal fit (solid line) fits the 
data very well. The fit range was chosen to [65 GeV, 135 GeV]. In order to allow for a variation 
of the limits for the evaluation of the systematic uncertainty on the fit stability, the upper limit 
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10 Background events with electrons niisidentiBed as photons 
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Figure 10.3: Invariant mass distributions for the two electrons from the Z — )• ee selection (left) 
and for the electron and the photon from the Z — )• e7fake selection (right) in data. 
The distributions were fitted with a convolution of a Breit-Wigner distribution 
and a Crystal-Ball function for the signal contribution (solid line), and an expo- 
nential for the background modelling (dashed line). 



was chosen to be slightly lower than the maximum value used in the selection (141 GeV). In 
the low-mass region (< 60 GeV), however, the signal distribution does not follow the CB*BW 
shape, because of the minimum requirements on the electron and photon E-^ of 25 GeV and 
15 GeV, respectively. In order to avoid a bias due to this effect, the lower bound was chosen at 
65 GeV. 

While the fit in the Z ^ ee sample yields a contribution of 0.97 it 0.04 % from multijet events, 
the fit in the Z — )• e7fake sample yields a multijet contribution of 6.5 ± 0.3 %. The uncertainties 
are the statistical ones from the fit only. As expected, the multijet contribution is larger in the 
Z — )• e7fake sample, because no isolation criterion was applied to photons. 

The fraction of multijet events is expected to decrease with the object's £'t, because the 
spectrum for jets drops exponentially, while the i^T spectrum for electrons from Z boson decays 
drops less strongly due to the underlying Breit-Wigner distribution. In addition, a dependence on 
r/ can be expected, because the varying detector geometry changes the probability to misidentify 
objects from jet fragmentation as electrons. Hence, the multijet background in the Z — )• e7fake 
sample was estimated in four bins of i?T and four bins of |??| of the probe photon: the bins in |r/| 
were the same as used in Ch.|8]and[9| which read [0,0.6), [0.6, 1.37), [1.52, 1.81), and [1.81, 2.37). 
The following bins were chosen: [15 GeV, 20 GeV), [20 GeV, 30 GeV), [30 GeV, 40 GeV), and 
[40 GeV, 45 GeV). Events where the photon had an £'t of more than 45 GeV were not taken 
into account, because in these events also the tag electron had an E^^ of at least 45 GeV. Z 
decays in which both decay products have a very large i?T are rare and also feature an invariant 
mass distribution which differs from the form of the CB*BW signal model. 

Fig. 10.4 shows the fits for the Z — )• e7fake selection for photons in the region < I??! < 0.6 and 
15 GeV <Et <2Q GeV (left), and in the region 0.6 < < 1.37 and 30 GeV < E^ < GeV 
(right). The resulting fraction of multijet events in the Z — ?• e7fake sample are summarised in 
Tab. |10.1 including the statistical uncertainty from the fit. 

Tab. 10.1 also shows systematic uncertainties on the multijet fraction. They were estimated 
by comparing the result from the fit as described above with a fit where the signal model was 
replaced by template distributions retrieved from Z — )• ee simulations. Both, the tag electron 
and the probe photon, were required to be matched to one of the generated electrons from 
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Figure 10.4: Invariant mass distributions for the electron and the photon from the 
Z — )• e7fake selection in data for photons in the region < \ri\ < 0.6 
and 15 GeV < E^t < 20 GeV (left), and in the region 0.6 < |r?| < 1.37 and 
30 GeV < -Et < 40 GeV (right). The distributions were fitted with a convolution 
of a Breit-Wigner distribution and a Crystal-Ball function for the signal contri- 
bution (solid line), and an exponential for the background modelling (dashed 
line) . 




Figure 10.5: Invariant mass distributions for the electron and the photon from the 
Z — )• e7fake selection in data for photons in the region < \ri\ < 0.6 
and 15 GeV < E^t < 20 GeV (left), and in the region 0.6 < |7?| < 1.37 and 
30 GeV < Et < 40 GeV (right). The distributions were fitted with templates 
derived from MC simulations (solid line), and an exponential for the background 
modelling (dashed line). 



the Z decay within a cone of 0.2 in rj-4>-space. Fig. |10.5| shows two example fits in the same 
regions as shown in Fig. 10.4 Generally, the estimated multijet fractions are of the same order 
in both approaches, but sizable differences of up to 12% were observed in the low-£^T-region. 
The difference was symmetrised and taken as a systematic uncertainty. 

Additionally, the upper and lower limits of the fit range were varied by it 2 GeV in order 
to evaluate the stability of the fit, and a systematic uncertainty was estimated from the maxi- 
mum deviation from the central value observed. The differences were found to be smaller than 
the statistical uncertainties from the fit in 15 out of the 16 E'Tc-rj-regions. Only for Z — )• e7fake 
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Table 10.1: Overview of the estimation of the multijet background in the Z — )■ ee and 
Z —7- e7fake samples. The estimated fraction of the multijet background is given. 
In the case of the Z — )• e7fake sample, the statistical uncertainty from the fit as well 
as the estimated systematic uncertainty are presented. For the Z — )• ee sample, 
the systematic uncertainty was estimated conservatively to be as large as 100%. 



events with photons in the region 0.6 < jryj < 1.37 and 15 GeV < < 20 GeV, a systematic 
uncertainty of 15% was found given a statistical fit uncertainty of 9%. The systematic uncer- 
tainties from the limit variation were added in quadrature to the uncertainties retrieved from 
the comparison with the alternative modelling. 

In the Z — )• ee sample, the multijet background was estimated using the template fit method 
only. The fraction of multijet events is low, even in the high- and low-mass regions, so that 
the kinematic bias of the Z peak due to the minimal requirements on the of the electrons 
dominates and does not allow for a fit with the CB*BW signal model. This is evident from the 



right plot in Fig. 10.6, which shows the invariant mass distribution for the events from Z — )• ee 
sample with an E^ of the probe electron between 40 GeV and 45 GeV. Clearly, the lower side 
of the spectrum cannot be described by a CB*BW signal model. The template fit, however, 
describes the distribution reasonably well. The fraction of multijet events was estimated in 
four regions of E^, similar to the estimate in the Z — )• e^f^ke sample. An ry-dependence was 
not accounted for, however, but a possible dependence is believed to be covered by the large 



systematic uncertainty assigned. The left plot in Fig. 10.6 shows another example fit in the 
region [15 GeV, 20 GeV). 



The results of the multijet estimate in the Z — )• ee sample are presented in Tab. 10.1 for 
low-i?T probe electrons, a fraction of 2.3% was estimated. It decreases with increasing electron 
St- The fraction of multijet background was found to be small in the Z — t- ee sample. However, 
also the effect on the overall estimation of fe^-y was evaluated to be small (see below), and a 
conservative systematic uncertainty of 100% was assigned, which turned out to be negligible 
with respect to other systematic uncertainties affecting /e-5.7. 
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10.1 Estimate of the misidentification rate 




Figure 1 0. 6: Invariant mass distributions for the two electrons from the Z — )• ee selection in 
data for the probe electron in the region 15 GeV < Et < 20 GeV (left), and in 
the region 40 GeV < E^ < 45 GeV (right). The distributions were fitted with a 
template derived from MC simulations (solid line), and an exponential for the 
background modelling (dashed line). 
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Figure 1 0. 7; Ratio of the number of events in the Z — )• e7fake and in the Z — )• ee sample in 
four bins of r] and four bins of E't in data. The yields were corrected using the 
estimate for the multijet contribution. 
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Figure 10.8: Scale factors for the misidentification rate fe^-y in 16 bins in E't and |r/|. The 
error bars show the sum of statistical and systematic uncertainties as presented 



in Tab. 10.2 
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10 Background events with electrons misidentiEed as photons 
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10.1 Estimate of the misidentification rate 



Derivation of the scale factors 



The scale factors were calculated according to Eq. (10.1 ): the number of events in the Z — )• e7fake 
and Z — )• ee samples in data were corrected using the multijet background estimate derived 
above. Fig. 10. 7| shows the ratio of the number of events in the Z — )• e7fake and the Z — )• ee 



samples Ng^/Nee in four bins of rj and four bins of Ei[ in data. 

For the MC simulations, this ratio was calculated using only events in which the tag electron 
and the probe electron or photon, respectively, were matched to the generated electrons from 



the Z decay within a cone of size 0.2 in r/-i;^-space. Fig. 10.8 shows the SFs in the 16 bins in 
and \r]\ together with the combined statistical and systematic uncertainty as evaluated in the 
next section. 

The SFs vary with the Et and the rj of the probe object, and it was hence necessary to derive 
them as a function of the probe kinematics. In general, the SFs are larger for lower E't and close 
to unity for transverse energies between 30 GeV and 45 GeV. The SFs in the lowest E'x-bin 
for < Ir^l < 0.6 and 0.6 < \7]\ < 1.37 are much larger than 1, but they are also subject to large 
uncertainties. 

The kinematic range of the probes was limited to 45 GeV in E^. However, the SFs in the 
ST-bins [30 GeV, 40 GeV) and [40 GeV, 45 GeV) are compatible within the uncertainties in all 
r/-regions, and, hence, it is well-motivated to use the SFs for the bin [40 GeV, 45 GeV) also for 
higher transverse energies. 



Evaluation of systematic uncertainties 

Different sources of systematic uncertainties were considered for the /e-s.7 SFs, an overview of 



which is given in Tab. 10.2 together with the statistical uncertainty and the quadratic sum of 
all uncertainties. 

The uncertainty on the multijet background in the Z — )■ ee and Z — t- e7fake were considered as 
well as several electron- and photon-specific uncertainties, the evaluation of which is described 



in Ch. 12 uncertainties on the electron SFs for triggering, reconstruction and identification, 
uncertainties on the electron energy scale and resolution, and uncertainties on the photon iden- 
tification, energy scale and resolution were considered. 



The SFs for electrons detailed in Sec. 5.1 were only derived for electron transverse energies 
larger than 25 GeV. However, SFs for the range 15 GeV to 25 GeV were derived for electrons 
without isolation requirement. Since the efficiency of the isolation cut is as large as roughly 97% 
for electrons, an additional uncertainty of 3% was added to the SFs below 25 GeV in order to 
conservatively account for the missing isolation cut. 

The dependence on the pile-up conditions was evaluated by comparing Nee and Ne'y in a 
low pile- up region with 1 ~ 5 primary vertices in the event to a high pile- up region with 6-10 
primary vertices. The differences with respect to the nominal SFs are of the same order of mag- 
nitude as the statistical uncertainties, so the observed effect could well originate from statistical 
fluctuations. The largest difference between the low pile-up and high pile-up regime with respect 
to the nominal SFs was taken as a symmetric systematic uncertainty, which is believed to be 
very conservative, because the MC simulations were reweighted to the pile-up conditions present 
during data taking. 

Apart from the uncertainty on the pile-up conditions and the statistical uncertainty, the 
following systematic uncertainties dominate the total SF uncertainty: the estimate of the multijet 
background in the Z — )• e7fake sample, and the SFs for electron and photon identification, which 
feature large uncertainties in particular for electrons and photons with low E't- 
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10.2 Application to processes with two leptons in the final state 



The SFs for fe^-y were applied to simulations of background events with a final state electron 
which was misidentified as a photon and an additional real lepton. Processes of interest are 
dileptonic decays of top quark pairs, Z — )• ee and Z — >• rr events with additional jets, single top 
production in association with a leptonically decaying W boson, and diboson production with 
at least two final state leptons. 
Tab. 



10.3 



shows the resulting expectations for 1.04 fb . The uncertainties on the expectations 



include all systematic uncertainties as discussed in Ch. 12 , in particular the uncertainty on /e^^ 
and uncertainties due to detector modelling. Uncertainties on the MC normalisation from the 
uncertainties on the cross sections and the luminosity as well as uncertainties due to limited 
MC statistics. Uncertainties which affect /e-s>7 as well as the selection efficiencies in the MC 
simulations were treated correlated. 

The negative estimate in the muon channel from single top Wt production is due to negative 
weights in the MC@NLO generation and the value of —0.1 is compatible with zero within the 
uncertainties. The largest contribution is from dileptonic tt decays, which features a larger 
yield in the muon channel due to the looser event selection with respect to the electron channel 



(Sec. 6.1). 



Process 


e+jets 


^+jets 


Dileptonic ti 
Z+jets 

Single top VFt-channel 
Diboson 


6.8 ± 2.3 

17 + 

- 1.7 

n 99 + 0.25 
'-'•^^ - 0.22 

04 + 

^■^^ - 0.04 


9.6 ± 2.7 

7 + 

- 0.7 

-0.10 ± 0.10 

00 + 

^•'-'^ - 0.00 



Table 10.3: Expected yields for processes with an electron misidentified as a photon for 
1.04 fb~^ from MC simulations using the SFs for the electron-to-photon misiden- 
tification rate. The uncertainties quoted are the total systematic uncertainties. 
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11 Background events with prompt photons 
in the final state 



The treatment of hadrons and electrons misidentified as photons has been discussed in the 
previous chapters. The remaining backgrounds are processes with prompt photons in the final 
state. 



Sec. 11.1 discusses the small background contribution from ttj production which does not 
fulfil the signal phase space cuts. Background estimates of multijet and VF+jets production 
with an additional prompt photon in the final state were derived with data-driven approaches 
as described in Sec. 11.2 and |11.3 respectively. The remaining background contributions from 
Z+jets, single top and diboson production in association with a prompt photon were estimated 
using MC simulations (Sec. 11.4). 



11.1 tt^ production outside of the signal phase space 

The tt'j signal is only well-defined within a certain phase space, as for example defined by the 



invariant mass cuts required in the signal modelling in this analysis (Sec. 4.1). In rare cases 



however, tt^ events outside of this signal phase space may fulfil the event selection. 



This is illustrated in Fig. 11.1, which shows a schematic representation of the angular phase 
space in r] and (j) for photons in tt events. Regions around charged massless fermions, such as 
light quarks and charged leptons, are forbidden for true photons, because the invariant mass 
cuts applied in the WHIZARD generation effectively translate into minimal angular distances: 



2 



m 



{Pf + P-ff = 2- Ef Ey ■ {1 - cosa) > ml^^ 



a > arccos | 1 - 

I- Ef ■ E^ 

with a the angle between the charged massless fermion and the photon. The minimal value for 
a depends on the fermion and photon energies. 

True photons and the corresponding reconstructed photons are indicated by open squares and 



crosses in Fig. 11.1, respectively. Two cases are depicted: case (1) shows a photon from the tfy 
signal phase space, while case (2) shows a photon from a background tt'j event: the true photon 
lies outside of the signal phase space, because it is too close to a light quark. However, the 
position in rj-(j)-space of the reconstructed photon differs from the position of the true photon 
and lies actually within the measured signal phase space. 

In order to estimate the contribution from background ttj events, tt events simulated with 
MC@NLO were used: the tt sample also contains photon radiation in the phase space not 



generated with WHIZARD through QED corrections provided by the PHOTOS package 120 
However, interference effects are not correctly taken into account and it is hence indicated to 
assign a conservative systematic uncertainty to the tt'j estimates obtained from this sample. 

The MC@NLO sample does not only contain tt'j events outside of the signal phase space, but 
also events within the latter. In order to avoid double-counting, the overlap with the WHIZARD 
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signal phase space 
(WHIZARD) 



true Y / 'I \ \ 

l.y) \ 

''^x reco. y 




V^/ X reco. Y 
truey/g^ 






charged H 


^ light quark 


lepton 



Figure 11.1: Illustration of photons within and outside of the tt'y signal phase space: the 
plot shows the angular phase space in photon r] and (j). Regions around charged 
massless particles, such as light quarks and charged leptons, are forbidden for 
true photons. True and the corresponding reconstructed photons are indicated 
by open squares and crosses, respectively. Case (1) shows a photon from the tt7 
signal phase space, while case (2) shows a photon from a background event. 



tt^ sample was removed from the MC@NLO tt simulation, as discussed in Sec. 4.2 After the 
overlap removal, an estimate of 0.8 and 1.3 background ti^ events was derived for 1.04 fb^^ of 
data in the single electron and single muon channels, respectiveljj^ 

The origin of these events was found to be due to collinear radiation from leptons and quarks. 
In the following, it is first shown that the event selection applied in this analysis is well-suited to 
select only events from the phase space generated in the ^7 simulation, and that the invariant 
mass cuts applied in the event generation were hence not chosen too loosely. Then, it is shown 
that background it7 events are indeed due to collinear radiation from leptons and quarks. 

Check of the invariant mass criteria 

The contribution of background ti'^ events was found to be as large as roughly 4% of the yields 



predicted by the WHIZARD tt^ simulation (Sec. 6.3). This could indicate that the signal phase 



space generated with WHIZARD is slightly too small for the considered event selection. 



Fig. 11.2 shows the smallest distance in 7/-i?i>-space between light quarks and true photons for 
diff'erent invariant mass cuts m(g, 7) of 5 GeV, 10 GeV, and 25 GeV. For all three invariant 
mass cuts, a turn-on at a certain typical Ai? is observed. The turn-on has a finite width because 
of the energy spectra of quarks and photons. As expected, the turn-on value increases for larger 
invariant mass cuts. 

In order to study the effect of the generation level cuts on the angular distances of recon- 
structed objects, the angular resolution needs to be considered. Jets are expected to feature the 



^Only reconstructed photons were considered which were matched to a true photon with a minimal pT of 10 GeV. 
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11.1 tt"f production outside of the signal phase space 




Figure 11.2: Smallest distance in r/-</)-space between light quarks and true photons for different 
invariant mass cuts m{q, 7) in WHIZARD ti^ simulations: 5 GeV (thin solid line), 
10 GeV (dashed line), and 25 GeV (thick sohd line). 




Figure 11.3: The left plot illustrates the angular resolution of jets: the smallest distance in 
7?-0-space between a light quark and the closest jet is shown. The right plot shows 
the smallest distance in ry-^-space between jets matched to a light quark and true 
photons for different invariant mass cuts m(g, 7) in WHIZARD ti^ simulations: 
5 GeV (thin solid hue), 10 GeV (dashed line), and 25 GeV (thick sohd line). The 
dotted vertical line shows the minimal Ai? applied between jets and reconstructed 
photons in the event selection. 



worst angular resolution of the considered objects: the left plot in Fig. 11.3 shows the distance 



between a light quark and the closest jet. The angular resolution is of the order of only 0.15 in 



The right plot in Fig. 11.3 shows the smallest distance in r/-i;^>-space between jets and true 
photons, again for invariant mass cuts of 5 GeV, 10 GeV and 25 GeV. The jets were required 
to be matched to a light quark within Ai? < 0.3. For large values of Ai?, the same trend as in 



Fig. 11.2 is observed with slightly broadened turn-on curves due to the angular resolution. At 
values of Ai? < 0.4, a peak is observed which is more pronounced for the lower invariant mass 
cuts than for the cut at 25 GeV. This peak is due to true photons which are close to a quark, 
which causes the photon and the jet to merge and form one single jet. At values of Ai? > 0.4, 
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photons and jets are well separated, which is consistent with the distance parameter of 0.4 used 
for the anti-kt jet algorithm (Sec. 5.3). 

Generally, the direction of reconstructed photons does not differ significantly from the direction 
of the corresponding true photon. Hence, in order to avoid biases due to jet-photon merging, 
reconstructed photons were required to be separated by more than 0.5 in rj-cp-space (Sec. 5.6), 
as indicated by a dotted vertical line in the right plot of Fig. 11.3 As seen in the plot, the cut 
value of 0.5 is well above the turn-on for an invariant mass cut of 5 GeV. 

Hence, the angular cuts applied in the event selection appear to be sufficient to select only 
events from the signal phase space. This does not only hold for the distances between jets 
and photons, but also for electrons and muons, which are separated from photons by minimal 
distances in r]-(p-space (Sec. 5.6). It is therefore concluded that the origin of the background tt^ 
events is not due to inappropriately high invariant mass cuts in the signal definition. 



Photon radiation from charged leptons 



Fig. 11.4 shows the distance in ?7-0-space between true photons and the closest true lepton in ttj 
background events in the electron (left plot) and in the muon channel (right plot), respectively. 
Photons at low AR values were radiated from the corresponding lepton, as checked with MC 
truth information. 

Experimentally, photons and leptons have a minimal distance in r]-(j)-space (Sec. 5.6). Hence, 



events with collinearly radiated photons can only pass the ttj event selection if the corresponding 
lepton was not identified as a lepton object according to Sec. 5.1 or 5.2 These events are 



mostly dileptonic tt decays with one well-identified lepton and one lepton with a coUinear photon 
radiation. They are a natural part of the tt^ background, because such events were not simulated 
with WHIZARD due to collinear divergences in the ME calculation. 



0.4 



0.2 



ATLAS work in progress 
simulation 



0.4- 



0.2- 



ATLAS work in progress 
simulation 



0.2 



0.4 



0.6 



0.8 1 
A R(y, e) 



0.2 



0.4 



0.6 



0.8 1 

A R(Y, ii) 



Figure 11. 4: Distance in rj-(/)-space between true photons and the closest lepton in ttj back- 
ground events from MC simulations with MC@NLO in the electron channel (left) 
and in the muon channel (right). In both plots, the last bin includes the overflow 
bin. 



A summary of the contributions from radiation from electrons, muons and r- leptons is pre- 
sented in Tab. 11.1 most photons are radiated from electrons. The radiation from muons has 
a smaller contribution. Radiation from r-leptons contributes less than 0.1 events in both chan- 
nels. Altogether, collinear radiation from charged leptons represents 0.7 and 1.1 background 
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11.2 Multijet production with a prompt photon 



tt'j events in the electron and muon channel, respectively. The systematic uncertainty on these 
contributions was conservatively estimated to be as large as 100%, i.e. 0.7 and 1.1 events. 

Photons from jet fragmentation 

Background tt'j events in single lepton tt decays are mostly due to photon production in jet 
fragmentation processes, called bremsstrahlung in the following. Bremsstrahlung is partly cov- 
ered by the WHIZARD simulation in cases where the final state photon fulfils the signal phase 
space criteria. Contributions from outside of the signal phase space were estimated using photon 
radiation in the MC@NLO tt simulation. For 1.04 fb~^, a contribution of about 0.3 events in 
both lepton channels together was estimated (Tab. [TlT| ). 



In Sec. 4.2 ambiguities in the removal of the signal phase space from the MC@NLO simulation 
due to bremsstrahlung processes were discussed: there are several possibilities for the choice of 
the light quark to be used in the requirement on the invariant mass of light quarks and final state 
photons. In order to assign a systematic uncertainty arising from these ambiguities, different 



possibilities for the choice of the light quark were considered (see Fig. 4.1 for an illustration of 
the particles with the difi'erent status codes): 

• Variation (a) : light quark with status code 123/124 

• Variation (b) : light quark with status code 143/144 

• Variation (c) : light quark with status code 2 



The yields for variations (a) - (c) are shown in Tab. 11.2 together with the default case, for 
which the four-momenta of all particles in the parton shower were added up except for the photon 
and the combined four-momentum was used as an estimate for the quark four-momentum after 
photon radiation. The expected yields for 1.04 fb~^ vary, and the largest yield was considered 
as a conservative estimate of the systematic uncertainty. 

Estimate and systematic uncertainty 

The total prediction of background tt'j events reads 0.8 events in the electron and 1.3 events 
in the muon channel. The systematic uncertainties on the difi'erent contributions are added 
linearly to conservatively account for correlations between the different sources of background 
tt'j events. The uncertainties add up to 1.1 events in the electron and 1.9 events in the muon 



channel, respectively, leading to background ttj estimates of 0.8^Q g and 1.3^^'^ 



11.2 Multijet production with a prompt photon 

Multijet events do not feature prompt electrons or muons, but may involve the production of 
prompt photons such as in 7+jet production, or may feature real photons originating from jet 
fragmentation. The background from multijet production with prompt photons was hence esti- 
mated in two steps: in a first step, the matrix method was used to estimate the contribution from 



processes which lead to good electron and good muon objects after the preselection (Sec. 6.1), 
so-called fake leptons. Examples for such processes are leptonic decays of heavy flavour mesons, 
and hadrons from jet fragmentation, which may be misidentified as electrons if a large fraction 
of the energy is deposited in the electromagnetic calorimeter. 

The matrix method, discussed in detail in the following, provided a sample of events which was 
used to estimate the yield as well as differential distributions of the multijet background. Hence, 
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e+jets 


/i+jets 


Radiation from electrons 
Radiation from muons 
Radiation from r-leptons 
Bremsstrahlung 


0.50 ± 0.50 
0.16 ± 0.16 
0.04 ± 0.04 

0.08 t rjs 


0.74 ± 0.74 
0.23 ± 0.23 
0.09 ± 0.09 
20 + 

^■'^^ - 0.20 


Sum 


0.8 t li 


1 o + 1.9 
J^.-J - 1.3 



Table 11.1: Overview of the different contributions from background ti^ events estimated in 
ti events simulated with MC@NLO for 1.04 fb""*^ together with the associated 
systematic uncertainties. 



Bremsstrahlung definition 


e+jets 


^+jets 


Variation (a) 


0.37 


0.69 


Variation (b) 


0.05 


0.16 


Variation (c) 


0.15 


0.38 


Default 


0.08 


20 +°-^^ 

^■'^^ -0.20 



Table 11.2: Overview of the systematic variations for the background tt^ contribution from 
bremsstrahlung for 1.04 fb~^ estimated in events simulated with MC@NLO. For 
the default estimate, the resulting systematic uncertainty is shown. 



in a second step, the final event selection (Sec. 6.2) was applied to this sample, which yielded 
a sample of events with a fake lepton and a photon candidate. The fraction of real photons 
within these candidate photons was estimated using a template fit to the p^^^'^^ distribution as 
described in Ch. [7j The templates for prompt photons and hadrons misidentified as photons as 
derived in Ch. [8] and [9] were used for this fit. 



Fake lepton contribution estimated with the matrix method 

The matrix method is a data-driven technique to estimate the amount of fake leptons with a 
minimised dependence on MC simulations. It was successfully used for top quark analyses in 
the single lepton decay channel at ATLAS, cf. for example Ref. [3| [44|[55|[58l i [62|[67] . The method 
is based on an additional loose lepton definition compared to the tight lepton definition used 
for the actual analysis. By replacing the tight by the loose definition, a loose event selection is 
obtained, which defines a loose data sample. 



The tight definitions for electrons and muons were those presented in Sec. 5.1 and 5.2 For 
the loose electron definition, the tight electron menu was replaced by the medium menu with 
looser shower shape requirements. Additionally, a hit in the Pixel fe-layer was required to reject 
backgrounds from converted photons. The electron isolation cut was loosened to 6 GeV instead 
of 3.5 GeV. For the loose muon definition, the requirements on the track and calorimeter 
isolations were disregarded. 

The number of events in the loose sample Nl consists of Nj^'^^ events with real leptons and 
jyfake gyents with fake leptons: 

Nl = Nl^""^ + Ni^^\ (11.1) 
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The same holds for the number of events in the tight sample Nt'- 

Nt = A^t''^ + N^t"^"" = e ■ iV£^^i + / • iVf^'^" , (11.2) 

where e and / are the probabilities for real and fake leptons in the loose sample to also fulfil 
the tight lepton definition. If e and / are known in addition to and Nt^ the number of 



events with fake leptons in the tight sample N^^^ can be calculated by solving Eq. (11.1) and 



Eq. (11.2). The solution reads: 



/ ' e-f 



with 



non— tight jirnon— tight . at 

w^ • ^ + wt • jSt , 



non— tight j p ^ 

, and wj^ = —J • 



£-/' ' e-f 

Not only the total yield of multijet events was estimated, but also kinematic distributions of 
the multijet background were extracted by weighting events in the loose sample: the weight w^^ 
was applied to events which fulfilled the tight event selection; the weight w^°'^~*^^'^* was applied 
to events which fulfilled only the loose criteria, but not the tight criteria. Typically the following 
relations hold: 

, n non— tight ^ n I I non— tight 

Wt < U , Wj^ " >U, \Wt\ <^ W^ 

The efficiency for real leptons e was measured in Z — t- events using the tag-and-probe 
method for electrons and muons. The fake efficiency / was estimated in a CR dominated by 
multijet production: 

5 GeV < < 20 GeV (e + jets) , 
< 20 GeV and + "1t < 60 GeV (M+jets) . 

Contributions from processes with real leptons, as for example from VF+jets and Z+jets pro- 
duction, were subtracted from the yields in the CR. 

In the electron channel, real and fake efficiencies were parametrised as a function of the electron 
rj in order to account for the varying detector geometry. In the muon channel, the efficiencies 
were parametrised as a function of the muon rj, and as a function of the highest jet present 
in the event: a high px indicates a lot of hadronic activity, which - in turn - is correlated to the 
muon isolation. 



Fig. 11.5 shows distributions in the respective CRs for the electron and the muon channel after 
the 5-tagging requirement. The CRs are largely enhanced in multijet production and the data 
distributions as well as the expectations from MC simulations for various processes are shown 
together with the multijet estimate from the matrix method. The two plots in the upper row 
show the electron distribution (left) and the number of jets (right) in the electron channel. 
The two plots in the lower row show the pT distributions of the muon (left) and of the jets 
(right) in the muon channel. 

The agreement of the sum of the expectations with data is reasonable. Small disagreements 
can be explained by the inclusion of the photon in the definition for this analysis, which was 
not included in the definition of the CRs used in the derivation of the lepton misidentification 
rates. In principle, this small change in the CR is expected to have a small effect on the overall 
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multijet estimate, which is confirmed by the reasonable agreement observed in Fig. 11.5 These 
small differences between data and the sum of the expectations are covered by the systematic 
uncertainty of 100% assigned to the multijet background estimate in the signal region. 

The yields obtained after the preselection with and without the 6-tagging requirement were 
shown in Tab. 6.1 The uncertainty on the yield was estimated to 50% for events without the 



6-tagging requirement and 100% in events with at least one 5-tagged jet. 



Prompt photon fraction estimated with the template fit 

In order to identify multijet events with an additional prompt photon, only events with a photon 



object according to the definition in Sec. 5.6 were selected from the weighted loose sample. Since 
loose electron objects may also be identified as photons, a minimal distance in jy-c/i-space was 
required between loose electrons and photons. 

This yielded a sample correctly normalised for the presence of a fake lepton and a photon 
candidate. The photon candidate could either be a prompt photon or a hadron misidentified 
as a photon, while the aim was to estimate only the contribution featuring prompt photons. 
As discussed in Ch. [7| p^^^"^^ is a good discriminating variable between prompt photons and 
hadrons misidentified as photons. The p™'^®^'^ distributions for the photon candidates in the 



weighted loose sample are shown in Fig. 11.6 for the electron channel (left) and for the muon 
channel (right), respectively. Both distributions feature a low number of expected events and 
show contributions in the low-p^'^^^'^-region, which is dominated by prompt photons, as well as 
in the high-pJj?'^^^''-region, which is largely dominated by hadrons misidentified as photons. The 
total number of events reads 2.7 events in the electron and 2.2 events in the muon channel. 

In order to isolate the contribution from prompt photons, a template fit to the p^^^'^^ dis- 
tribution, as introduced in Ch. [7| was performed in both lepton channels separately. However, 
a template fit cannot be performed to data distributions which feature fractional numbers of 
events for which the Poissonian probability density function (pdf) is not well-defined. This issue 
was solved by noting that the estimate presented in Fig. 11.6| was dominated by the events with 



positive weights, while events with negative weights constitute only a small correction. Hence, 
the following procedure was applied: the purity of prompt photons within the photon candidates 
was estimated using only the events with a positive weight, and the result was then rescaled 
to the weighted yields (2.7 and 2.2 events, respectively). This procedure is valid, because the 
photon isolation and the loose/tight classification of the electron object are uncorrelated. 

For the template fit, the prompt photon template was taken as derived in Ch. [Sj because it 
depends only marginally on the photon kinematics. The hadron fake template, however, needed 
to be rederived with respect to Ch. |9]by estimating the E't spectrum of hadron fakes as well as 
the fraction of hadron fakes with \7]\ > 1.81. For the muon channel, the same procedure for the 
reweighting of the hadron fake template was applied as presented in Ch. |9j 

In the electron channel, the statistics in the CR was too low to estimate the spectrum 
and the fraction of photons with \r]\ > 1.81 for the hadron fakes. The template estimated for 
the ti^ topology was used instead (Ch. |9]) and for the evaluation of systematic uncertainties, the 
extreme templates in E't and rj were used, as presented in Fig. |9.2[ 



Fig. 11.7 shows the resulting template fits in the electron channel (left) and the muon channel 



(right), respectively. In the electron channel, the fit yields 5.7 out of 13 events with a real 
photon. In the muon channel, the result reads 1.6^'j^ g real photon events out of 7 events. The 
mode of the marginalised pdf was used as estimator and the statistical uncertainty was estimated 
using the smallest interval containing 68% of the marginalised pdf. 
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Figure 11.5: Distributions in control regions largely enhanced in multijet background for data 
and the expectations from MC simulations as well as the multijet estimate from 
the matrix method: the two plots in the upper row show the electron E'y dis- 
tribution (left) and the number of jets (right) in the electron channel. The two 
plots in the lower row show the pT distributions of the muon (left) and of the jets 
(right) in the muon channel. In all plots, the last bin includes the overflow bin. 
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Figure 11.6: p^'^^'^^ distribution for photon candidates in events with a fake lepton in the 
electron channel (left) and the muon channel (right), respectively. In both plots, 
the last bin includes the overflow bin. 
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Figure 11.7: Template fit to the photon p^'^^'^^ distribution in events with a fake lepton in the 
electron channel (left) and the muon channel (right), respectively. Only events 
with a positive weight from the matrix method are shown and their respective 
weights were not used in the distribution. The results of the fit yielded the 
contributions from real photons and hadrons misidentified as photons. In both 
plots, the last bin includes the overflow bin. 
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Figure 11.8: Calibration curves in the electron (left) and the muon channel (right): pseudo- 
data were created with a variable amount of real photons (x-axis). The amount 
of hadrons misidentified as photons was set to the value from the fit to data 



(Fig. 11.7). The mean output of real photon events from 5000 sets of pseudo- 
data is shown in the y-axis. A parabola was used to parametrise the resulting 
curves. 



In order to check if the fit result correctly estimated the real underlying number of multijet 
events with prompt photons, the following test was performed: pseudo-data were created by 
fluctuating each bin of the p^'^'^'^^ distribution individually around the expectation value ac- 
cording to the Poissonian pdf. While the expectation for the hadron fake contribution was fixed 
to the global mode of the actual fit, the expectation for the prompt photon contribution was 
varied between and 10 events. 



For each input value for the number of prompt photons events. Fig. 11.8 shows the correspond- 
ing mean output of the fit from 5000 pseudo-data distributions. For larger input values, the fit 
output corresponds to the input value. For lower input values, however, the fit was found to 
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be biased, because the Poissonian pdf is asymmetric around its mode for small expected values, 
and the mean of the pdf is larger than the expected value. A parabola was used to parametrise 
these calibration curves. 



Source 




Uncertainty 






e+jets 


^+jets 


Yield 


1.2 


0.3 


Statistical 


+ 


0.7 
0.9 


+ 0.9 
0.3 


Fake lepton normalisation 




1.2 




0.3 


Electron to photon extrapolation 




0.04 




0.02 


Fraction of converted prompt photons 




0.04 




0.02 


Pile-up dependence of the signal template 




0.03 




0.02 


Reweighting of the background templates (-Et) 




0.35 




0.05 


Reweighting of the background templates (r/) 




0.68 




0.12 


Fraction of converted hadron fakes 




0.07 




0.05 


Pile-up dependence of the background template 




0.03 




0.03 


Sum (stat. + syst.) 


+ 


1.6 
1.2 


+ 1.0 
0.3 



Table 11.3: Different sources of uncertainties on the estimate of multijet events with an addi- 
tional prompt photon: the statistical uncertainty from the template fit, the uncer- 
tainty on the amount of fake leptons, and the systematic uncertainties associated 
to the template fit are presented. 



The results retrieved from the fit were corrected using the parabolas from Fig. 11.8 and 



rescaled to the actual weighted multijet estimate. The final estimates read 1.2 ^[J g and O.Sl^Qg 
in the electron and muon channel, respectively. The uncertainty is the statistical uncertainty 
only as retrieved from the template fit. 

Systematic uncertainties were evaluated for the fake lepton estimate as well as for the amount 
of prompt photons within the photon candidates. An overview of all uncertainties considered is 
shown in Tab. 111.31 

The systematic uncertainty on the amount of fake leptons was estimated to 100% of the 
final yield. The systematic uncertainties on the amount of prompt photons were estimated 
by replacing the nominal prompt photon and hadron fake templates in the template fit by 
systematically varied templates. Pseudo-data were constructed and 5000 p^'^^'^^ distributions 
were fitted to evaluate the impact of each systematic effect. The procedure is described in detail 
in Ch. 12 where also the different templates from the systematic variations are presented. 

The following sources of systematic uncertainties were considered for the prompt photon tem- 
plate: the uncertainty on the extrapolation from electron to photon templates; the uncertainty 
on the fraction of converted photons used for the derivation of the template; the uncertainty on 
a possible dependence of the template on the pile-up conditions. 

For the hadron fake template, the following sources of systematic uncertainties were consid- 
ered: uncertainties on the reweighting of the template in E't and r]] the uncertainty on the 
fraction of converted photons used for the derivation of the template; the uncertainty on a pos- 
sible dependence of the template on the pile-up conditions. The variations for the reweighting 
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of the hadron fake templates were increased with respect to the description in Ch. 12 as stated 
above. 

Dominant uncertainties were found to be the statistical uncertainty from the template fit 
as well as the uncertainty on the amount of fake leptons. In the electron channel, also the 
uncertainty due to the reweighting of the fake hadron template was found to be sizable, which 
is due to the low statistics in the CR used to estimate the reweighting parameters. The final 
estimate in the electron channel reads: 1.2 ^}'2 (stat. + syst.). The estimate in the muon channel 
reads: 0.3 (stat. + syst.). 



11.3 VF+jets production with a prompt photon 



The background contribution from PF+jets production with an additional photon in the final 
state was estimated using a data-driven approach in order to reduce the dependence on MC 
simulations. The event selection from Ch. [6] was modified in order to define a CR with an 
enhanced VK+jets+7 contribution: one to three jets were required instead of four jets and none 
of the jets needed to be 6-tagged. The yield in the signal region was then extrapolated from the 
CR using MC simulated Vl^+jets+7 events. 

The event selection in the CR yielded 1446 events in the electron and 2468 events in the muon 
channel. However, background contributions were due to processes with prompt photons in the 
final state, but also from events with electrons and hadrons misidentified as photons. 



The upper part of Tab. 11.4 shows the expected background contributions for events with a 
prompt photon in the final state: ti, Z+jets, single top, and diboson production with an addi- 
tional photon in the final state were estimated using MC simulations. Systematic uncertainties 
include uncertainties on the expected cross sections for the individual processes, as well as un- 



certainties due to detector modelling effects as discussed in Sec. 12.3 and uncertainties due to 



limited MC statistics (cf. Sec. 12.2). For the contributions from Z+jets+7, single top+7 and 



diboson+7 production, an additional uncertainty of 100% was added as discussed in Sec. 11.4 



Background contribution 


e+jets 


/x+jets 


ttj 


36 


± 


8 


49 


± 


11 


Z+jets+7 


80 


+ 


100 

80 


260 


+ 


290 
260 


single top+7 


6.6 


± 


6.6 


8.1 


+ 


8.2 

8.1 


diboson+7 


6.8 


+ 


7.2 

6.8 


12.4 


+ 


12.6 
12.4 


multijet+7 


80 


± 


40 


75 


± 


38 


ttj (e — )• 7) 


1.4 


± 


0.4 


1.2 


± 


0.3 


tt (e 7) 


48 


± 


9 


70 


± 


13 


Z+jets (e — )• 7) 


220 


± 


130 


14 


± 


9 


single top (e — 7) 


5.8 


± 


1.3 


6.8 


± 


1.3 


diboson (e — )• 7) 


9.1 


± 


2.4 


8.4 


± 


2.2 



Table 11. 4: Expected yields for different background contributions in the VF+jets CR in the 
electron (left) and in the muon channel (right). The upper part shows the expec- 
tations from events with prompt photons in the final state. The lower part shows 
the expectations from events with electrons misidentified as photons. 
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The contribution from multijet events with prompt photons was estimated using the approach 
with the matrix method and a template fit to the p^^^'^^ distribution. 



11.2 



discussed in Sec. 

Estimates of 80 it 40 and 75 it 38 events were obtained in the electron and muon channel, re- 



spectively. In contrast to the estimate derived in Sec. [11.2 for the signal region, the reweighting 
of the hadron fake templates in E't and r] was successfully performed in both lepton channels, 
because of increased statistics in the CR enhanced in hadron fakes. Moreover, the result of 
the template fit turned out to be unbiased given the larger expected yields, and therefore a 
correction as depicted in Fig. 11.8 was found to be not necessary. 

The lower part of Tab. 11.4 shows the expected contributions from background events with an 
electron which was misidentified as a photon. The yields were estimated using the scale factors 



for the electron-to-photon misidentification rate /e->7 as described in Ch. 10 The quoted un- 



certainties include the systematic uncertainties on the expected cross sections for the individual 
processes, the uncertainties due to detector modelling effects, uncertainties due to limited MC 
statistics and the uncertainties on the fe^-y scale factors. 

The largest background contributions in the CR with either prompt photons or electrons 
misidentified as photons were found to be from Z+jets and multijet production. These con- 
tributions also feature the largest absolute uncertainties on the background predictions. As 
expected, the Z+jets production with a misidentified electron contributes significantly less in 
the muon compared to the electron channel: the only contribution in the muon channel is from 
Z — )• rr events with one r-lepton decaying to a muon and one r-lepton decaying to an electron, 
which is then misidentified as a photon, while in the electron channel Z — )• ee events contribute. 
Additional sizable background contributions are due to tt'j production and tt production with 
a misidentified electron. 

The amount of events in the CR with hadrons misidentified as photons was estimated using 
a similar approach as used for the multijet+7 estimate: a template fit to the photon p^^^"^^ 
distribution was used to distinguish prompt photons and misidentified electrons from hadron 
fakes as described in Ch. [7} While the signal template was used as derived in Ch. [Sj the template 
for misidentified hadrons needed to be rederived following the approach in Ch. [9| due to the Et 
and Tj dependence of the p^^^"^^ distribution from misidentified hadrons, the templates needed 
to be reweighted in i?T and r] using the CR dominated by hadron fakes as introduced in Ch. [9| 
The resulting hadron fake template represents the p^^^"^^ distribution for hadron fakes in the 
W"+jets+7 CR. 

Fig. 11.9 shows the resulting template fits in the VF+jets+7 CR in the electron (left) and 
the muon channel (right). The prompt photon template represents the VK+jets+7 process as 
well as all background processes with prompt photons and misidentified electrons as listed in 
Tab. 11.4 After subtraction of the background expectations, the resulting VF+jets+7 yields in 
the CR read 480 and 1190 itso in the electron and the muon channel, respectively, where the 
uncertainties represent just the statistical uncertainty of the template fit. 

These yields were extrapolated using VF+jets+7 simulations to the signal region with four 
jets out of which at least one had to be 6-tagged. The yields in the signal region read 1.84;toi3 
and 3.72 in the electron and the muon channel, respectively, where, again, the uncertainties 
are only the statistical uncertainties from the template fit. 

Tab. |11.5| shows the resulting Ty+jets+7 yields in the signal region using different control 
regions (CRs) in both lepton channels. While the default CR features one to three jets, CRs 
with only one, two or three jets were investigated. The uncertainties represent the statistical 
uncertainties from the fit only and the results in the different CR regions were found to be 
consistent within the fit uncertainties. Only the estimate in the muon channel using the three- 
jet-CR differs from the estimates from the other CRs, but the discrepancy is still only of the 
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Figure 11.9: Template fit to the photon p'^'^^'^^ distribution in events in the TF+jets+7 control 
region in the electron (left) and the muon channel (right), respectively, using the 
templates for true photons and for hadrons misidentified as photons. In both 
plots, the last bin includes the overflow bin. 



Control region 


e+jets 


/U+jets 


Default (1-3 jets) 

1 jet 

2 jets 

3 jets 


1.84 1 

1 70 -1- 0.18 
^•'"J - 0.18 

^•'-'O - 0.32 

1 70 -1- 0.59 
J^-'^ - 0.57 


q 79 -1- 0.18 

"^•'^ - 0.15 

Q Qv -l- 0.22 
3-»' - 0.18 

or,o + 0.32 
•5-^° - 0.36 

9 c:i + 0.60 
^•^^ - 0.70 



Table 11.5: Resulting VK+jets+7 yields in the signal region using different control regions 
(CRs) in the electron (left) and the muon channel (right): while the default CR 
features one to three jets, CRs with only one, two or three jets were investigated. 
The uncertainties represent the statistical uncertainties from the fit only. 



order of two standard deviation^ The result is hence considered to be stable with respect to 
the choice of the CR. 



Tab. 11.6 shows the contributions from different sources of uncertainties on the Vl^+jets+7 
yields in the signal region in both lepton channels. Uncertainties are due to the shape of the 
prompt photon and the hadron fake template, as well as due to various detector modelling ef- 



fects as described in Ch. 12 Additional sources of systematic uncertainties are the background 
estimates: for estimates derived from MC simulations, uncertainties originate from the uncer- 
tainties on the production cross sections, from the uncertainty on the luminosity, and from the 
uncertainties on the /e-s.7 scale factors. Moreover, the amount of background events with prompt 



photons in the final state was found to be subject to large uncertainties as described in Sec. 11.2 



for multijet production and in Sec. 11.4 for Z+jets+7, single top+7 and diboson+7 production. 

The largest sources of systematic uncertainties were found to be due to the shapes of the 
templates, the cross sections for the processes estimated from MC simulations, the background 
estimates from processes with prompt photons in the final state, and the modelling of jets, 
and 6-tagging. The uncertainty on the MC normalisation is much larger in the electron than in 



^The probability that at least one out of six measurements deviates by at least two standard deviations from 
the expectation value is roughly 24%. 
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Source of uncertainty 


e+jets 


/i+jets 


Statistical 


+ 


9.3 % 
7.1% 




+ 


4.9 % 
4.0% 


Template shapes 




12.6 ^ 


^0 




8.8 % 


Background MC cross sections and statistics 




15.0 ^ 


^^0 




5.6 % 


Prompt 7 background 




18.7 ^ 


^0 




22.0 % 


/e— >-7 




4.6 ^ 


^0 




0.8 % 


Luminosity 




3.2 ^ 


^^0 




1.3 % 


Electron modelling 




8.0 ^ 


^^0 




0.7 % 


Muon modelling 


± 


0.1 ^ 


^^0 




1.0 % 


Jet and modelling 




15.0 ^ 


1^0 




21.2 % 


Photon modelling 


± 


4.1 ^ 






1.4 % 


6-tagging modelling 




11.3 ^ 


^ 




17.8 % 


Modelling of LAr readout issues 




5.1 ^ 






4.2 % 


Total systematics 




36.3 ^ 




± 


37.2 % 


Total 




36 ^ 


^0 




37 % 



Table 11.6: Sources of uncertainties on the PF+jets+7 yield in the signal region in the electron 
(left) and in the muon channel (right). 



the muon channel. This can be understood since the electron channel features larger background 
contributions (Tab. 11.4) with respect to the H^+jets+7 expectation and is hence subject to 
larger uncertainties. The uncertainty on the jet and ^1 modelling is largely dominated by the 
uncertainty on the jet energy scale (JES). 

Most systematic uncertainties due to the detector modelling affect the expected yields in the 
CR as well as the yields in the signal region and are therefore expected to have a reduced impact 
on the extrapolation to the signal region using the VF+jets+7 simulation. However, this does 
not hold for the uncertainty due to the modelling of the 6-tagging performance, because a 6-tag 
was not required in the CR, which explains that the systematic uncertainty due to 6-tagging 
is of the same order of magnitude as the uncertainties due to the JES modelling, while for the 
final cross section (Ch. 13), the JES uncertainties dominate. 

The final estimates for VF+jets+7 in the signal region including statistical and systematic 
uncertainties read 1.8 it 0.7 and 3.7 it 1.4 in the electron and muon channel, respectively. From 
MC simulations 2.5 it 1.2 and 4.1 it 2.0 events were expected, where the uncertainty is just the 
48% uncertainty on the W^+jets normalisation (Sec. 4.2) and hence does not include additional 
systematic uncertainties. The data-driven estimate is hence consistent with the expectation from 
MC simulation within the uncertainties. Moreover, the uncertainty was reduced with respect to 
the uncertainty associated to the MC expectation. 



11.4 Additional backgrounds with a prompt photon 

The remaining background sources from Z+jets, single top and diboson production with an 
additional prompt photon in the final state were estimated using MC simulations as described 



in Sec. 4.2 In Tab. 11.7, the expected yields for 1.04 fb are presented. 
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Since the Z+jets, single top and diboson simulations only included photon radiation as QED 



corrections provided by the PHOTOS package 120 and not at ME level, the predicted yields 



are only approximations and hence subject to large systematic uncertainties. The uncertainties 
were estimated conservatively to 100% of the yields. Additionally, all systematic uncertainties 



as discussed in Ch. 12 were evaluated, in particular the uncertainties on the detector modelling, 
but also the uncertainties on the cross sections of the different processes, on the luminosity, and 
uncertainties due to limited MC statistics. 



Process 


e+jets 


^+jets 


Z+jets+7 
Single top+7 
Diboson+7 


1 o + 2.4 
'^■'^ - 1.3 

6 + o-'^ 

"-'•O - 0.6 

16 + 

'-'•^0 - 0.16 


16 + 

- 1.6 
2 + 

'^■^ - 0.2 

04 + 

^■^^ - 0.04 



Table 11.7: Expected yields for Z+jets, single top and diboson production with an additional 
prompt photon in the final state for 1.04 fb~^ from MC simulations. 



11.5 Summary of background processes with prompt photons 



Tab. 11. 8| gives a summary of the estimates for background processes with prompt photons in 
the final state for 1.04 fb~^ including the total uncertainties. 



Process 


e+jets 


;U+jets 


Background ti'^ 


0.8 


+ 


1.1 

0.8 


1.3 


+ 


1.9 
1.3 


Multijet+7 


1.2 


+ 


1.6 
1.2 


0.3 


+ 


1.0 
0.3 


V7+jets+7 


1.8 


± 


0.7 


3.7 


± 


1.4 


Z+jets+7 


1.3 


+ 


2.4 
1.3 


1.6 


+ 


2.3 
1.6 


Single top+7 


0.6 


+ 


0.7 
0.6 


0.2 


+ 


0.3 
0.2 


Diboson+7 


0.16 


+ 


0.34 
0.16 


0.04 


+ 


0.18 
0.04 



Table 11.8: Expected yields for background processes with prompt photons in the final state 
for 1.04 fb"^ 
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12 Systematic uncertainties 



Various effects may lead to systematic biases of the measurement of the tt'j cross section. The 
different sources of systematic uncertainties are discussed in the fohowing: Sec. 12.1 describes the 
treatment of systematic uncertainties due to the modelhng of the tt'j signal. Sec. 12.2 discusses 
systematic uncertainties arising from the modelling of the various background processes. Effects 
due to the modelling of the detector are presented in Sec. 12.3 and Sec. 12.4 discusses the 
uncertainty due to the measurement of the integrated luminosity. The different sources of 
systematic uncertainties were combined, as described in Sec. 12. 5[ 

Systematic effects may change the acceptance and efficiency of the ttj signal and the ex- 
pectations for the background contributions with misidentified electrons or prompt photons. 
They may also influence the shape of the p^^'^'^^ distributions used for the estimation of the 
background from hadrons misidentified as photons. 

For the evaluation of the systematic uncertainties, ensembles of pseudo-data were created 
according to the expected yields for ttj and the different background contributions. The ex- 
pected yields for the ttj processes were taken from the WHIZARD MC simulations (Sec. 6.3). 
The expectations for the background processes with electrons misidentified as photons and with 
prompt photons were derived in Ch. 10 and 11, respectively. For the contribution from hadrons 
misidentified as photons, the result of the final template fit (Ch. 13 ) was taken as expectation, 
because this was the most reliable estimate available. 

For each systematic uncertainty, 5000 ensembles of pseudo-data were created by fiuctuating 
each bin in the p^^^"^^ distribution according to the Poissonian uncertainty around the expec- 
tation. The pseudo-data were then fitted with modified efficiencies, background predictions 
and modified template shapes where applicable. The mean difference with respect to ensemble 
tests with the nominal efficiencies, background predictions and template shapes was taken as 
a measure of the systematic uncertainty. The intrinsic statistical uncertainty of the tests with 
5000 ensembles was found to be less than 0.4%, and hence negligible. Sources of systematic 
uncertainties which were found to yield an uncertainty smaller than the intrinsic uncertainty 
from the ensemble test, were conservatively estimated to be as large as 0.4%. 

All systematic uncertainties were symmetrised by taking the largest systematic variation as a 
conservative estimate of the upper and lower uncertainty. 



12.1 Signal modelling 

The choice of the MC generator, the impact of including higher-order QCD effects, as well as 
the choice of the parton showering model, the impact of the ISR and FSR modelling, and un- 
certainties on the choice of the PDFs were evaluated as sources of systematic uncertainties for 
the modelling of the ttj signal. All of these effects have an impact on the signal acceptance and 
efficiency only, and not on shape of the p^'^^'^^ templates. Since samples generated with alterna- 
tive MC generators, a different parton showering and varied settings for ISR and FSR were only 



available for tt production (Sec. 4.2) and not for ttj production, these effects were estimated by 



comparing the changes in acceptance and efficiency for tt events for the preselection (Sec. 6.1). 
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The relative change found in tt production was assumed to be an estimate of the relative sys- 
tematic uncertainty on the ti'j acceptance and efficiency. Similarly, the relative uncertainty due 
to the PDF choice in tt events was taken as an estimate for the relative uncertainty on the tt'j 
signal. 

Additionally, the dependence of the signal acceptance and efficiency on the pile-up conditions 
was studied. 



Monte Carlo generator: The acceptances and efficiencies for the preselection from tt samples 
generated with MC@NLO and POWHEG were compared. Both generators were interfaced to 
HERWIG for parton showering. The relative difference was taken as systematic uncertainty on 
the choice of the MC generator for ti'j production. 



Finite order calculation: The acceptances and efficiencies for the preselection from tt sam- 
ples generated with POWHEG and AcerMC were compared. While POWHEG is a NLO gener- 
ator, AcerMC generates ti production in LO. Both generators were interfaced to HERWIG for 
parton showering. The relative difference was taken as systematic uncertainty on the dependence 
of the ti'j modelling, generated in LO, on higher-order QCD corrections. 



Parton shower: The acceptances and efficiencies for the preselection from two tt samples 
generated with POWHEG were compared, where the HERWIG cluster model and the PYTHIA 
Lund string model for parton showering were compared. The relative difference was taken as 
systematic uncertainty due to the choice of the parton shower model. 



Initial and final state radiation: The acceptances and efficiencies for the preselection from 
different tt samples generated with AcerMC interfaced to PYTHIA for parton showering were 
compared, where the parameters which describe the amount of ISR and FSR in PYTHIA were 
varied in a range comparable to the Perugia Soft/Hard tune variations jl23 124 . The largest 



deviation with respect to the nominal AcerMC sample was taken as systematic uncertainty on 
the modelling of ISR/FSR in ttj production. 



Parton density functions: The impact of the choice of the set of PDFs was found to be 
of the order of 2% in previous cross section measurements in the single lepton tt decay chan- 
nel [136j . The procedure applied for the evaluation of the systematic uncertainty was based on 
a reweighting of the events in the MC@NLO tt MC sample according to different NLO PDF 



sets as described in [123 . Half of the total envelope of all variations was taken as symmetric 
systematic uncertainty. 



Since this uncertainty is small with respect to the contributions from other sources (Ch. 13), 
it was not evaluated separately for this analysis. Twice the uncertainty found in [136| was taken 
as a conservative estimate. 



Pile-up conditions: The stability of the signal modelling with respect to the pile-up condi- 
tions was checked by studying a possible dependence of the selection efficiency of ti'y simulations 
on the number of interactions per bunch crossing. This is shown in Fig. [12. 1| for the electron 
(left) and for the muon channel (right). The selection efficiency increases slightly with increas- 
ing pile-up. However, as shown in Fig. |3.2[ nearly all data were taken with an average number 



of interactions between three and eight as indicated by the grey regions in Fig. 12.1, In this 



regime, a stable selection efficiency was observed in the simulations, which is illustrated by the 
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linear fits in Fig. |12.1 Also, MC simulations were reweighted to the pile-up conditions present 
in data, so that pile-up effects were safely ignored as a source of systematic uncertainty on the 
signal modelling in this analysis. 
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Figure 12.1: Total selection efficiency of simulated ti^ events as a function of the number of 
interactions per bunch crossing in the electron channel (left) and in the muon 
channel (right). The shaded regions indicate the number of interactions per 
bunch crossing present in the data analysed in this thesis. Additionally, a linear 
fit to these regions is shown. 



12.2 Background modelling 

The number of events from processes with hadrons misidentified as photons was estimated using 
a fit to the p^j?"*^^*^ distribution of the photon candidates (Ch. [7|. Hence, effects which change 
the shape of either the prompt photon template or the hadron fake template were considered as 
sources of systematic uncertainties. 

Moreover, uncertainties on the background contributions with electrons misidentified as pho- 
tons and on the backgrounds with prompt photons were considered. Some of the contributing 
processes were modelled with MC simulations, which are subject to several sources of systematic 
uncertainty: the cross sections for the respective processes are only known with a certain pre- 



cision (Sec. 4.2) and the luminosity measurement was subject to uncertainties (Sec. 12.4). The 



uncertainties on the cross sections read Ij^g^ for tt, 48% for Z+jets, and 5% for diboson produc- 
tion. For single top production, the uncertainties are as large as ±4% and ±7% in the t-, s- 



and l^t-channels, cf. Sec. 4.2 Additionally, uncertainties in the detector simulation (Sec. 12.3) 



affect all simulated processes. The luminosity and the detector modelling uncertainties were 



treated as correlated between all MC samples as described in Sec. 12.3 and 12.4 

Additionally, the effect of limited available MC statistics was evaluated for all background 
contributions estimated from MC simulations. The upper (lower) limit was estimated by the 
expectation value of a Poissonian distribution giving a probability smaller than 16% for observing 
less (more) events than estimated in the simulations. 

The uncertainties on the cross section predictions were also treated as correlated between the 
different processes which were estimated from the same MC sample, as for Z+jets+7 production 
and Z+jets production with an electron misidentified as a photon. Template fits were performed 
with cross section predictions varied within their respective uncertainties and the largest dif- 
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ference with respect to the nominal template fit was taken as systematic uncertainty on the fit 
result. 

pcone20 template for prompt photons: The derivation of the prompt photon template from 
electron distributions in Z — )■ ee events was discussed in Ch. [8| In particular, a MC correction 
for the differences between the electron and photon distributions was applied. The nominal 
photon template is shown together with the uncorrected electron template in the upper plot in 



Fig. 12. 2[ The uncertainty due to the MC correction was estimated by comparing template fits 



with the two different templates. The difference between the two fits was taken as systematic 
uncertainty. 

The shape of the prompt photon template for unconverted and converted photons only is 



shown in the lower left plot of Fig. 12.2|in comparison to the nominal template. The dependence 



on the fraction of converted photons is only due to the MC correction applied to the electron 
templates from Z — )• ee data. Since this fraction is unknown, the nominal fit was compared to 
fits with the unconverted- and converted-only templates. The largest difference with respect to 
the nominal fit was taken as systematic uncertainty. 



The number of primary vertices per event is a measure of the amount of pile- up. Fig. 12.3 
shows the distribution of the number of primary vertices for the selected events in the electron 
and the muon channel. Additionally, the distribution for the Z — )• ee sample which was used for 
the derivation of the electron templates is shown as well as the distribution for the control region 
(CR) used for the derivation of the hadron fake template as defined in Ch. |9j The distribution 
in the Z — )• ee sample is consistent with those in the selected ti'j candidate events in both lepton 
channels given the limited statistics, and hence the prompt photon template represents the same 
pile-up conditions as present in the candidate events. For the hadron fake CR, a slightly higher 
average pile-up is observed than in the Z — )• ee sample, which is discussed in the next section 
treating the hadron fake template. 

Nevertheless, the dependence of the prompt photon template on pile-up was studied and 



the lower right plot in Fig. 12.2 shows templates derived for a low pile- up regime with 1-5 
primary vertices per event and a high pile-up regime with 6-10 primary vertices. These two 
templates with extremely different pile-up conditions were tested in the template fit and the 
largest difference with respect to the nominal sample was taken as a very conservative estimate 
for the systematic uncertainty due to a mismodelling of the pile-up conditions. 

pcone20 template for hadrons misidentified as photons: As discussed in Ch. |9| the tem- 
plate for hadrons which were misidentified as photons showed a dependence on the object's 
and 7], and the templates were reweighted accordingly. 

For the reweighting in Et, an exponential fit to the distribution in a CR was used, which 
was subject to significant uncertainties due to the limited amount of hadron fake candidates in 
the CR. The upper left plot in Fig. |12.4l shows the nominal hadron fake template together with 
templates which were reweighted with exponential curves corresponding to variations of the mean 
lifetime within the uncertainties from the nominal exponential fit. An additional systematic 
uncertainty of 25% was added on the mean lifetime of the exponential due to differences between 
true hadron fakes and hadron fake candidates observed in ti simulations (Ch. |9]). The largest 
difference of the template fits with the two alternative templates with respect to the nominal fit 
was taken as systematic uncertainty. 

The reweighting in rj was done in two regions < |7/| < 1.81 and 1.81 < \rj\ < 2.37, where 
the fraction of high-|r/| hadron fakes was also estimated from the CR. The upper right plot in 



Fig. 12.4 shows the nominal hadron fake template together with templates which were reweighted 
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Figure 12.2: Comparison of the nominal prompt photon p^^^"^^ template (thin solid line) with 
systematically varied templates (thick solid and thin dashed lines) due to differ- 
ent effects: the difference between photon and electron templates (upper plot), 
the fraction of unconverted and converted photons (lower left plot), and the de- 
pendence of the p^^^'^^ distribution on the pile- up conditions (lower right plot). 
Below each plot, the difference between the systematically varied templates and 
the nominal template is shown. In all plots, the last bin includes the overflow 
bin. 
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Figure 1 2. 3: Normalised distribution of the number of primary vertices of the selected events 
in the electron channel (solid circles) and in the muon channel (open squares). 
Additionally, the distribution for the Z — )• ee sample which was used for the 
derivation of the electron templates is shown (thick solid line) as well as the 
distribution for the control region used for the derivation of the hadron fake 
template (thin dashed line). The last bin includes the overflow bin. 
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Figure 12.4: Comparison of the nominal p^j?"''^'' template for hadron fakes (thin solid line) with 
systematically varied templates (thick solid and thin dashed lines) due to different 
effects: the uncertainties on the derivation of the spectrum of the hadron fakes 
(upper left plot) and the amount of hadron fakes with large \ri\ (upper right plot), 
the fraction of unconverted and converted hadron fakes (lower left plot), and the 
dependence of the p^^^"^^ distribution on the pile-up conditions (lower right plot). 
Below each plot, the difference between the systematically varied templates and 
the nominal template is shown. In all plots, the last bin includes the overflow 
bin. 



according to varied fractions of high-|T/| hadron fakes within the uncertainties derived in Ch. [oj 
The largest difference of the template fits with the two alternative templates with respect to the 
nominal fit was taken as systematic uncertainty. 

As for the prompt photon templates, the fraction of unconverted and converted photons from 



misidentified hadrons was unknown. The lower left plot in Fig. 12.4 shows the templates for 
unconverted and converted hadron fakes only in comparison to the nominal template. The 
nominal fit was compared to fits with the unconverted- and converted-only templates, and the 
largest difference with respect to the nominal fit was taken as systematic uncertainty. 

The average number of primary vertices is slightly larger for the events in the CR used for the 
derivation of the hadron fake template compared to the Z — )• ee samples from which the prompt 



photon template was derived (Fig. 12.3). This is expected, since with increasing pile-up, the 



number of events with at least one photon which is faked by a hadron from fragmentation may 
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increase because of the increasing number of jets, while the number of Z events is not expected 
to increase significantly. 

Although the pile-up conditions are in good agreement with those present in the ti'^ candidate 
events in both lepton channels, the dependence of the hadron fake template on the amount of pile- 
up was studied. The lower right plot in Fig. 12.4 shows templates derived for a low pile-up regime 
with 1-5 primary vertices per event and a high pile-up regime with 6-10 primary vertices. 
These two templates with extremely different pile-up conditions were tested in the template fit 
and the largest difference with respect to the nominal sample was taken as a very conservative 
estimate for the systematic uncertainty due to a mismodelling of the pile-up conditions. 



Estimate of background contributions with electrons misidentified as photons: The 

contribution from processes with electrons which were misidentified as photons was derived in 
Ch. 10 Li particular, scale factors (SFs) for the misidentification rate were derived in order to 
correct MC simulations. The resulting yields were presented in Sec. 10.2 together with systematic 
uncertainties. 

The SFs for the electron-to-photon misidentification rate were varied within their uncertain- 
ties. Moreover, the background estimate was varied within the uncertainty due to the limited 
MC statistics. In both cases, the largest difference with respect to the nominal template fit 
result was taken as systematic uncertainty. 

Uncertainties due to detector modelling, the predictions for the cross sections of the different 
processes, and due to the luminosity measurement, were treated separately in order to account 
for the correlations between the processes estimated from MC simulations, as stated above. 



Estimate of background contributions with prompt photons: The contribution from 
background processes with prompt photons was discussed in Ch. 11 The systematic uncertain- 
ties on the predictions were found to be sizable: they are at least as large as 100% for processes 
which were believed to be only modelled with limited precision in MC simulations. This holds 
for tt^ production outside of the signal phase space as well as for backgrounds from Z+jets+7, 
single top+7 and diboson+7 production. The limited MC statistics was taken into account in 
the uncertainties. 

The estimates for Ty+jets+7 production and for multijet processes with prompt photons 
were derived with data-driven techniques involving template fits to the p™"*^^*^ distribution of 
the photon candidates in the respective CRs. Both estimates feature large uncertainties from the 
limited statistics in the CRs. The sources of systematic uncertainties on the prompt photon and 
hadron fake templates as described above were evaluated also in these template fits. Additionally, 
the M^+jets+7 estimate was found to be subject to uncertainties from the extrapolation to the 
signal region. The multijet+7 estimate features large uncertainties due to the amount of fake 
leptons. 

In order to estimate the systematic uncertainty from backgrounds with prompt photons, the 
estimate of each contributing process was varied within its uncertainty in the template fit. The 
largest deviation with respect to the nominal fit result was taken as systematic uncertainty in 
each case. 

As stated above, uncertainties on the cross sections, the luminosity, and the detector modelling 
were treated separately in order to account for correlations between the processes estimated from 
MC simulations. 
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12.3 Detector modelling 

The MC simulations used in this analysis included a detailed simulation of the ATLAS detec- 
tor. The modelling of the detector is subject to systematic uncertainties which translate into 
uncertainties on the description of the different physics objects: uncertainties on the electron, 
muon, jet and photon modelling were evaluated as well as uncertainties on the description of 
the measurement and the 6-tagging performance. An additional systematic uncertainty was 
evaluated for the modelling of the broken optical links in the LAr calorimeter (Ch. [s]). 

Systematic uncertainties on the detector modelling do not influence the shape of the p^^^"^^ 
templates, but only change the respective acceptances and efficiencies of the signal and back- 
ground processes. They were evaluated for all processes estimated from MC simulations, in- 
cluding the ti'j signal as well as the background MC simulations. The variations were treated 
consistently between all processes in order to account for the correlations between the samples 
due to the different systematic effects. 

Electron modelling: MC simulations were corrected by SFs for the mismodelling of elec- 
tron trigger, reconstruction and identification efficiencies, which were derived from Z — >• ee and 
W ^ eu data with tag-and-probe methods [103| (Sec. 5.1). The SFs were varied within their 



respective uncertainties and the largest difference obtained with respect to the nominal template 
fit was considered as systematic uncertainty. 



The electron energy scale was corrected in data based on measurements in Z — t- ee events 103 
In order to evaluate the associated systematic uncertainty, the electron scale in MC was varied 
within the relative scale uncertainty. The electron energy resolution was also measured in Z — )• ee 



events [103 , and the resolution in MC simulations was corrected to match the resolution observed 



in data. The systematic uncertainty was estimated by comparing template fits with larger and 
smaller energy resolutions according to the uncertainty on the resolution measurement. 

The largest deviation with respect to the nominal template fit was taken as systematic un- 
certainty for the energy scale as well as for the energy resolution. 

Muon modelling: As for electrons, also for muons, SFs for the mismodelling of the trigger. 



reconstruction and identification efficiencies were applied to MC simulations [138J (Sec. 5.2 ). The 
SFs were derived with tag-and-probe methods from Z — )• events. They were varied within 

their respective uncertainties and the largest difference obtained with respect to the nominal fit 
was considered as systematic uncertainty. Since trigger matching was not applied in data in the 
muon channel, an additional systematic uncertainty was considered, which was estimated to be 
as large as 2.2% for events with zero or one 6-tags, and 1.5% for events with at least two 6-tags. 
The muon momentum scale and resolution were corrected in MC simulations based on mea- 



surements in Z — )• events 139 . In order to evaluate the associated systematic uncertain- 

ties, the template fit was performed with varied scales and resolutions within their respective 
uncertainties. For both effects, the largest deviation with respect to the nominal template fit 
was taken as systematic uncertainty. 

Jet modelling: The jet energy scale (JES) was derived using data from test beams and LHC 
collisions, and from simulations [104 | (Sec. |5.3[ ). The total JES uncertainty is of the order of 
4% (6%) for low-pT jets in the range 0.3 < \t]\ < 0.8 (2.1 < \r]\ < 2.8) and decreases to roughly 



2% (2.5%) for jets with 60 GeV < px < 800 GeV (Fig. |2.11[ ). Additional uncertainties were 
added due to the different flavour composition and the higher amount of close- by jets present 
in tt events with respect to the dijet topologies used for the derivation of the JES. In order to 
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account for a possible pile-up dependence, an uncertainty was considered of up to 5% in the 
central and up to 7% in the forward region {\rj\ > 2.1). Jets originating from b-quarks may have 
a slightly different energy scale than light quark and gluon jets. This was accounted for by an 
additional uncertainty on the 6-JES, as large as 2.5% for a of 20 GeV, and decreasing to 
0.0076% for a pT of 600 GeV. The uncertainties on the JES, its pile-up dependence and the 
JES for 6-jets were varied within their respective uncertainties and for each effect, the largest 
difference with respect to the nominal fit was considered as systematic uncertainty. 

The jet energy resolution ( JER) was measured in dijet events [143] and was found to agree with 
MC simulations within 10%. In order to account for a possible bias due to the JER modelling, the 
resolution was worsened in the simulations according to the uncertainty of the JER measurement. 
The difference with respect to the nominal fit was symmetrised and considered as systematic 
uncertainty. 

The jet reconstruction efficiency was measured by comparing jets built from calorimeter clus- 



ters to jets built from tracks 104 . Data and simulations were found to be in agreement within 
2%. The associated systematic uncertainty was estimated by randomly disregarding jets in MC 
simulations based on the uncertainty on the reconstruction efficiency. The difference of the 
template fit with respect to the nominal fit was taken as systematic uncertainty. 

Photon modelling: The photon identification efficiency was measured by shifting the MC 



shower shapes so that they matched the shower shapes in data 151 153] (Sec. 5.6). The as- 



sociated relative systematic uncertainties were taken from the diphoton cross section measure- 



ment [151 . The evaluation of the uncertainties included the intrinsic precision of the method, 
the choice of the photon candidate sample, the knowledge of the material in front of the calorime- 
ter, the admixture of fragmentation photons and the classification of unconverted and converted 
photon candidates. The efficiency was varied within its uncertainty and the largest difference 
with respect to the nominal result was taken as systematic uncertainty. 

Differences in reconstruction efficiencies between data and simulations were found to be smaller 
than 1% for electrons with uncertainties smaller than 0.7% for |r/| < 2.37. Since photons were 



reconstructed within the same algorithm (Sec. 5.6), differences in the photon reconstruction 



efficiency between data and simulations can be expected to be of the same order of magnitude. 
The uncertainties due to the photon identification efficiency are significantly larger, and the 
additional uncertainty due to the mismodelling of the reconstruction efficiencies was believed to 
be covered by the uncertainty on the identification efficiencies. 

The photon energy scale in data was corrected based on the electron scale measurement in 
Z ^ ee events [103 . In order to evaluate the associated systematic uncertainty, the scale in 



MC simulations was varied within the relative scale uncertainty. The electron energy resolution 
calibration was used to correct also the photon energy resolution in MC simulations. The 
systematic uncertainty was estimated by comparing template fits with larger and smaller energy 
resolutions according to the uncertainty on the resolution measurement. The largest deviation 
with respect to the nominal template fit was taken as systematic uncertainty due to the energy 
scale as well as due to the energy resolution, respectively. 

modelling: All systematic variations of the energy scales and resolutions of leptons, jets 
and photons were propagated consistently into the calculation of the missing transverse energy. 
Additional systematic uncertainties were due to a possible mismodelling of the soft jet and cell 



out terms in the cf- Sec. 5.4 and Eq. (5.1). Moreover, a 10% uncertainty was added in 



order to cover a possible pile-up dependence of the measurement not modelled in the MC 
simulations. All of these sources of systematic uncertainties were treated as correlated. 
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For each effect, the measurement was varied within the respective uncertainties and the 
largest difference with respect to the nominal template fit was considered as systematic uncer- 
tainty. 



Modelling of the b-tagging performance: The 6-tagging efficiencies and mistag rates were 
measured in data 147, 148] (Sec. 5.5), and SFs for the correction of the MC simulations were 
derived. The uncertainties on the SFs vary between 0.06 and 0.15 for the 5-tagging efficiencies 
and between 0.11 and 0.22 for the mistag rates. Template fits with varied efficiencies according 
to the uncertainties on the SFs were compared to the nominal fit result, and the largest deviation 
was taken as systematic uncertainty. 



Modelling of LAr readout issues: A large fraction of the data analysed featured a broken 
optical link in the LAr barrel calorimeter (Ch. [3]). In data, events with a jet with a of at least 
20 GeV which was closer to the affected calorimeter region were disregarded. In simulations, 
events with such a jet were reweighted according to the fraction of the integrated luminosity 
which featured the broken optical link. The systematic uncertainty on this procedure was 
evaluated by varying the px-threshold in MC simulations by it 4 GeV. Template fits with 
these variations were compared to the nominal fit result, and the largest deviation was taken as 
systematic uncertainty. 



12.4 Luminosity measurement 



The luminosity of the LHC was measured with van der Meer scans |155| , and the uncertainty 
on the luminosity measurement was estimated to 3.7% 156 ,157 . It was considered as common 



systematic uncertainty to all signal and background processes modelled with MC simulations, 
and was treated as correlated between all of these contributions. 



12.5 Combination of systematic uncertainties 

The systematic uncertainties due to signal, background and detector modelling were assumed to 
be uncorrelated and were hence summed in quadrature. This holds also for the uncertainty on 
the luminosity measurement. The only exception is the combination of the uncertainty on the 
pile-up dependence of the JES and the uncertainties, which were treated fully correlated. 
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In this chapter, the results of the final template fit to the p™"*^^^ distribution of the photon 
candidates are presented. The results were translated into a tt'j cross section estimate. Moreover, 
the significance of the measured tt'j signal was estimated. At the end of this chapter, the 
evolution of the expected significance with the integrated luminosity is estimated and discussed. 

In order to estimate the amount of tt^ signal events and of background events with hadrons 
misidentified as photons, the template fit described in Ch. [7]was performed including the back- 
ground estimates for processes with prompt photons and electrons misidentified as photons 
(Ch. 10 and 11 ). Fig. 13.1 shows the result of the combined fit in both lepton channels together: 
the left plot shows the resulting distributions in the electron channel, and the right plot shows 
the distributions in the muon channel. 
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Figure 13.1: Result of the combined template fit to the photon p™"^^'' distribution in both 
lepton channels. The left plot shows the distributions in the electron channel, 
the right plot shows the distributions in the muon channel. In both plots, the 
last bin includes the overflow bin. 



The signal and th e diff erent background contributions as depicted in Fig. 13.1, are shown in 
an overview in Tab. 



13.1 



-11 8 events out of the 122 observed candidate 



The fit yielded 43 

events. The uncertainties quoted in the table are the statistical uncertainties from the fit only. 
While the background contributions with prompt photons and electrons misidentified as photons 
were fixed in the fit, for the signal contribution and the contributions with hadrons misidentified 
as photons, the smallest interval containing 68% of the marginalised probability density function 
is quoted. 

In the combined fit, the efliciency with respect to the whole tt'j signal phase space was ac- 
counted for, as described in Sec. [Tj The fit yielded 1930 ^ggj] ti'y events in the whole signal 
phase space before object and event selection including true photon transverse momenta down 
to 8 GeV. The uncertainty was again extracted from the smallest interval containing 68% of the 
marginalised pdf, which is shown in Fig. 13.2 The dashed line shows the global mode of the fit. 



103 



13 Results 



Contribution 


Fit result / yield 


tt'j signal 

Background with prompt photons (e+jets) 

Background with prompt photons (/x+jets) 

Background with electrons misidentihed as photons (e+jetsj 

Background with electrons misidentihed as photons (^u+jets) 

Background with hadrons misidentihed as photons (e+jets) 

Background with hadrons misidentihed as photons (^+jets) 


^'^■^ - 11.8 

5.9 
7.1 

O '7 

8.7 
10.3 

20.5 + li 

26.6 t e:? 


Data (e+jets and ;U+jets) 


122 



Table 13.1: Result of the combined template fit showing the different contributions together 
with the statistical uncertainty for the contributions which were varied in the fit. 
The amount of observed candidate events in data is also shown. 
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Figure 13.2: Posterior probability density for the number of tfy events for the combined tem- 
plate fit. The dashed line shows the global mode of the fit. The marked region 
around the mode represents the smallest interval containing 68% of the distribu- 
tion. 




As a cross-check, the template fit was performed in the electron channel and muon channel 
separately. The results of these fits are shown in Fig. |13.3 for the two channels. They yielded 
1640 ^'^y^Q and 1990il^y2o events in the whole signal phase space, respectively. The fits in the 
single channels are hence consistent with each other within statistical uncertainties. 



Systematic uncertainties were evaluated as described in Ch. 12 An overview of the systematic 



uncertainties is presented in Tab. 13.2 in terms of fraction of the measured number of tt^ 
events. The different sources of systematic uncertainties were summarised in several categories 
as indicated in the table. 

The largest systematic uncertainty is due to the jet and modelling (30%), which is dom- 
inated by the nominal uncertainty on the jet energy scale (JES), and the additional JES un- 
certainty due to a possible pile-up dependence of the JES. Both sources of uncertainty change 
the acceptance of the tt7 signal by roughly ±10%. The acceptance changes for the background 
processes are of the same order, or even larger, as for example for Z+jets events, which feature 
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Source 


Uncertainty 


Source 


Uncertainty 


MC generator 

Finite order calculation 

Showering 

Initial / final state radiation 

Parton density functions 


± 4.8 % 
± 3.2 % 
± 3.1 % 
± 8.0 % 
± 4.0 % 


Signal modelling 


± 11.1 % 


MC correction (plioton template) 
Converted photons (photon template) 
Pile-up (photon template) 
?7-reweighting (hadron fake template) 
i?T-reweighting (hadron fake template) 
Converted photons (hadron fake template) 
Pile-up (hadron fake template) 


± 6.4 % 
± 2.6 % 
± 0.5 % 
± 2.5 % 
± 7.4 % 
± 7.3 % 
± 2.0 % 


Template shapes 


± 12.9 % 


Electron trigger efficiency 
Electron reconstruction efficiency 
Electron identification efficiency 
Electron energy scale 
Electron energy resolution 


± 1.0 % 
± 1.4 % 
± 4.5 % 
± 1.5 % 
± 1.7 % 


Electron modelling 


± 5.3 % 


Muon trigger efficiency 

Muon reconstruction efficiency 
Muon identification efficiency 
Muon momentum scale 
Muon momentum resolution 


± 3.2 % 
± 0.5 % 
± 0.4 % 
± 0.4 % 
± 0.7 % 


Muon modelling 


± 3.4 % 


Jet energy scale 
6-jet energy scale 
Jet energy scale and (pile-up) 
Jet energy resolution 
Jet reconstruction efficiency 
(cell-out contribution) 


± 18.4 % 
± 2.8 % 
± 23.3 % 
± 4.2 % 
± 0.4 % 
± 1.9 % 


Jet and modelling 


± 30.2 % 


Photon identification efficiency 
Photon energy scale 
Photon energy resolution 


± 4.8 % 
± 1.4 % 
± 2.2 % 


Photon modelling 


± 5.5 % 


6-tagging performance 


± 12.5 % 


^-tagging 


± 12.5 % 


Modelling of LAr readout issues 


± 2.4 % 


LAr readout issues 


± 2.4 % 


/e-»-y scale factor 


± 4.4 % 


fe^j scale factor 


± 4.4 % 


Background tt^ estimate 

Z-|-jets-|-7 estimate 
Single top-|-7 estimate 
Diboson-1-7 estimate 
Multijet-l-7 (fake leptons) 
Multijet+7 (prompt photons) 
H^+jcts+7 estimate 


± 7.1 % 
± 5.9 % 
± 2.1 % 
± 1.0 % 
± 1.1 % 
± 5.0 % 
± 1.2 % 


Prompt photon background 


± 10.9 % 


ttj fc-factor and MC statistics 

tt cross section and MC statistics 
Z+jets cross section and MC statistics 
Single top cross section and MC statistics 
Diboson cross section and MC statistics 


± 0.6 % 
± 4.4 % 
± 9.1 % 
± 0.4 % 
± 0.6 % 


Background cross sections 
and MC statistics 


± 10.2 % 


Total systematic 


± 40.9 % 


Luminosity 


± 6.1 % 


Statistical 


+25.1% 
-27.6% 


Total 


+48.3% 
-49.6% 



Table 13.2: Statistical and systematic uncertainties from the combined fit. 
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Figure 1 3. 3: Result of the two separate template fits to the photon p'^'^^'^^ distribution in the 
electron channel (left) and in the muon channel (right). In both plots, the last 
bin includes the overflow bin. 



a softer jet pT spectrum, i.e. more jets close to the threshold. The resulting uncertainties 
on the cross section due to the JES and the JES pile- up dependence are of the order of 20%, 
which can be understood as follows: a smaller (larger) background prediction results in a larger 
(smaller) number of estimated signal events. Moreover, the change of the tt'j acceptance means 
that the measured amount of signal corresponds to a smaller (larger) phase space and therefore 
is subject to a larger (smaller) acceptance correction. The acceptance variations of the signal 
and the background modelling hence have a positive correlation and lead to the observed large 
systematic uncertainties on the cross section of the order of 20%. 

In contrast, the modelling of electrons, muons and photons results in only moderate systematic 
uncertainties (3.4 - 5.5%). The largest contributions are the electron and photon identification 
efficiencies and the muon trigger efficiency, where the latter was increased to cover effects due 
to trigger matching (Sec. 12. 3| ). The uncertainty due to the modelling of the 6-tagging perfor- 
mance was found to be sizable (13%). Uncertainties due to the modelling of the readout issues 
in the LAr calorimeter and due to the systematic uncertainty on the SF for the electron-to- 
photon misidentification rate were found to have only moderate contributions (2.4% and 4.4%, 
respectively) . 

The change of the ttj acceptance only is of the order of 5% for both the variations of the 
photon identification efficiency and of the 6-tagging performance. In the case of the 6-tagging 
performance, the background estimate also changes significantly, which leads to an uncertainty 
on the cross section of the order of 10% due to the positive correlation as discussed for the JES 
uncertainty. For the photon identification efficiency, however, the background estimate only 
varies by a small amount: the main background estimated from MC simulations are events with 
electrons misidentified as photons. Although the estimate for these background contributions is 
expected to increase (decrease) for a larger (smaller) photon identification SF, it does not vary 
significantly, because the electron-to-photon misidentification rate /e^7 changes in the opposite 
direction. This is due to the fact that the number of events in the Z — >■ e7fake sample in MC 
simulations appears in the denominator in Eq. (10.1). The uncertainty on the cross section is 
therefore only of the order of 5%, and hence smaller than the uncertainty due to the 6-tagging 
modelling. 

Uncertainties due to the background modelling arise from the shape of the templates (13%), 
the estimates of the background contributions with prompt photons (11%) and the uncertainties 
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due to the predicted cross sections and limited MC statistics for various background processes 
(10%). The uncertainty due to the template shapes was found to be dominated by the MC 
correction applied to the electron templates from Z — )• ee data and the fractions of unconverted 
and converted photon candidates used for the hadron fake template, as well as the spectrum 
of the hadron fakes. The first two uncertainties were estimated conservatively: the uncertainty 
on the MC correction was estimated by comparing the result with the correction to the result 
without any correction. For the fractions of unconverted and converted photon candidates, the 
nominal fit was compared to fits in which it was assumed that all hadron fakes were either 
unconverted or converted, respectively. The uncertainty due to the hadron fake spectrum, 
in turn, is dominated by the low statistics in the control region in which the spectrum was 
estimated (Ch. [9]). 

The dominant contributions to the uncertainty from the prompt photon background were 
found to be due to the data-driven multijet-l-7 background and the estimate of the contribution 
from Z-|-jets-|-7 events, but also due to the estimate of the amount of background ttj events. All 
of these estimates feature conservative uncertainties larger than 100% (Sec. |11.T| |11.2 and 11.4). 



The uncertainties on the tt and Z+jets cross sections are the dominant uncertainties among 
the predictions of cross sections for the background processes as well as the limited statistics of 
the corresponding MC samples: dileptonic tt events with an electron misidentified as a photon 



result in a large background contribution (Sec. 10.2 ). For the Z+jets process, the same simulated 



sample was used to estimate the contributions from electrons misidentified as photons as well 
as from Z-l-jets-l-7 events. Although the resulting background yields are not very large, the 
uncertainty on the Z+jets cross section still has a non- negligible effect, because it is only known 
to 48%. A very small uncertainty was found to be due to the uncertainty of the /c-factor for ti'j 
production, which originates from a small corrections from tt'j events in the control region used 
for the estimate of the VF-l-jets-l-7 background. 

Uncertainties due to the modelling of the ttj signal add up to 11%, where the largest contri- 
bution is from the modelling of the initial and final state radiation. All systematic uncertainties 
add up to a total systematic uncertainty of 41%. 

From the measured number of tt'j events, the tt'j cross section in the single lepton and dilepton 
channels was calculated using the integrated luminosity of 1.04 fb~^ of the data set analysed. 
The uncertainty on the tt'j cross section due to the luminosity measurement was found to be as 
large as 6.1%. It is larger than the intrinsic uncertainty on the luminosity measurement because 
of the positive correlation between the effects on the background estimates and on the number 
of tt^ signal events. The measurement of the tt'j cross section times branching ratio (BR) into 
the single lepton and dilepton final state reads: 

au ■ BR = 1.9 ± 0.5 (stat.) ± 0.8 (syst.) ± 0.1 (lumi.) pb . 



This result is compatible with the SM expectation (Sec. 4.1 ) of 2.1 it 0.4 pb within uncertainties. 

In order to estimate the significance of the measured ttj signal, the probability of the back- 
ground giving rise to at least 122 observed candidate events was estimated. The total back- 
ground yield in both lepton channels together reads 79l"'^g (stat.) it 13 (syst. + lumi.) = 79lj;g. 
The total uncertainty was interpreted as the standard deviation of a Gaussian pdf. The total 
background pdf was constructed by a convolution of this Gaussian pdf with a Poissonian pdf 
to allow for statistical fluctuations. This pdf is shown in the left plot of Fig. |13.4[ The p- value 
for the background to produce at least 122 observed events is 1.4%, which corresponds to a 
significance of 2.5 a. 
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The right plot in Fig. 13.4 shows the expected significance given the background estimate 
as well as the ti'j expectation from the SM calculation, which reads 49 ± 7(stat.) it 10 (syst.). 
The systematic uncertainty is due to the uncertainty on the ttj /c-factor. The mean of the 
distribution is 2.8 a. However, the root mean square is as large as 1.0 o" given the statistical and 
systematic uncertainties on the predicted number of tt^ events. The observed significance of 
2.5 a is hence consistent with the expected significance within this uncertainty. 



0.02- 



0.01- 





— bkg. only 




p-value = 1 .4 % 




ATLAS work in progress 




|Ldt= 1.04 fb"' 



50 100 150 200 250 

number of observed events 



0.15 



0.1 



0.05 



mean = 2.8± 1.0o 

expected significance 

measured significance 



ATLAS work in progress 
|Ldt=1.04fb"' 



6 8 10 

expected significance [a] 



Figure 13-4: Estimation of the significance: the left plot shows the probability distribution 
of the observed number of events for the background-only hypothesis given the 
statistical and systematic uncertainties. The p-value for observing at least 122 
events is 1.4%, which corresponds to a significance of 2.5 a. The right plot shows 
the distribution for the expected significance. The solid line shows the mean of the 
distribution at 2.8 a. The root mean square is 1.0 a, as indicated by the marked 
region around the mean. The dashed line shows the measured significance. 



Evolution of the expected significance with increasing integrated luminosity 

It is interesting to investigate how the expected significance of the ti'j signal will evolve when data 



sets increase in size. Fig. 13.5 shows the evolution of the expected significance for integrated 
luminosities between 1 and 15 fb"^. Four different scenarios are shown: the thick solid line 
represents the systematic uncertainties as presented in this analysis. The thin solid line shows 
the significance for a scenario where the uncertainty due to the JES modelling was reduced by a 
factor of two. This is a likely scenario for the near future corresponding to a JES uncertainty of 
the order of 1%. The dotted and dashed-dotted lines show the significance for scenarios where 
the contributions from background processes were reduced by 30% and 40%, respectively and 
the relative systematic uncertainty on the background estimations was assumed to be the same 



as in this analysis. In Fig. 13.5, the uncertainty band is shown for the first scenario only. The 



uncertainty corresponds to the root mean square of the distribution of the expected significance 



(see right plot in Fig. 13.4), which is large and mainly due to the uncertainty on the expected 
tt'j cross section. 

The third and fourth scenarios assume a significant reduction of the background contributions, 
while the main background processes feature hadrons or electrons misidentified as photons. 
Background contributions from hadrons misidentified as photons could be further reduced by 
tightening the requirements on the photon shower shapes or by combining several shower shapes 
in a multivariate discriminant exploiting their correlations. The background contribution from 
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Figure 13.5: Evolution of the expected significance with the integrated luminosities for four 
different scenarios: systematic uncertainties as in this analysis (thick solid line), 
uncertainties on the jet energy scale reduced by a factor of two (thin solid line), 
and with all background contributions reduced by a factor of 30% and 40% (dot- 
ted and dashed-dotted lines). The uncertainty band is only shown for the first 
scenario. It corresponds to the root mean square of the distribution of the ex- 
pected significance. 

The horizontal dashed lines indicate the 3 a threshold for an evidence of ti'^ 
production and the Sex threshold for an observation. The vertical solid line 
indicates the maximum available integrated luminosity at ^/s = 7 TeV of 5 fb~^. 



electrons misidentified as photons could be reduced by considering only unconverted photon 
candidates, because a large fraction of the misidentified electrons results in converted photon 
candidates. This requirement would reduce the efficiency for photons from ti^ production by 
roughly 50%, but would clearly enhance the signal-over-background ratio. 

Generally, the expected significance increases with an increasing integrated luminosity. At 
5 fb~^, which is the size of the whole data set available at 7 TeV, the measurement is largely 
dominated by systematic uncertainties. A significance of more than 3 a, corresponding to an 
evidence of the ify signal, can be expected for systematic uncertainties as presented in this thesis. 
Deviations from the predicted A;-factor for ti'^ production may however lead to significantly 
smaller or larger values of the significance. From a decreased JES uncertainties, an increase of 
the significance of roughly 0.5 a can be expected. 

A significance of 5 cr, corresponding to an observation of the tfy signal, will only be possible 
if the contribution from background processes can be significantly reduced. With a reduction 
of 30 - 40%, a 5o" observation is in reach with 5 fb~^ if the relative systematic uncertainties 
remain of the same order and the signal efficiency is not strongly reduced. 

Data sets larger than 5 fb~^ will only be available at ^/s = 8 TeV or at larger energies. The 
cross sections for the different processes change with respect to 7 TeV, but to first order the 
extrapolations shown in Fig. 13.5 are expected to hold also for 8 TeV. At 8 TeV, the detector 
modelling will need to be studied, and a reduction of systematic uncertainties will only be 
possible after thorough studies of the performance of the different objects. The increasing pile- 
up may will be a particular issue. 
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Top quark pair production with an additional photon in the final state (^7) is a process which 
provides direct sensitivity to the strength and structure of the top quark's electromagnetic 
coupling. Since the top quark decays before producing bound states, ttj production is a unique 
process for the direct measurement of the electromagnetic coupling of a quark and hence an 
important test of the Standard Model. Additionally, tt'y production is sensitive to modifications 
of the t7-vertex beyond its Standard Model structure. 

In this thesis, the first measurement of the ttj cross section at the LHC was presented using 
1.04 fb~^ of data taken with the ATLAS detector in proton-proton collisions at a centre-of-mass 
energy of ^/s = 7 TeV. It is hence an important milestone towards tests of the electromag- 
netic coupling of the top quark. The treatment of the various background contributions to ttj 
production was a particular challenge, and a strategy was set up which strongly reduced the 
dependence on Monte Carlo simulations. 

The measurement was performed in the single electron and single muon channels and was 
hence subject to background contributions from W and Z boson production with additional jets 
(^+jets and Z-|-jets), single top, multijet and diboson production, where all of these processes 
featured an additional photon in the final state. The main background contributions, however, 
were due to processes where a hadron from jet fragmentation or an electron were misidentified 
as a photon. 

Since the misidentification of hadrons and electrons as photons is known to be poorly modelled 
in Monte Carlo simulations, a strategy which minimised the dependence on simulations was set 
up for the estimation of these main background contributions: the contribution from hadrons 
misidentified as photons was estimated from data with a template fit to the photon isolation 
distribution. The probability for electrons to be misidentified as photons was measured in 
Z ^ ee events in data. 

Data-driven techniques were developed for the estimation of the background processes from 
W^-|-jets and multijet production with an additional real photon in the final state. These con- 
tributions were estimated in control regions in data enhanced with the respective background 
process and then extrapolated to the signal region. The remaining background contributions 
were found to be small and were estimated using Monte Carlo simulations. 

A total of 122 ti^ candidate events were identified in the data set and a combined measurement 
in the single electron and single muon channel yielded a cross section of 

attj ■ BR = 1.9 ± 0.5 (stat.) ± 0.8 (syst.) ± 0.1 (lumi.) pb , 

where BR is the branching ratio into the single lepton and dilepton tt^ decay channels. This 
measurement is in agreement with the expectation from Standard Model calculations. The 
statistical significance of the tt^ signal was estimated to 2.5 a in agreement with the expected 
significance. 

The dominant systematic uncertainty was found to be due to the modelling of jets, in par- 
ticular due to the measurement of the jet energy scale. Other large systematic uncertainties 
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were found to be the modelling of the 6-tagging performance, the shape of the photon isolation 
templates, and the modelling of initial and final state radiation. 

Based on the techniques established in this thesis, a 3 a evidence for ttj production can be 
expected with the 5 fb~^ of data taken at ^ys = 7 TeV. If background contributions can further 
be reduced, a 5 a observation is in reach. Moreover, a measurement of the ratio of the ttj 
and tt production cross sections may reduce the impact of systematic uncertainties due to the 
detector modelling. With larger data sets, measurements of differential tt^ cross sections may 
enhance the sensitivity to the t7-vertex. tt'j analyses at the LHC in the next years will hence 
provide further milestones towards a measurement of the strength and the structure of the 
electromagnetic coupling of the top quark and will therefore contribute to our picture of the 
heaviest elementary particle known to date. 
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Figure A.l: pJj?'^'^^^ distributions for electrons from simulated Z — t- e~^e~ decays (upper part 
of each plot) in different bins of £^t for < \ri\ < 0.60 (four upper plots) and 
0.60 < \rj\ < 1.37 (two lower plots) normalised to unity. The lower part of each plot 
shows the difference of the distribution of photons from simulated tt'y events (solid 
line) with respect to the electron distribution. Additionally, the distributions for 
unconverted (dotted grey line) and converted photons (dashed black line) from 
tt'y simulations are depicted. In all plots, the last bin includes the overflow bin. 
Photons closer than 0.2 in r]-(p-spa,ce to a true lepton were not considered and the 
agreement with the electron distribution is improved with respect to the plots in 
Fig. 



8.3 and 8.4 These plots complete the example plots shown in Fig. 8.5 
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Figure A. 2: pjp'^'^^o (distributions for electrons from simulated Z — t- e~^e~ decays (upper part 
of each plot) in different bins of E't for 1.52 < \rj\ < 1.81 (four upper plots) and 
1.81 < \r]\ < 2.37 (four lower plots) normalised to unity. The lower part of each 
plot shows the difference of the distribution of photons from simulated tij events 
(solid line) with respect to the electron distribution. Additionally, the distribu- 
tions for unconverted (dotted grey line) and converted photons (dashed black line) 
from tt'y simulations are depicted. In all plots, the last bin includes the overflow 
bin. 

Photons closer than 0.2 in r]-(j)-space to a true lepton were not considered and the 
agreement with the electron distribution is improved with respect to the plots in 



Fig. 8.3 and 8.4 These plots complete the example plots shown in Fig. 



^.5 
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